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1.  ACCOMPLISHMENTS  DUE  TO  AFOSR  SUPPORT  : 


PRINCIPAL  INVESTIGATOR:  M.  R.  Flannery 
School  of  Physics,  Georgia  Institute  of  Technology 
Grant  AFOSR-89-0426,  Period  7/1/90  -  6/30/91 
Project  Task:  2301/A4 

1.1  RESEARCH  OBJECTIVES  AND  DEVELOPMENTS 

The  objective  of  the  present  research  program  is  to  formulate,  develop  and  implement  new  theoretical 
descriptions  of  the  following  atomic  and  molecular  processes  (A-C)  below. 


(A)  Laser- Assisted  Collisions:  A  new  theory  of  laser-assisted  electron- (excited)  atom  collisions, 

eT  4-  A  +  nhv  — ♦  4-  A*  4*  mhu  (1) 

in  which  the  dressed  states  of  the  atom  A  in  the  laser  field  are  closely  coupled  and  the  Volkov  states  of  the 
projectile  electron  in  the  laser  field  are  included  has  been  developed  by  M.  R.  Flannery  and  P.  H.  G.  Smith. 
One  paper  has  already  been  published,  a  second  is  in  process  of  publication  and  a  third  has  been  submitted 
for  publication.  See  Appendices  (A)-(C). 


(B)  Termolecular  Recombination:  The  transport-collisional  set  of  Master  equations  for  Teimolec- 
ular  Recombination, 


A4-B4-  M  -  AB  +  M  (2) 

as  a  function  of  gas  density  has  been  developed  by  M.  R.  Flannery.  See  Appendix  (D). 


(C)  Angular  Momentum  Changes  in  Collisions  with  excited  atoms:  Work  is  progressing  on 
the  cross  sections  for  angular  momentum  changes, 

A  +  B(n*)-A  +  B+4-«TM') 

in  heavy-particle  and  electron-atom  (e-B)  collisions  where  the  target  atom  is  initially  in  an  excited  state. 

The  cross  sections  for  ni  — »  €?  collisional  transitions  increase  as  t  is  increased  until  a  maximum  ^ax 
is  attained  after  which  the  cross  sections  decrease  preciptously.  M.  R.  Flannery  and  A.  Haffad  have  shown 
that  this  effect  not  only  can  be  explained  by  a  previous  quantum  description  of  Flannery  and  McCann,  but 
also  by  classical  scattering. 
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1.2  RESEARCH  COMPLETED  DURING  CURRENT  PERIOD  (7/1/90  -  6/30/91) 


Within  this  yearly  period,  the  following  projects  were  investigated,  completed  and  written  up  for  pub¬ 
lication: 

A.  Ph.  D.  Thesis:  *A  Semiclassical  Treatment  of  Laser  Assisted  Collisions  in  a  Soft-Photon  Weak-Field 
Regime’,  by  P.  H.  G.  Smith,  Georgia  Institute  of  Technology  (Ph.D.  awarded  6/3/91) 

B.  Papers  in  Press  and  submitted  for  publication: 

1.  ‘Electron- Atom  Collisions  in  a  Laser  Field’,  by  P.  H.  G.  Smith  and  M.  R.  Flannery,  Nucl.  instr. 
Meths.  Phys.  Res.  B  56/57  (1991)  166-9,  Appendix  A. 

2.  ‘ELectron-Hydrogen  Collisions  in  a  Laser  Field’,  by  P.  H.  G.  Smith  and  M.  R.  Flannery,  J.  Phys.  B: 
At.  Mol.  Opt.  Phys.  1991  (zn  press),  Appendix  B. 

3.  ‘Electron-Hydrogen  Collisions  with  Dressed  Target  and  Volkov  Projectile  States  in  a  Laser  Field’, 
by  P.  H.  G.  Smith  and  M.  R.  Flannery,  J.  Phy.  B:  At.  Mol.  Opt.  Phys.  1991  (submitted  for  publication ), 
Appendix  C. 

4.  ‘Transport-Collisional  Master  Equations  for  Termolecular  Recombination  as  a  function  of  Gas  Den¬ 
sity’,  by  M.  R.  Flannery,  J.  Chem.  Phys.  1991  (in  press  Oct.  issue),  Appendix  D. 

1.8  RESEARCH  COMPLETED  DURING  PREVIOUS  PERIOD  (7/1/89  -  6/30/90)  AND 
WRITTEN  UP  AND  PUBLISHED  DURING  CURRENT  PERIOD 
(7/1/90  -  6/30/91)* 

1.  ‘Recombination  Processes’,  M.  R.  Flannery  in  Molecular  Processes  in  Space,  ‘Physics  of  Atoms  and 
Molecules’  series,  edited  by  T.  Watanabe  ,  I.  Shimamura,  M.  Shimizu  and  Y.  Itikawa  (Plenum  Press,  London, 
1990)  Chapter  7, 

2.  ‘Electron  Collision  Cross  Sections  Involving  Excited  States’  E.  J.  Mar  sky  in  Nonequilibrium  Processes 
in  Partially  Ionized  Gases ,  NATO  ASI  series  B:  Physics  vol.  220,  edited  by  M.  Capitelli  and  J.  N.  Baidsley, 
Plenum  Press  1990,  pages  349-55. 

3.  ‘The  Issue  of  Basis  Set  Size  in  +  If(l«  — ♦  2s,  2p)  Collisions’,  E.  J.  Mansky  and  M.  R.  Flannery,  J. 
Phys.  B:  At.  Mol.  Opt.  Phys.  23  L501-507  (1990). 

4.  ‘Polarization  Ftactions  for  the  21P,  3lP  and  3 lD  states  of  Helium’  E.  J.  Mansky  and  M.  R.  Flannery, 
J.  Phys.  B:  At.  Mol.  Opt,  Phys.  23  3987-02  (1990). 

5.  ‘The  Multichannel  Eikonal  Theory  of  Electron-Hydrogen  Collisions  I.  Excitation  of  H(ls)’  E.  J. 
Mansky  and  M.  R.  Flannery  3.  Phys.  B:  At.  Mol.  Opt.  Phys.  23  4549-72  (1990). 

6.  ‘The  Multichannel  Eikonal  Theory  of  Electron-Helium  Collisions  I.  Excitation  of  He(l15)’  E.  J. 
Mansky  and  M.  R.  Flannery,  J.  Phys.  B:  At.  Mol.  Opt.  Phys.  23  4573-4604  (1990). 

7.  ‘Electron-Metastable  Helium  Differential  and  Integral  Cross  Sections’  E.  J.  Mansky  and  M.  R. 
Flannery,  J.  Phys.  B:  At.  Mol.  Opt.  Phys.  in  press 

8.  ‘Indirect  Coupling  Mechanisms  and  Stokes  Parameters  foi  Electron- Atom  Scattering’  E.  J.  Mansky 
and  M.  R.  Flannery,  J,  Phys.  B:  At.  Mol.  Opt.  Phys.  in  press 

Above  publications  nos.  1  and  5-8  were  included  as  Appendices  B-F  of  the  previous  Annual  Technical 
Report  GIT-89-001  for  the  period  7/1/89  -  6/30/90  .  Six  (6)  reprints  of  the  first  six  papers  (section  1.3 
nos.  1-6  above)  are  enclosed  separtely  with  this  annual  report  to  AFOSR  under  report  numbers  GIT-89-002, 
-003,  -004,  -005,  -006  and  -007,  respectively.  Reprints  of  the  remaining  papers  nos.  7  and  8  will  be  sent  to 
AFOSR  when  available. 
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1.4  SUMMARY:  PAPERS  PUBLISHED  AND  IN  PRESS 


A  total  of  twelve  (12)  papers  have  either  been  already  published  (as  detailed  in  sections  1.2  no.  1  and 
R3  nos.  1-6)  or  are  currently  in  press  (as  in  section  1.3  nos.  7,8  and  in  section  1.2  nos.  2,3  and  4),  or  have 
been  submitted  for  publication  (section  1.2  no.  3)  during  the  two  years  (7/1/89  -  6/30/91)  of  the  current 
AFOSR  Grant.  Reprints  of  all  of  the  above  papers  will  be  sent  to  AFOSR. 

In  addition,  two  Ph.D.  thesis  : 

1.  Termolecular  Ion- Atom  Association  of  Rare  Gase  Ions  in  Rare  Gases  by  M.  S.  Keenan 
(Ph.D.  awarded  3/17/90) 

and, 

2.  A  Seraiclassical  Treatment  of  Laser  Assisted  Collisions  in  a  Soft-Photon  Weak  Field 
Regime  by  P.  H.  G.  Smith  (Ph.D.  awarded  6/3/91) 

have  been  accomplished  due  to  AFOSR  support.  Copies  of  these  theses  will  also  be  sent  to  AFOSR. 
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2.  PAPERS  PRESENTED  AT  SCIENTIFIC  MEETINGS  (7/1/90  -  6/30/91) 


1.  ‘Electron- Atom  Collisions  in  a  Laser  Field’,  P.  H.  G.  Smith  and  M.  R.  Flannery,  Bull.  Amer.  Phys. 
Soc.  36  No.  2  (1991)  188 

2.  ‘Angular  Momentum  Changes  in  Collisional  Ioniza.ion’,  A.  Haffad  and  M.  R.  Flannery,  Bull.  Amer. 
Phys.  Soc.  36  No.  2  (1991)  188 

The  above  two  papers  were  presented  at  the  43rd  Annual  Gaseous  Electronics  Conference,  16-19  October, 
1990,  Champaign-Urbana,  Illinois. 

3.  ‘Functional  Parallelism  and  Atomic  Scattering  Theory*  E.  J.  Mansky,  to  appear  in  the  Proceedings 
of  the  Fifth  SIAM  Conference  on  Parallel  Processing  for  Scientific  Computing  (SIAM  Press  1992),  25-27 
March  1991,  Houston  Texas. 

2*1  ABSTRACTS  OF  PAPERS  PRESENTED 


Abstract  of  Contributed  Poster  Paper  presented  at  the  Fifth  SIAM 
Conference  on  Parallel  Processing  for  Scientific  Computing 
March  25-27,  1991,  Houston  ,  Texas 


Functional  Parallelism  and  Atomic  Scattering  Theory 

The  parallelizability  of  the  numerical  solution  of  systems  of  N  cou¬ 
pled  first-order  linear  partial  differential  equations,  which  arise  in  the 
solution  of  Schrodinger’s  equation  in  electron-atom  scattering,  is  in¬ 
vestigated.  In  particular,  an  optimal  strategy  is  outlined  for  paralleliz¬ 
ing  the  solution  of  systems  of  coupled  l3*-order  PDE’s  by  balancing 
the  competing  demands  of  scheduling  (ie.  load  balancing),  granularity 
and  computational  intensity  of  the  algorithm.  In  this  regard,  the  al¬ 
gorithmic  phase  diagram  of  Hockney  proves  instrumental  in  choosing 
which  type  of  numerical  technique  (ie.  rational  extrapolation,  Runge- 
Kutta,  predictor-corrector)  is  "best”  depending  on  the  number  of  cou¬ 
pled  equations  N  and  global  error  tolerance  chosen.  Algorithmic  phase 
diagrams  and  Hockney  numbers  (nx/2,  ^i/2j /1/2)  will  be  preserted  for 
the  solution  of  N  coupled  PDE’s  (N  =  20  -  100)  which  arise  in  the 
semi-  classical  multichannel  eikonal  theory  of  inelastic  electron-atom 
scattering.  Work  supported  bv  AFOSR  under  grant  no.  AFOSR-89- 
0426. 

E.  J.  Mansky 
School  of  Physics 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332-0430. 


6 


Papers  Presented  at  the  43rd  Annual  Gaseous  Electronics  Conference, 
Urbana-Champaign,  Illinois,  16-19  October,  1990 


D-7  Electron* Atom  Collisions  in  a  Laser  Field.*  Philip. H.G 
Smith  and  M.R. Flannery,  Georgia  Institute  cl Technc'.'^.  A  semi- 
classical  Floquet  approach  is  used  to  solve  exactly,  the  Schrodinger  equa¬ 
tion  for  the  laser/hydrogen  interaction  in  a  soft  photon  weak-field  limit, 
to  give  dressed  states  of  the  atom  in  the  laser  field.  Perturbative  dress¬ 
ing  is  shown  to  provide  an  incomplete  description,  and  cannot  predict 
the  distinctive  features  of  the  Floquet  approach.  Electron-hydrogen 
collisions  in  a  laser  field  are  then  described  via  a  multichannel  eikonal 
treatment,  in  which  the  dressed  states  are  closely  coupled.  Cross  sec¬ 
tions  for  1S-2S  and  lS-2Po  excitations  are  presented  as  a  function  of 
field  strength  and  impact  energy,  and  compared  with  the  Born-wave 
result. 

*  Research  supported  by  AFOSR-89-0426. 


D-8  Angular  Momentum  Changes  in  Collisional  Ionization* , 

A.  Haffad  and  M.  R.  Flannery,  Georgia  Institue  of  Technology  - 
Single  and  Double  Differential  cross  sections  for  ionization  in  e  -  H(nl) 
and  H(ls)  -  H(nl)  collisions  are  reported  as  a  function  of  impact  energy 
E,  final  energy  <  and  angular  momentum  Y  of  the  ejected  electron. 

This  process  is  assumed  to  occur  via  an  energy-changing  and 
angular  momentum-changing  binary  collision  between  the  Rydberg 
electron  in  state  nl  and  the  projectile  e  or  H(ls). 

The  atomic  projectile  can  also  be  excited.  Systematic  trends  in  the 
variation  of  the  classical  cross  sections  with  final  angular  momentum  1# 
are  discussed  and  are  in  accord  with  a  previous  quanta!  treatment1. 


‘Research  supported  by  U.S.  Air  Force  Office  of  Scientific  Reasearch 
under  Grant  No.  AFOSR-89-0426. 

lM.  R  Flannery  and  K.  McCann,  Phys.  Rev.  19  (1979)2206. 


3.  PERSONNEL  INVOLVED 


1.  Professor  M.  R.  Flannery  -  Principal  Investigator 

2.  Dx.  E.  3.  Mansky  -  Research  Scientist  II 

3.  Mr.  P.  H.  G,  Smith  -  Graduate  Student  (Ph.D.  completed  6/3/91) 

4.  Mr.  A.  Haffad  -  Graduate  Student  (Ph.D.  completed  8/16/91) 

5.  Mr.  X.  Qi  -  Graduate  Student 


4.  SPECIAL  HIGHLIGHT:  LASER  ASSISTED  ELECTRON-ATOM  COLLISIONS 


A  new  theory  of  laser  assisted  electron-atom  collisions  has  been  formulated  ,  developed  and  applied  to, 

e“  +  H(ls)  +  ATAu?  if(25,2po,  2p±i)  -f  N'hut 

The  laser  can  perturb  both  the  bound  atomic  electrons  as  well  as  the  incident  projectile  electron.  The  first 
effect  is  acknowledged  by  a  semiclassical  Floquet  approach  used  to  dress  the  excited  states  of  hydrogen  in 
the  laser  field.  This  approach  is  compared  to  dressing  by  the  traditional  Perturbative  approach  which  is  then 
shown  to  provide  an  imcomplete  description  of  the  laser  interaction  and  which  cannot  predict  the  distinct 
features  provided  by  the  Floquet  approach. 

The  second  effect  of  the  laser  interaction  on  the  projectile  electron  is  acknowledged  via  Volkov  dressed  states 
for  the  projectile.  These  states  are  shown  to  exert  significant  influence  on  the  cross  sections  for  individual 
state-to-state  transitions  which  involve  absorption  or  emission  of  a  specified  number  of  photons,  They 
however  have  only  a  negligible  effect  on  the  cross  sections  obtained  by  summing  over  all  absorptions  and 
emissions. 
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APPENDIX  A 


Electron^ Atom  Collisions  in  a  Laser  Field 
P.  H.  G.  Smith  and  M.  R.  Flannery 
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Electron-atom  collisions  in  a  laser  field 

Philip  H.G.  Smith  and  M.R.  Flannery 

St  faun  ft  Pin  \ti  r.  i/eor&a  Institute  of  Tethtwln  Atlanta.  (/A  .W.O-'J-W.  iSA 


Cross  sections  for  the  1S-2S  and  1S-2P,,  transitions  in  laser  assisted  e-HMS)  collisions  are  calculated  in  both  the  multichannel 
eikonal  and  the  Born-wave  treatments  as  a  function  of  impact  energy  and  laser  field  intensity  and  phase.  The  laser  considered  i>  a 
monotonic,  plane  polarized  CO:  laser  (photon  energy  =*0.117  eV),  with  the  polarization  direction  parallel  to  the  initial  projectile 
\elocttv  Floquet  dressing  of  the  hvdrogen  atom  in  the  soft-photon  weak-field  limit  reveals  a  concise  description  of  the  laser  assisted 
electron-atom  collision.  This  model  also  links  the  microscopic  detail  of  the  individual  collisions  with  the  macroscopic  considerations 
of  experimental  analysis. 


1.  Introduction 

The  work  reported  in  this  paper  is  a  study  of  the 
influence  of  the  laser  field  on  electron-atom  collisions 
in  the  soft-photon  weak-field  regime.  In  this  regime  the 
photon  energy  is  a  lot  less  than  the  energy  required  to 
ionize  the  atom,  and  the  field  strengths  can  always  be 
considered  as  a  perturbation  to  the  field  of  the  nucleus 
on  the  bound  electrons.  The  effect  of  coupling  a  laser 
field  to  a  projectile  electron  in  this  regime  during  a 
collision  with  an  atom  has  been  well  explained  in  a 
number  of  studies  (1).  The  effect  of  a  laser  field  on  the 
target  atom,  however,  has  met  with  a  lot  less  success. 
This  is  due  to  the  off-diagonal  elements  introduced  into 
the  Schrddinger  equation  for  an  atom  in  a  laser  field, 
which  not  only  provide  couplings  between  eigenstates  of 
the  isolated  atom,  but  simultaneously  involve  the  ab¬ 
sorption  or  emission  of  a  photon.  In  this  model  the 
off-diagonal  elements  are  dealt  with  by  creating  dressed 
atomic  states  for  the  atom  tn  the  laser  field  by  a 
semiclassical  Floquet  approach  [2].  These  dressed  states 
are  then  used  to  solve  the  electron- hydrogen  collision 

e'+  H(1S)  +  .VAw  -e"  +  H(2S,  2P0,  2P.,)  +  Mu. 

U) 

This  model  can  include  the  laser-projectile  interac¬ 
tion  by  using  the  well-known  Volkov  states.  As  is  shown 
in  ref.  [3],  the  Volkov  states  for  the  range  of  impact 
energies  considered  in  this  work  -  50  eV.  100  eV  and 
200  eV  -  provide  a  wide  range  of  high  order  multipho¬ 
ton  transitions.  But  as  is  also  shown  in  ref.  (3],  when  the 
state-to-state  cross  sections  are  summed  over  an  increas¬ 
ingly  wider  range  of  projectile  energies,  the  full  treat¬ 
ment  which  includes  Volkov  dressing  is  found  to  con¬ 
verge  on  the  treatment  that  neglects  the  laser- projectile 
interaction.  This  finding  has  been  used  to  justify  pre¬ 


senting  state-to-state  cross  seettons  in  this  paper,  calcu¬ 
lated  when  neglecting  the  laser-projectile  interaction. 


2.  The  laser-atom  coupling 

The  starting  atomic  Hamiltonian  is  (4) 

H,  d(r.  t)  -  +  yjr)  +  £0 •  r  sm( ur  +  SJ  (2) 

where  H0  is  the  field-free  atomic  Hamiltonian  and  £0  is 
the  electric  field  strength  *  A0v/c  expressed  in  terms  of 
the  laser  frequency  u.  The  phase  shift  is  which  is 
explained  in  section  4  of  this  paper.  This  Hamiltonian 
can  be  seen  to  be  time  dependent.  By  expanding  into 
Floquet  states  and  diagonalizing  a  restricted  version  of 
the  resultant  infinite  matrix  [3,5],  time  independent 
dressed  states  can  be  obtained  that  exactly  solve  eq.  (2). 
The  time  dependence  has  been  removed  by  allowing  for 
absorption  and  emission  of  photons  of  energy,  so  that 
the  new  dressed  atomic  states  can  be  used  to  solve 
collisional  problems  in  a  time  independent  manner  di¬ 
rectly  analogous  to  field-free  atomic  states.  The  cost  of 
removing  the  time  dependence  is  an  infinite  set  of 
atomic  slates  and  energy  levels 

+  (3a) 

I  n+m).  (3b) 

a  tn 

where  £a  is  a  summation  over  field-free  atomic  ‘nates 
<pa%  and  is  a  summation  over  their  Floquet  expan¬ 
sion.  It  is  emphasized  that  n  vanes  across  the  range 
-  x  to  x,  and  so  provides  an  infinite  set  of  periodic 
dressed  states. 
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3.e'+H(lS)+.VAw  collision 

Now  that  dressed  atomic  states  have  been  attained, 
they  are  used  to  solve  the  time  dependent  laser  per¬ 
turbed  collisionai  Hamiltonian 

A, 

<  £  Jm  exp(i/n<$_  ).  (4) 

’!)  *  -  TC 

in  a  time  independent  manner.  The  summation  over 
Bessel  functions  appears  from  the  Volkov  dressing  of 
the  projectile  electron.  As  mentioned  earlier,  this  has 
been  well  explained  by  others  and  so  will  not  be  dealt 
with  m  detail  in  this  paper.  Any  experiment  detects  the 
products  of  the  collision  outside  the  region  of  the  laser 
field  where  //a  d  =  //„.  To  obtain  results  comparable 
with  experimental  observations  this  model  calculates 
the  state-to-state  cross  sections  from  the  “measured” 
transition  amplitude  taken  with  respect  to  the 
field-free  atomic  states.  Eq.  (4)  provides  a  solution  to 
the  “dressed”  transition  amplitude  Tpq  taken  with  re¬ 
spect  to  the  Floquet  dressed  atomic  states.  By  projecting 
the  dressed  states  Fp  on  to  the  field  free  states  <t>p  it  is 


possible  to  write  the  “measured”  transition  amplitude 
in  terms  of  the  “dressed”  transition  amplitude 

=  EII[  °n . r,  r. : o„ •  r„ ( ./„, e\p( i ) j . 

r  o  ”» 

<5) 

The  "measured”  transition  amplitude  can  he  Jiuded 
into  two  distinct  parts  (l)  <t\  -t.  Aw)  F  J.„  de¬ 
scribes  the  collisionai  transition  between  dressed  states 
of  the  atom:  (2)  (  4.  Aw)  or  the  projection  o,  /*. 
describes  the  laser-atom  interaction  betore  and  after 
the  collisionai  event  ( e.  4.  Aw)  There  is  no  analogous 
laser-projectile  interaction  ( e ,  Aw)  in  this  model  Both 
these  parts  allow  for  absorption  and  emission  of  pho¬ 
tons  through  interactions  with  the  laser  field.  This  model 
can  be  seen  to  give  a  very  concise  description  of  the 
laser  assisted  electron-atom  collision,  with  all  three 
interactions  (e%  A ,  Aw).  {,4.  Aw)  and  (e.  Aw)  as 
described  recently  [6].  included  exactly. 


4.  The  dipole  approximation 

This  model  is  based  upon  a  typical  crossed  beam 
experiment,  where  a  projectile  electron  beam,  a  target 


Fig.  1.  Cross  section  vs  laser  field  strength,  for  third  order  Floquet  dressed  hydrogen  atoms,  in  laser  assisted  e  +H(1S)  collision*, 
using  a  Bom-wave  treatment,  (a)  1S-2S  and  (b)  1S-2P0  cross  sections  Cross  sections  are  compared  over  a  range  of  phases  and 

impact  energies. 
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Fig.  3  IS-2S  and  1S-2P0  cross  sections  vs  laser  field  strength,  for  third  order  Floquet  dressed  hydrogen  atoms,  in  laser  assisted 
e  +H(1S)  collisions,  using  (a)  Bom-wave  and  (b)  multichannel  Eikonal  treatments.  The  cross  sections  are  averaged  over  in  the 

range  0  °  to  90  °  and  compared  over  a  range  of  impact  energies. 
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atomic  beam  and  a  laser  beam  all  cross  at  one  point  in 
space.  with  the  vector  potential  for  the  laser  field  ex¬ 
pressed  as 

A  =  € . /t0  COS(  &  T-raj/  +  Su)).  (6) 

At  the  atomic  level,  a  dipole  approximation  can  be  used 
in  all  calculations  so  that  eq.  (6)  can  be  reduced  to 

A  =c.A0  cos(wr  +  0.  (7) 

However,  at  the  macroscopic  level,  the  dipole  approxi¬ 
mation  does  not  hold  across  the  full  width  of  either  the 
projectile  or  the  atomic  beams.  This  can  be  dealt  with 
by  varying  the  phase  factor  8^  in  A.  The  wavelength  for 
the  electromagnetic  radiation  considered  in  this  work  is 
of  the  order  of  10" 5  m.  Hence  any  experiment  will 
observe  a  range  of  collisions  over  many  wavelengths 
and  time  periods  of  the  laser.  To  predict  the  experimen¬ 
tal  cross  sections  it  is  thus  necessary’  to  take  an  average 
over  a  range  of  from  0°  to  360°. 

^experiment  =  JZ  f <K$U)-  (8) 


5.  Results 

As  explained  in  the  introduction,  the  Volkov  dress¬ 
ing  of  the  projectile  states  provides  little  additional 
information  to  the  summed  $tate«to-state  transitions. 
Hence  the  results  presented  here  neglect  the  laser-pro¬ 
jectile  interaction,  and  are  still  considered  to  be  a  good 
approximation.  All  Floquet  dressed  states  are  taken  to 
third  order  in  approximation,  which  allows  for  multi¬ 
photon  couplings  of  up  to  three  photons  in  height. 
Third  order  is  used  since  it  is  shown  that  convergence  in 
the  cross  sections  is  reached  by  this  time  (5). 

Numerical  calculations  of  the  dressed  transition  am¬ 
plitudes  of  eq.  (5)  are  based  on  a  Bom-wave  treatment 
and  a  multichannel  Eikonal  treatment  (7).  The  restricted 
basis  set  -  IS,  2S.  2P0  and  2P±l  -  is  used  in  the 
dressing  of  the  atomic  states,  which  in  turn  provides  the 
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basis  set  used  in  the  multichannel  eikonal  approxima¬ 
tion.  The  axis  of  quantization  is  along  the  direction  of 
incidence  of  the  projectile  electron,  and  is  also  taken  as 
the  direction  of  polarization  of  the  laser  field. 

Fig.  1  shows  the  1S-2S  and  IS-2R,  cross  sections 
calculated  by  the  Born-wave  treatment  for  a  ranae  of  8 
from  0°  to  40°  The  same  set  of  curves  is  repeated  over 
the  next  270°  Hence  the  predicted  experimental  cross 
section  need  be  averaged  only  over  the  reduced  ran^e 
6.-0°  to  90°.  Fig.  2  shows  the  same  1S-2S  and 
IS-2P0  state-to-state  cross  sections  calculated  bv  the 
Eikonal  treatment.  As  can  be  seen,  the  two  treatments 
give  very  similar  results.  Fig.  3  shows  the  averaged 
1S-2S  and  1S-2R}  cross  sections,  (a)  calculated  m  the 
Bom-wave  approximation,  and  (b)  calculated  in  the 
multichannel  eikonal  approximation. 
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Georgia  Institute  of  Technology,  Atlanta,  GA,  30332-0430,  U.S.A. 

ABSTRACT:  The  non-perturbative  Floquet  method  is  used  to  provide  the  dressed  states 
of  a  hydrogen  atom  in  a  laser  field  in  the  soft-photon  weak-field  regime.  These  dressed 
atomic  states  then  provide  a  basis  set  expansion  for  use  within  a  consistent  semiclassical 
Multichannel  Eikonal  Treatment  of  laser  assisted  e~-H(lS)  collisions.  The  variations  with 
field  strength  of  the  1S-2S  and  lS-2Po  state-to-state  cross  sections  are  presented.  Special 
attention  is  employed  in  correlating  the  time  frame  of  the  laser  field  with  the  time  frame 
of  the  relative  orbit  of  the  collisional  species,  and  this  is  shown  to  require  the  inclusion  of 
a  phase  shift  6 w  within  the  vector  potential  of  the  laser  field.  This  inclusion  is  important 
when  comparing  with  experimental  results. 

1. INTRODUCTION 

The  Floquet  treatment  has  already  been  successfully  employed  (Chu  1985,  Potvliege 
and  Shakeshaft  1991)  in  calculations  of  laser  induced  multiphoton  ionizations,  where  it 
provides  dressed  states  for  an  atom  in  a  laser  field.  That  (perturbative)  dressing  of 
target  states  can  have  important  consequences  in  laser-assisted  scattering  was  illustrated 
by  Byron  and  Joachain  (1984).,  These  dressed  states  are  useful,  not  only  for  laser  induced 
phenomenon,  but  also  as  a  collisional  basis  set  for  laser  assisted  collisions.  In  this  role  they 
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are  in  fact  very  appealing,  since  the  Floquet  treatment  naturally  lends  itself  to  a  time- 
independent  analysis,  and  hence  are  compatible  with  present  field  free  scattering  theories. 
Despite  the  apparent  applicability  of  this  approach,  work  along  these  lines  has  only  just 
recently  appeared  (Sharmaand  Mohan  1990,  Smith  and  Flannery  1991a,  Burke  et  al  1991). 
Byron  and  Joachain  (1984)  have  illustrated  that  perturbative  dressing  of  the  target  states 
can  have  important  consequences  in  laser-assisted  scattering.  Floquet  dressing  however 
provides  a  more  complete  description  (Smith  and  Flannery  1991b). 

This  short  paper  provides  an  outline  of  the  research  (Smith  1991)  being  conducted 
by  the  authors  on  e~-H(lS)  collisions,  using  a  Floquet  dressed  basis  set  for  a  hydrogen 
atom  in  a  laser  field.  A  short  paper  has  previously  been  published  (Smith  and  Flannery 
1991a)  but  the  range  of  the  electric  field  strength,  over  which  collisional  cross  sections 
are  calculated,  has  now  been  extended  to  provide  greater  insight  into  the  role  of  the  laser 
field  in  the  collision.  A  more  detailed  account  is  presently  under  preparation.  The  work 
described  here  centers  on  the  laser  perturbation  of  the  atom  alone,  and  attempts  to  probe 
its  exclusive  effect  by  neglecting  the  laser  perturbation  of  the  projectile  electron.  In  a  later 
paper,  the  laser  perturbation  of  the  projectile  will  also  be  included  via  the  use  of  Volkov 
dressed  states  but  these  will  be  shown  to  have  a  negligible  effect  on  the  state-to-state  cross 
sections,  within  the  range  of  electric  field  strengths  studied  in  this  work. 


2. THEORY  OF  LASER  ASSISTED  COLLISIONS 

The  starting  point  of  a  discussion  of  laser  assisted  collisions  is  the  laser  perturbed 


2 


Schrodinger  equation 


TvPr  +  £  5(P'i  +  A/c)2  +  v(r)  +  V'(R,r) 


(1) 


The  channel  coordinates,  (Pfl,R),  represent  the  momentum  operator  and  the  position 
vector  for  the  projectile.  The  internal  coordinates,  (p,  r)  represent  the  momentum  operator 
and  the  position  vector  for  each  bound  electron.  The  electronic  coordinates  will  be  denoted 
collectively  by  r,  c  is  the  speed  of  light,  n  is  the  projectile-target  reduced  mass,  V(r)  is  the 
internal  potential  for  the  bound  atomic  electrons  and  V(R,  r)  is  the  external  projectile- 
target  interaction  potential. 

After  the  usual  dipole  approximation,  the  vector  potential  for  a  monochromatic  plane 
polarized  laser  of  frequency  u>  is 


A  =  A0  cos(u >/  +  6W)  (2) 

It  is  customary,  especially  in  multiphoton  ionization,  to  omit  the  phase  shift  Sw  and  to  just 
write  the  vector  potential  in  the  form  A0  cos  u>f  or  A0sinu>t.  However  for  semiclassical 
collision  theory  there  is  already  a  time  frame  of  reference  imposed  upon  the  collision, 
namely  the  time  along  the  relative  trajectory  of  the  collisional  species.  The  time  i  =  0 
is  usually  defined  at  the  point  of  closest  approach  between  the  projectile  and  the  target, 
ie  the  orbit’s  passage  through  the  periapsis.  Thus,  if  the  vector  potential  were  used  in 
the  form  A0cosu>f,  it  would  indicate  for  each  e-atom  collision  that  the  laser  field  was  at 
a  maxima  of  its  cycle  at  the  point  of  closest  approach,  as  indicated  by  the  dotted  curve 
in  the  diagram  of  Fig.l.  This  is  in  general  not  true,  so  that  it  is  necessary  to  adopt  the 
phase  shift  8W  in  order  to  synchronize  the  time-zero  of  the  laser  field  with  the  time-zero  of 
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the  collisional  e-atom  orbit.  The  vector  potential  of  eq.(2),  used  in  these  calculations,  is 
therefore  shifted  in  phase  by  Sw,  from  A0cosu as  shown  by  the  solid  line  in  the  diagram 
of  Fig.l.  The  phase  shift  varies  between  0  and  2n.  Bachau  and  Shakeshaft  (1984)  have 
explicitly  acknowledged  the  effect  of  the  phase  shift  for  excitation  in  H+  —if (Is)  scattering 
in  a  nearly  resonant  laser  field. 

On  omitting  the  projectile  dependent  terms  from  the  laser  perturbed  Schrodinger 
equation  eq.(l),  the  time-independent  Floquet  prescription  (Shirley  1965,  Sambe  1973, 
Chu  1985)  yields  the  matrix  equation. 


E 


(^a  4"  ^hw)^3Q6n— ,7,  H*  ^n  —  m  —  1  H pa&n=m  +  l 


Ko  =  Q,K 


(3) 


The  field  free  atomic  energy  levels  eQ  and  the  laser  induced  off-diagonal  couplings 
fip0  =  —i~(4>^ |e0  •  r|</>a),  between  field  free  atomic  states,  are  known  so  that  the  matrix 
can  be  diagonalized  to  obtain  the  new  dressed  basis  set.  This  involves  determining  the 
“Quasi”  energies  Qq  and  the  dressed  atomic  states 


=  E  E  -c.  i*»>  i")  (4) 

rt  a 

which  are  written  in  terms  of  the  above  combination  of  field  free  atomic  states  |^»a)  and 
periodic  states  \n)  =  emwl .  These  new  dressed  atomic  states  form  an  orthogonal,  time- 
independent  basis  set,  which  provides  a  medium  in  which  any  time-independent  field  free 
scattering  theory  may  be  applied  to  laser  assisted  collisions.  However  the  price  paid  for  this 
formal  simplification  is  the  introduction  of  an  infinite  set  of  periodic  solutions,  demanded 
by  the  Floquet  prescription. 

Using  the  Lippman  Schwinger  equation,  it  is  relatively  straightforward  (Smith  1991) 
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to  obtain  the  “dressed”  transition  amplitude 


Tpq  =  (i i/>p\V(R,r)\4'q )  exp(i(k,-  —  k/)  -R)  6(QP  +  j^k)  +  (n-m)hw  -Qq-  (5) 

taken  with  respect  to  the  dressed  atomic  states  4'q-  This  “dressed”  transition  amplitude 
describes  collisional  transitions  between  the  dressed  atomic  states,  and  includes  all  of 
the  possible  photon  absorptions  and  emissions  that  can  occur  during  such  transitions. 
However  the  “dressed”  transition  amplitude,  eq.(5),  pertains  to  the  wrong  physical  basis 
set  for  comparison  with  field  free  collisions  and  experimental  results.  Rather  a  transition 
amplitude  taken  with  respect  to  the  field  free  atomic  states,  <j)a ,  is  required.  Since  the  set  of 
dressed  atomic  states,  obtained  via  the  Floquet  prescription,  is  assumed  to  be  a  complete 
and  normalized  set  it  is  possible  to  form  the  projection  operator  Y^q  IV’gXV’gl  =  T  On 
using  this  projection  operator,  the  “measured”  transition  amplitude 

Tga  =  <«MV(R,r)|*„)  »££  T„  (6) 

9  9 

taken  with  respect  to  the  held  free  atomic  states,  may  be  written  in  terms  of  the  “dressed” 
transition  amplitude  Tpg,  taken  with  respect  to  the  dressed  atomic  states.  Writing  the 
“measured”  transition  amplitude,  T^a,  in  terms  of  the  “dressed”  transition  amplitude,  Tvq , 
in  this  manner  effectively  allows  for  a  transition  between  the  dressed  and  the  field  free  basis 
sets.  Hence  semiclassical  scattering  theories  can  be  applied  in  a  time-independent  manner 
and  the  calculated  transition  amplitudes  then  transformed  to  the  field  free  “measured” 
basis  set  in  order  to  yield  the  state-to-state  transition  cross  sections.  It  can  be  shown 
that  the  calculated  cross  sections  allow  for  laser  photon  absorption  and  emission  and  also 
Raman  photon  emission  (Smith  and  Flannery  1991b).  Both  the  probability  amplitudes 
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and  the  photon  absorption  and  emission  locations  during  the  collisional  sequence  of  events 
are  tracked  by  eq.(6). 

Since  the  vector  potential  (2)  can  only  be  defined  to  within  a  phase  shift  6^,  as  in 
Fig.l,  it  is  therefore  necessary  to  take  the  following  average 

i  r2* 

^experiment  =  X  /  (7) 

^  Jo 

by  <r  over  the  range  of  8W  from  0°  to  360°.  This  proves  to  be  an  important  consideration. 

RESULTS:  e”  +  H(1S)  +  Nhw  — >  e"  +  H(2S,2P0)  +  MTuo 

State-to-state  cross  sections  for  the  1S-2S  and  lS-2Po  transitions,  in  laser  assisted 
e~-H(lS)  colhsions,  are  calculated  by  the  Multichannel  Eikonal  Treatment  <*s  a  function 
of  impact  energy  and  laser  field  intensity.  All  cross  sections  reflect  a  summation  over 
the  state-to-state  cross  sections,  for  all  possible  final  projectile  energies,  consistent  with  a 
specified  initial  relative  energy.  Changes  in  the  final  projectile  energy  arise  from  photon 
absorption  and  emission.  The  laser  considered  is  a  monotonic,  plane  polarized  CO2  laser 
(hu>  =  O.llTeV),  with  the  polarization  direction  parallel  to  the  initial  projectile  velocity. 
The  restricted  basis  set  -  IS,  2S,  2Po  and  2P±]  -  is  used  in  dressing  the  states,  which 
in  turn  provide  the  basis  set  used  within  the  Multichannel  Eikonal  Treatment.  The  axis 
of  quantization  is  along  the  direction  of  incidence  of  the  projectile  electron  and  is  also 
taken  as  the  direction  of  polarization  of  the  laser  field.  Convergence  in  the  cross  sections  is 
reached  (Smith  and  Flannery  1991b)  by  a  third  order  approximation  to  the  dressed  atomic 
states  and  justifies  the  use  of  a  third  order  approximation  here. 

Cross  sections  versus  field  strength  for  1S-2S  and  1S-2P<>  transitions  are  presented  in 
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Fig-2,  (a)  and  (b)  respectively,  for  Floquet  dressed  hydrogen  atoms.  These  cross  sections 
are  compared  over  a  range  of  phase  shifts  from  Su  =  0°  to  90°,  since  cross  sections  for 
6^  between  90°  and  360°  are  found  to  be  a  repeat  of  the  curves  shown.  It  can  be  seen 
that  the  cross  sections  exhibit  a  very  dramatic  phase  dependence  over  the  field  strength 
shown.  The  Sw  =  0°  curve  exhibits  two  stationary  points,  while  the  Su  =  90°  curve 
exhibits  only  one  stationary  point.  A  model  purporting  predictions  of  experimental  results 
must  take,  as  explained  earlier,  an  average  of  the  cross  sections  over  the  range  of  phases 
from  =  0°  to  360°,  or  the  reduced  range  of  8W  =  0°  to  90°.  Cross  sections  after  such 
averaging  are  presented  in  Fig. 3a  where  it  is  immediately  apparent  that  the  1S-2S  cross 
section  never  crosses  the  lS-2Po  cross  section.  This  non-crossing  is  very  striking  when 
compared  to  the  percentage  2Po  component  of  the  dressed  state  S  (dressed  2S  state)  and 
the  dressed  state  P  (dressed  2Po  state)  presented  in  Fig.3b.  It  is  evident  that  the  dressed 
2S  state  is  gaining  an  increasingly  2Po  character  as  the  electric  field  strength  increases, 
thereby  raising  the  1 S-2S  cross  section  at  the  expense  of  the  lS-2Po  cross  section.  But 
because  of  the  infinite  set  of  periodic  Floquet  states  the  IS  to  dressed  2S  cross  section  is 
increasingly  contributing  to  the  “measured”  lS-2Po  cross  section,  and  hence  prevents  the 
“measured”  1S-2S  cross  section  from  exceeding  the  “measured”  lS-2Po  cross  section.  The 
reasons  behind  this  phenomenon  are  explained  in  greater  depth  in  a  later  paper  (Smith 
and  Flannery  1991b). 

4. CONCLUSION 

The  Floquet  prescription  used  to  obtain  the  dressed  collisional  basis  set,  has  been 
found  to  yield  a  very  concise  description  of  the  role  of  a  laser  field  in  collisional  processes. 
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The  overall  semiclassical  approach  is  consistent  in  that  both  the  laser  and  the  electron 
perturbations  are  described  semiclassically.  This  description  includes  both  the  probability 
amplitudes  and  locations  of  photon  absorption  and  emission  during  the  collisional  sequence 
of  events,  and  are  summarized  by  the  descriptions  (A ,hu>)  and  (e,A,hu>)  as  given  by  Newell 
(1990). 

This  work  has  demonstrated  the  significant  dependence,  of  the  collisional  cross 
sections,  on  the  phase  shift  6#  included  within  the  vector  potential  for  the  laser  field. 
This  phase  has  been  found  necessary  to  synchronize  the  time-zero  of  the  laser  field  with 
the  time-zero  of  the  collisional  orbit.  It  is  suggested  that  any  theoretical  model  which 
attempts  to  predict  experimental  results,  must  take  an  average  of  the  phase  dependent 
cross  sections  over  a  range  of  phases  from  6 w  =  0°  to  360°.  With  this  averaging  it  is 
shown,  in  Fig.3a,  that  the  1S-2S  cross  section  rises  up  to,  but  does  not  exceed,  the  lS-2Po 
cross  section  in  laser  assisted  e~-H(lS)  collisions.  This  feature  should  be  experimentally 
observable  with  present  day  technologies. 
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FIGURE  CAPTIONS 


Figure  (1):  Vector  potentials  _  A0cos(a>f  +  Su)  and  .  A0  cos  u>t,  for  a 

monotonic  plane  polarized  laser  field  of  frequency  u>,  are  shown  positioned  in  time  relative 
to  the  orbit’s  pericenter  at  time  t  =  0. 

Figure  (2):  Cross  sections  vs  field  strength  for  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e~-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  and  a  Multichannel  Eikonal  Treatment.  The  cross  sections  are  presented 

over  a  range  of  phase  shifts  f>w _ =0°, _ =22.5°,  . =45°, _ 

=67.5°  and - =90°,  at  a  photon  energy  of  hw  =  0.0043au. 

Figure  (3).  Cross  sections  in  e--H(lS)  collisions  are  compared  against  the  percentage  2Po 
component  of  the  dressed  atomic  states,  (a)  and  (b)  respectively,  over  a  range  of  electric 
field  strengths  for  third  order  Floquet  dressed  hydrogen  atoms  in  a  laser  field.  The  cross 

sections  ( -  1S-2S  and _ lS-2Po)  are  an  average  of  the  phase  dependent 

cross  sections  over  a  range  of  phases  from  =  0°  to  90°,  using  a  Multichannel  Eikonal 

Treatment.  For  the  percentage  2Po  component, _  represents  the  dressed  state  S 

and _ represents  the  dressed  state  P. 
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ELECTRON-HYDROGEN  COLLISIONS  IN  A  LASER  FIELD 


Philip. H.G. Smith  and  M.R.Flannery 
School  of  Physics, 

Georgia  Institute  of  Technology,  Atlanta,  GA,  30332-0430,  USA. 


Abstract:  Cross  sections  for  the  1S-2S  and  1S-2P0  transitions  in  laser  assisted  e~-H(lS) 
collisions  are  calculated  in  both  the  Multichannel  Eikonal  Treatment  and  the  Born- Wave 
approximation,  as  a  function  of  impact  energy  and  laser  field  intensity.  The  laser  considered 
is  a  monotonic,  plane  polarized  CO2  laser  (photon  energy  =  0.117eV)  with  the  polarization 
direction  parallel  to  the  initial  projectile  velocity.  The  first  part  of  this  paper  confines  the 
laser  perturbation  to  the  bound  electrons  of  the  atom.  A  semiclassical  Floquet  approach 
is  used  to  dress  the  hydrogen  atom  in  this  soft-photon  weak-field  regime,  and  is  shown 
to  reveal  a  concise  description  of  the  laser  assisted  collision.  The  Floquet  dressing  is 
compared  to  dressing  by  the  traditional  time-dependent  perturbation  theory,  showing  that 
the  Perturbative  approach  gives  an  incomplete  description  of  the  laser  interaction,  and 
cannot  predict  the  distinct  features  provided  by  the  Floquet  approach.  The  second  part 
of  this  paper  extends  the  laser  perturbation  to  the  projectile  electron,  and  the  familiar 
Volkov  dressed  states  are  used.  Although,  in  the  range  of  impact  energies  and  electric 
field  strengths  considered,  the  Volkov  dressed  states  exert  significant  influence  on  the  cross 
sections  for  individual  state-to-state  transitions,  which  involve  absorption  or  emission  of 
a  given  number  of  photons,  they  have  only  a  negligible  effect  on  the  cross  sections  when 
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summed  over  all  absorptions  and  emissions. 

Special  attention  is  employed  in  synchronizing  the  time  frame  of  the  laser  field  with 
the  time  frame  of  the  trajectory  of  the  collisional  species  orbit.  This  requires  the  inclusion 
of  a  phase  shift  <5W  within  the  vector  potential  of  the  laser  field.  This  inclusion  is  important 
when  comparing  theoretical  cross  sections  with  cross  section  measurements. 
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1. INTRODUCTION 


The  work  reported  in  this  paper  studies  the  influence  of  the  laser  field  on  projectile- 
atom  collisions,  in  the  soft-photon  weak-field  regime.  In  this  regime  the  photon  energy 
is  much  less  than  the  energy  difference  between  sub-levels  of  the  atom,  and  the  electric 
field  strength  of  the  laser  can  always  be  considered  as  a  perturbation  to  the  electrostatic 
interactions  between  the  nucleus  and  the  bound  atomic  electrons.  The  effect  of  coupling 
a  laser  field  to  a  projectile  electron,  in  this  regime,  during  a  collision  with  an  atom 
has  been  well  explained  (Volkov  1935,  Kroll  and  Watson  1973,  Geltman  and  Macquet 
1989)  via  the  introduction  of  Volkov  states.  Byron  and  Joachain  (1984)  have  illustrated 
that  perturbative  dressing  of  the  target  states  can  have  important  consequences  in  laser- 
assisted  scattering.  The  effect  of  a  laser  field  on  the  target  atom  however  involves  off- 
diagonal  elements  introduced  into  the  Schrodinger  equation  of  the  target  atom  by  the 
laser  field,  which  not  only  provide  couplings  between  eigen  states  of  the  isolated  atom,  but 
simultaneously  involve  the  absorption  or  emission  of  a  photon,  which  in  turn  demands  the 
use  of  a  non-perturbative  approach.  Little  work  has  been  done  on  this  aspect  although 
some  has  recently  emerged  (Smith  and  Flannery  1991,  Burke  et  al  1991). 

This  work  mainly  centers  on  the  laser  perturbation  of  the  atom,  and  initially  attempts 
to  probe  its  exclusive  effect  by  neglecting  the  laser  perturbation  of  the  projectile.  Before 
the  laser  perturbed  atomic  Schrodinger  equation  can  be  used  to  calculate  cross  sections  in  a 
laser  field,  it  is  essential  to  remove  the  time-dependence  and  to  re-instate  the  orthogonality 
of  the  atomic  states.  This  is  achieved  using  a  semiclassical  Floquet  approach  (Shirley  1965, 
Sambe  1973,  Chu  1985)  which  exactly  solves  the  laser  perturbed  Schrodinger  equation 
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to  produce  a  new  set  of  orthogonal  dressed  atomic  states,  which  can  then  be  treated 
in  a  time-independent  manner.  Once  these  new  states  have  been  determined,  they  can 
be  used  within  any  time-independent  scattering  theory  in  a  manner  directly  analogous 
to  field  free  collisions.  In  this  work  cross  sections  will  be  presented  in  the  Born- Wave 
approximation  and  the  Multichannel  Eikonal  Treatment  (Flannery  and  McCann  1975, 
Mansky  and  Flannery  1990)  and  are  compared  over  a  range  of  impact  energies  and 
electric  field  strengths.  It  will  be  shown  that  both  the  Born- Wave  approximation  and 
the  Multichannel  Eikonal  Treatment  produce  cross  sections  that  exhibit  the  same  essential 
dependence  on  the  electric  field  strength.  As  a  comparison  with  the  Floquet  approach, 
cross  sections  will  also  be  determined  via  the  more  traditional  time-dependent  Perturbation 
theory  (Bayfield  1979).  The  Floquet  and  the  Perturbative  cross  sections  are  strikingly 
different,  and  it  is  claimed  that  this  difference  originates  from  an  important  omission  in 
the  description  of  the  dressed  atomic  states  by  the  Perturbative  approach.  Floquet  dressing 
will  be  shown  to  provide  a  more  complete  description. 

Having  obtained  state-to-state  cross  sections  which  neglect  the  laser  perturbation  of 
the  projectile,  this  perturbation  will  then  be  included  within  the  collisional  Schrodinger 
equation  in  the  form  of  Volkov  dressed  states  (Volkov  1935.  Kroll  and  Watson  1973, 
Geltman  and  Macquet  1989).  For  a  single  state-to-state  transition,  with  a  specific  final 
projectile  energy,  the  Volkov  dressed  states  exercise  a  very  significant  effect  on  the  cross 
sections.  When  the  state-to-state  cross  sections  are  summed  over  all  possible  final  projectile 
energies  reached  through  absorption  and  emission  of  photons,  this  work  however  shows  that 
the  cross  sections  that  include  the  laser  perturbation  of  the  projectile,  are  only  marginally 


different  from  the  cross  sections  that  have  neglected  the  perturbation  of  the  projectile. 

In  order  to  synchronize  the  time  frame  of  the  relative  orbit  of  the  collisional  species 
with  the  time  frame  of  the  laser  field,  a  phase  shift  8W  must  be  introduced  into  the  vector 
potential  A  of  the  laser  field.  Apart  from  Smith  and  Flannery  (1991)  this  phase  shift 
has  been  ignored  in  all  the  above  work,  and  has  important  effects.  Bachau  and  Shakeshaft 
(1984)  have  explicitly  acknowledged  the  phase  shift  for  2 p  excitation  in  H*  —H(l$)  inelastic 
scattering  in  a  nearly  resonant  field. 

This  paper  has  also  been  able  to  describe  in  greater  detail  the  role  played  by  the 
Floquet  states,  so  as  to  give  a  very  concise  description  of  the  interaction,  in  terms  of 
the  probability  amplitudes  and  the  locations  of  photon  absorption  and  emission  during 
the  collisional  sequence  of  events.  In  section  4  the  photon  absorptions  and  emission  are 
discussed  in  the  form  of  (A,hu>)  and  (e,A,ftu>),  as  recently  described  (Newell  1990). 
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2.FLOQUET  DRESSED  STATES  OF  AN  ATOM  IN  A  LASER  FIELD 

This  theoretical  model  is  based  on  a  typical  crossed  beam  experiment,  in  which  the 
laser  beam,  projectile  beam  and  atomic  particle  beam  all  cross  at  one  point  in  space  at 
90°  angles  to  each  other.  The  collision  will  be  considered  in  the  center  of  mass  frame  of 
reference,  but  it  should  be  noted  that  the  projectile  will  have  a  much  greater  velocity  than 
the  target  atom  in  the  laboratory  frame  of  reference.  This  is  important  since  the  atom  will 
be  inside  the  laser-beam  for  the  duration  of  the  collision,  and  so  must  be  excited  from  one 
dressed  atomic  state  to  another. 

When  the  laser  perturbation  of  the  projectile  is  neglected  the  Hamiltonian,  in  the 
center  of  mass  frame  of  reference,  is  given  by 

B  =  +  E  !(P"  +  A/c)2  +  +  F<R’r>  (!) 

i 

The  channel  coordinates,  (P#,R),  represent  the  momentum  operator  and  the  position 
vector  for  the  projectile.  The  internal  coordinates,  (pr,,rt),  represent  the  momentum 
operator  and  the  postion  vector  for  each  bound  electron.  The  electronic  coordinates  are 
denoted  collectively  by  r,  c  is  the  speed  of  light,  /j  is  the  projectile-target  reduced  mass, 
V(r)  is  the  internal  potential  for  the  bound  atomic  electrons  and  V(R,  r)  is  the  external 
projectile-target  interaction  potential.  Atomic  units  are  used  throughout. 

The  monochromatic  laser,  of  frequency  u>  and  phase  6^,  should  be  exactly  described 
by  the  vector  potential 


A(w,rl_../)  =  A0  cos^.  •  i\.  tu )1  +  K-) 


(2) 


But  at  the  atomic  scale,  since  the  reduced  wavelength  Au.  =  J-  is  much  larger  than  atomic 

A  w 
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dimensions,  the  dipole  approximation  is  valid  so  that  the  vector  potential  reduces  to 


A  =  A0  cos(u;f  +  6U) 


(3) 


In  semiclassical  scattering  theories,  the  point  of  closest  approach  (the  periapsis),  between 
the  projectile  and  the  target,  is  generally  defined  as  occurring  at  time  t=0.  If  the  vector 
potential,  under  the  dipole  approximation,  was  written  as  A  =  A0  cos(wt),  then  this  would 
indicate  that  for  each  collision  the  electric  field  strength  was  at  a  maxima  of  its  cycle  at 
the  point  of  closest  approach,  as  indicated  by  the  dotted  curve  in  Fig.l.  Since  this  is  not 
in  general  true,  it  is  necessary  to  use  the  vector  potential  given  in  eq.(3),  where  is  the 
phase  shift  of  the  vector  potential  eq.(3),  from  A0  cos(u>f),  at  the  point  of  closest  approach 
as  indicated  by  the  solid  curve  in  Fig.l. 

To  dress  the  atomic  states  it  is  necessary  to  consider  the  effect  of  the  laser  field  on 
the  isolated  atomic  target.  This  is  done  by  removing  the  projectile  components  from  the 
Hamiltonian  eq.(l).  On  applying  the  unitary  transformation  U  =  exp (i  A*r;/c),  (Reiss 
1970),  the  Schrodinger  equation  for  an  atom  in  a  laser  field  is  then  expressed  as 


(4) 


where  e0  is  the  electric  field  strength  A0u;/c,  and  where  it  has  been  assumed  that  the 
periodicity  is  the  only  time-dependent  part  of  the  vector  potential.  The  solution  of  eq.(4) 
is  the  dressed  atomic  state  4'1  (at  present  unknown).  On  using  the  field  free  atomic 
Hamlitonian  H0  —  ^  +  F(r),  and  on  setting  /.if  ~  -ile0' ,  the  laser  perturbed 

atomic  Schrodinger  equation  eq.(4)  is  written  in  the  more  compact  form 


(5) 
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which  is  now  suitable  for  application  of  the  Floquet  prescription  (Shirley  1965,  Sambe 
1973,  Chu  1985).  Under  this  prescription,  the  solution  to  eq.(5)  is  of  the  form 

IV-')  =  K-;>  =  (6) 

where  Qq  is  the  energy  or  “quasi-energy”  of  the  new  dressed  state,  and  Gq  has  a  periodic 
time-dependence  of  period  Since  Gq  is  periodic  in  time,  it  is  expanded  by  the  Fourier 
series 

|G,)  =  £  |G»  (7) 

n 

where  the  Floquet  notation,  introduced  by  Shirley  (1965)  is  used.  Thus 

|V’,  n)  —  |V>)|«)  with  (f.|?i)  =  emwi  (8) 

and  \n)  will  be  called  a  “periodic”  state  for  clarity  later.  Using  this  Fourier  expansion,  the 
Schrodinger  equation  eq.(5)  can  now  be  rewritten  so  as  to  explicitly  include  the  periodic 
photon  dependence  nhu>  as 

Wff0  +  nKa>  +  X>;V“'  -  c-i"))e-««'/l|G;,n>  =  £  l<?»  (9) 

n  i  /  n 

In  this  weak-field  regime  the  dressed  atomic  states  can  be  approximated  by  a  combination 
of  field  free  atomic  states,  as 

W'q)  =  A?a  1^9’”)  (10) 

n 

where  the  periodic  time-dependence  of  the  coefficients  has  been  separated  by  using  the 
periodic  states  \n).,  On  comparing  this  with  the  Floquet  dressed  atomic  states 

K>  =  «-  ,9'm  £  in) 
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then  the  laser  perturbed  Schrodinger  equation  eq.(9)  can  be  written  in  the  matrix 
representation  as 


E 


(ca  '^^’^)^/3a^n=Tn 


+  Sr, 


■1  M/3a^7l=m+1 


A”q  =  QqA 


m 

90 


(12) 


The  field  free  atomic  energy  levels  ea  and  the  laser  induced  off-diagonal  elements  fipa  = 
—i^(4>/3\e0-r\(f)a)e±t6u  are  known,  so  that  the  matrix  defined  by  eq.(12)  can  be  diagonalized 
to  obtain  the  quasi  energies  Qq  and  the  dressed  atomic  states  V’g-  Because  the  Floquet 
prescription  allows  for  an  infinite  set  of  periodic  states,  from  n  =  —  oo  to  oc,  it  is  necessary 
to  truncate  the  matrix.  But  as  long  as  the  off-diagonal  couplings  fi^a  are  small,  this 
truncation  can  be  quite  severe  whilst  still  allowing  convergent  dressed  states  to  be  obtained. 
Fig.2  shows  a  specific  example  of  a  Floquet  matrix  for  a  hydrogen  atom  in  a  laser  field. 
Significant  couplings  only  occur  between  field  free  atomic  states  <f>p  and  <j>a  if  the  energy 
level  separation  is  of  the  order  of,  or  less  than,  the  photon  energy  of  the  laser 


|Ca  -  ep\  <  fcw 


(13) 


For  e  collisional  problem  that  considers  excitation  from  the  n=l  to  the  n=2  sub-levels, 
internal  couplings  within  the  n=2  manifold  need  only  be  considered  at  the  photon  energies 
used  within  this  work. 

To  a  first  order  approximation,  which  allows  only  single  photon  coupling,  the  matrix 
of  Fig.2  yields  the  dressed  atomic  states 

|  S)  =  A\2S,n)  +  J3f,{‘|2P0,r»-t-l)  +  Dc^UP^n-l)  (14a) 

| P)  -  A\2P0,n}  +  Be'e“\2S,n  +  l)  +  Be~l6*\2S,n-l)  (146) 
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It  is  emphasized  for  these  dressed  atomic  states,  that  n  varies  across  the  range  from 
n  =  -oo  to  oo,  giving  rise  to  an  infinite  set  of  periodic  Floquet  states.  In  higher  order 
approximations  to  the  dressed  states,  multiphoton  couplings  are  included,  so  that  eq.(14 
a,  b)  include  further  terms  for  n±2,  n±3  etc.  Results  of  section  5  shows  that  convergence 
is  achieved  by  a  third  order  approximation,  so  that  terms  up  to  n±3  provide  an  adequate 
description  of  the  dressed  atomic  states,  over  the  electric  field  strength  range  considered 
in  this  work. 

From  the  form  of  eq.(14  a,  fc),  it  can  be  seen  that  the  fractional  2S  and  2Po  components 
of  both  dressed  states  S  and  P  change  with  varying  electric  field  strength.  The  character 
of  the  dressed  states  S  and  P  will  then  change  with  increasing  e0.  This  changing  character 
of  the  dressed  states  is  clearly  seen  from  the  radial  distribution  functions,  where  the  radial 
distribution  functions  Ds(t’)  and  Dp(r)  are  defined  by  the  equations 


(S\S)  = 

poo  p2ir  pi r 

/  Ds(r)dr  /  /  dcosddcj) 

Jo  Jo  Jo 

(15o) 

(P\P)  = 

poo  p2i:  pn 

/  Dp(r)dr  /  /  d  cos  Od(f> 

Jo  Jo  Jo 

(156) 

The  radial  distribution  functions,  for  the  dressed  states  S  and  P,  are  shown  in  Fig. 3a 
-  Fig. 3c  as  the  electric  field  strength  e0  increases  from  0  au  to  3  x  10-3  au.  These 
radial  distribution  functions  correspond  to  third  order  Floquet  dressed  hydrogen  atoms. 
Convergence  in  the  state-to-state  cross  sections  of  section  5,  for  third  order  Floquet  dressed 
atomic  states,  is  used  to  justify  a  third  order  approximation  here.  It  is  seen  that  as  the  field 
strength  increases,  the  radial  distribution  function  Ds{r )  loses  its  intrinsic  2S  character, 
while  gaining  an  increasingly  2Po  character.  In  a  similar  fashion  the  radial  distribution 
function  Dp(r)  is  losing  its  intrinsic  2P0  character,  while  gaining  an  increasingly  2S 
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character.  The  plot  in  Fig.4  attempts  to  represent  this,  by  looking  at  the  percentage 
component  of  the  2Pq  field  free  radial  distribution  function,  D2P0{r),  in  both  Ds{r)  and 
Dp(r).  This  shows  very  clearly  the  play  off  between  the  dressed  states  S  and  P  indicated 
in  the  set  of  plots  in  Fig. 3.  As  the  electric  field  strength  rises,  both  states  S  and  P 
simultaneously  give  up  some  of  their  intrinsic  zero  field  character,  and  assume  more  of 
each  other’s  intrinsic  zero  field  character.  This  character  swapping  continues  past  the 
point,  where  the  dressed  state  S  has  more  of  a  2Po  character  than  the  dressed  state  P. 
If  the  statement  is  made,  for  the  collisional  energies  studied  in  this  work,  that  the  field 
free  lS-2Po  state-to-state  integral  cross  section  is  always  larger  than  the  field  free  1S-2S 
state-to-state  integral  cross  section,  ie  if 

^lS,2Po  >  OJ S.25  (16) 

then  from  the  plot  of  Fig.4  it  might  be  reasonable  to  expect,  as  the  electric  field  strength 
rises,  that  the  laser  assisted  1S-2S  state-to-state  cross  section  will  rise,  and  the  laser  assisted 
lS-2Po  state-to-state  cross  section  will  fall.  This  rise  and  fall  can  also  be  expected  to 
continue  past  the  point  where  the  laser  assisted  1S-2S  state-to-state  cross  section  is  greater 
than  the  laser  assisted  lS-2Po  state-to-state  cross  section. 
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3.TIME-INDEPENDENT  SCATTERING  CROSS  SECTIONS 


In  section  2  the  time-dependent  Schrodinger  equation  for  an  atom  in  a  laser  field, 
was  solved  using  the  Floquet  prescription  to  obtain  new  eigenstates,  or  dressed  states 
for  the  atom  in  a  time-independent  form.  The  time-dependence  has  been  removed  by 
allowing  for  absorption  and  emission  of  photons.  These  dressed  states  can  now  be  used 
to  solve  collisional  problems  in  a  manner  directly  analogous  to  field  free  collisions.  In 
this  work  the  scattering  cross  sections  for  the  laser  assisted  collisions  are  derived  from 
the  Lippman  Schwinger  equation  (Bransden  1970).  This  approach  is  valid  to  all  orders 
of  approximation  of  the  dressed  states,  given  the  necessary  assumption  that  the  dressed 
states  form  a  complete  and  normalized  set.  From  section  2  the  laser  assisted  collision  can 
be  described  by  the  Hamiltonian 

H  =  ±P2r  +  ^Pri  -  lbo  •  -  e-4*1)]  +  V(r)  +  F(R,r) 


*  (17) 

=  Hf  +  V{  R,r) 

When  the  interaction  potential  V(R,  r)  is  removed,  the  remainder  Hp  describes  the 
projectile  far  removed  from  the  target  atom,  which  remains  in  the  presence  of  the  laser 
field.  The  solution  $  of  the  Schrodinger  equation 


(18) 


is  the  combination  =  V’g  exp(ik9  •  R)  of  the  dressed  atomic  states  if*1  of  section  2,  with 
a  field-free  projectile  plane  wave.-  The  scattered  collisional  wave  function  in  the  presence 
of  V  then  satisfies 


I*?**)  =  e' 


■iEqt/h 


|$„)  -  lim 


1 


f—0  Eq  ±  i(  --  H0 


Y  (  —  tEqt/hxft/k 


(19) 
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The  probability  of  a  transition  between  the  scattered  collisional  wave  functions  and 
W"  is  given  by  the  S-matrix 

««(*)  =  (20) 

On  inserting  the  scattered  collisional  wavefunction  eq.(19)  into  the  S-matrix  eq.(20),  and 
expanding  the  dressed  atomic  states  ifr1  into  the  combination  of  field  free  atomic  states 
of  eq.(10),  the  S-matrix  is  then 

Spq{t)  =  (^pl^?)^p? 

m  (3 

X  ?  5  Aqa  Qq  +  Tjki  -  nhu:  +  U-QP 
x  e~ilQ*+-kk'-nhu>+i()i/hy\ii,+} 

Under  the  following  definition 

^exp[i(Qp  +  ±k2f  -  Qq  -  +  {n-m)hu  -  u)*/fi] 

«“*°  Q9  +  ^t2  -  rcfcw  +  u  -  -  j^k2f  +  mhu;  (22) 

=  2 7T  i«(Qp  +  ^jfc2  +  (n--m)fru>  -  Qq  -  ±k2x) 
of  the  Dirac  ^-function,  eq.(21)  becomes 

Spq(t)  =  p\^  q)^pq  +  *27r  7Pq  6(QP  +  ^hf  -f  (n— rri)fiu>  —  Qq  —  y-k, )  (23) 

Using  terminology  directly  analogous  to  field  free  collisions,  the  S-matrix  eq.(23)  has  been 
written  in  terms  of  the  transition  amplitude 

T„  =  E  E  E  E  e*p(i(k.  -  M  ■  R)  (24) 

rt  q  m  (3 


-  +  m^w 

(21) 
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and  includes  all  of  the  possible  photon  absorptions  and  emissions  that  can  occur  during 
collisonal  transitions  between  the  dresses  states.;  For  a  field  free  collision,  kq  and  kp  would 
be  defined  as  the  wave  numbers  for  the  projectile,  associated  with  the  initial  atomic  state 
4’'q  and  the  final  atomic  state  ,  and  the  difference  between  kq  and  kp  depends  on  the 
difference 

Qq  —Qp  =  2j;kp  ~  (25) 

between  the  energy  levels  of  ip'q  and  Vp-  The  initial  asymptotic  wave  number  for  a  laser 
assisted  collision  is  still  kq,  since  the  projectile  has  not  yet  interacted  with  the  atom, 
but  the  final  wave  number  kf  differs  from  kp  so  as  to  satisfy  the  conservation  of  energy 
requirement 

jfjj kf  —  Qg  +  277 hq  ~  Qp  ~  (n-m)hoj  (26) 

Since  kp  is  defined  by  eq.(25),  kf  can  be  defined  by 

~  2p^p  {n  —  Tn)hu}  (27) 

This  shows  that  kf  can  differ  from  kp  by  the  photon  energy  (n-m)huj,  which  corresponds 
to  the  number  of  photons  absorbed  or  emitted  in  the  state-to-state  transition. 

The  S-matrix  obtained  in  this  treatment  describes  scattering  within  the  dressed  basis 
set.  The  associated  scattering  amplitudes  which  can  be  derived  from  various  theoretical 
techniques,  only  describe  transitions  between  dressed  states  that  are  part  of  the  new 
basis  set.  To  obtain  results  that  can  be  directly  compared  with  field  free  collisions  and 
experimental  results,  it  is  necessary  to  calculate  the  laser  assisted  scattering  cross  sections 
with  respect  to  the  field  free  basis  set.  From  Bransden  (1970)  the  cross  section  for  a 
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transition  from  the  field  free  atomic  state  <j)a  to  field  free  atomic  state  <f>p  is 


(Tpa  —  2^3  J  dtl  J  dkf  hj  Wp a  (28) 

where  Wpa  is  the  transition  rate.  For  any  given  scattering  angle  and  state-to-state 
transition,  there  can  be  more  than  one  wave  number  for  the  projectile,  corresponding 
to  absorption  and  emission  of  energy  from  the  laser  field.  This  means  that  the  transition 
rate  is  defined  as 

we„  =  £  2V  |rfl„|2  6(te  -  +  ifcj  -  ±kl  +  6E)  (29) 

fcy 

where  it  is  written  in  terms  of  the  transition  amplitude  Tpq,  and  includes  the  unknown 
energy  change  8E  due  to  the  influence  of  the  laser.  In  eq.(29)  the  final  wave  number  kf 
must  be  written  in  terms  of  the  wave  number  kp,  associated  with  the  atomic  state  <f>p  in 
a  field  free  collision.  Using  the  definition  of  kp 

Yjjkp  =  tp  +  ±-h^  -  eQ  (30) 

in  the  delta  function  of  eq.(29),  the  resulting  equation 

hk)  =  ~ SE  (31) 

relating  kf  to  kp  is  very  similar  to  eq.(27),  and  includes  the  unknown  energy  change  6E.  It 
is  now  possible  to  insert  the  transition  rate  eq.(29)  into  the  integral  cross  section  eq.(28). 
Using  \kf\dkf  —  d(^kj)  so  that  the  integration  over  kf  can  be  evaluated,  the  integral  cross 
section  becomes 

»/>.  =  jTTJJu  j dSlY,kf  2V  ! W  *(  £*/  ~  frl  +  iE)  (32) 

hi 
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which  is  equivalent  to  the  standard  integral  cross  section  in  the  field  free  basis  set,  except 
that  allowance  has  been  made  for  the  absorption  or  emission  of  an  unknown  amount  of 
energy  6E  from  the  laser  field  during  the  collision.  The  set  of  dressed  atomic  states  obtained 
in  section  2  are  assumed  to  be  a  complete  and  normalized  set  so  that  Y ]q  \*l>q){4'q\  =  1*  On 
using  this  projection  operator  the  “measured”  transition  amplitude,  taken  with  respect  to 
the  field  free  basis  set,  is  written  in  terms  of  the  “dressed”  transition  amplitude,  taken 
with  respect  to  the  dressed  basis  set  as 

=  (*„|V|*„)  =  Y,Y  to.OIlM  Tn  (V-,16.,0)  (33) 

V  9 

The  “dressed”  transition  amplitude  Tpq  of  eq.(24)  has  been  determined  by  the  Lippman 
Schwinger  equation,  and  according  to  eq.(26),  kf  must  satisfy  the  conservation  of  energy 
requirement 

Qp  +  hk)  -  Q*  -  hkl  +  nhu  = 0  (34> 

where  n  has  replaced  m  —  n.  As  demonstrated  in  eq.(33)  the  “measured”  transition 
amplitude,  T$a,  can  be  obtained  from  the  “dressed”  transition  amplitude,  Tpq,  and  kf 
must  satisfy  the  second  conservation  of  energy  requirement 

e/3  +  2^/  -  Cor  -  =  0  (35) 

obtained  from  eq.(30).  The  final  wave  number,  kf,  of  the  scattered  projectile  must  be 
identical  for  both  equation  eq.(34)  and  eq.(35).  This  provides  an  opportunity  to  determine 
the  total  energy  absorbed  or  provided  by  the  laser  field.  It  may  be  a  natural  assumption 
that  this  energy  bE  is  an  integer  number  of  laser  photon  energies  nhuo.  But  on  diagonalizing 
the  Floquet  matrix  to  obtain  the  dressed  wave  function,  it  can  be  seen  for  some  transitions 
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to  occur,  that  the  energy  change  SE  cannot  be  an  integer  number  of  laser  photon  energies. 
In  molecular  spectroscopy  an  analogous  situation  occurs,  where  a  molecule  is  excited  to 
a  higher  energy  level  by  a  laser  photon,  which  has  more  energy  than  is  required  for  the 
transition  (Weissbluth  1989).  In  the  Raman  effect,  the  molecule  absorbs  the  laser  photon 
and  then  emits  a  new  Raman  photon  to  remove  the  excess  energy.  This  is  the  process  that 
is  assumed  to  be  occuring  whenever  necessary,  with  no  additional  restrictions  applied  to 
the  transition  because  it  must  proceed  via  the  Raman  effect.  A  discussion  of  this  effect  is 
given  in  section  4,  and  is  shown  to  play  a  significant  role  in  the  laser  dependent  changes 
in  the  integral  cross  sections. 

Since  the  vector  potential  eq.(2)  can  only  be  defined  to  within  a  phase  shift  6a,,  as  in 
Fig.l,  it  is  therefore  necessary  to  take  the  following  average 

1  f27r 

^experiment  =  J  d(6o,)  (36) 

of  a  over  the  range  of  6a,  from  0°  to  360°. 
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4.THE  ROLE  OF  THE  FLOQUET  DRESSED  STATES 

In  discussing  the  use  of  Floquet  theory  in  laser  assisted  collisions,  reference  will  be 
made  to  a  reduced  Floquet  system  of  a  hydrogen  atom  in  a  laser  field.  An  abbreviated 
energy  level  diagram  is  shown  in  Fig.5;  The  states  shown  are  only  a  limited  number  of 
the  infinite  set  of  periodic  Floquet  states.  In  a  first  order  approximation  which  allows  for 
single  photon  coupling  only,  the  system  is  adequately  described  by  the  following  7  states 


selected  in  the  diagram  of  Fig.5, 

|1)  =  |1S,0)  (37a) 

|2)  -  A\2S  -1}  +  Bci6“  j2Po,0)  +  Be~ib“  |2.P0,-2)  (376) 

|4)  =  A|2S,0)  +  Bei6“  |2P0,+1>  +  Be~i6“  |2P0,-1)  (37c) 

|6)  =  A|25,+l)  +  Bei6w  |2Po>+2)  +  Be”16*  |2Po-0)  (37 d) 

|3)  =  A|2P0,-1)  +  Be*-  \2 S, 0)  +  Bc~lb-  \2S-2)  (37e) 

|5)  =  A\2P0,0)  +  Bei6“\2S,+l)  +  Be~i6“\2S-l)  (37/) 

1 7)  =  A|2P0,+1)  +  Bei6“  125, +2)  +  Be~l6“  |2S,  0)  (375) 


For  illustration  purposes  a  first  order  approximation  has  been  used  here  in  the  atomic 
dressing  to  reduce  the  number  of  possible  transitions  available;  but  all  that  is  said  here  is 
valid  to  any  order  of  approximation 

For  the  purposes  of  numerical  calculations,  the  energy  levels  2S  and  2Po  can  be 
considered  to  be  degenerate,  though  it  must  be  kept  in  mind  that  this  is  not  so.  Even  if 
the  field  free  energy  levels  were  degenerate,  this  degeneracy  would  be  broken  when  in  the 
influence  of  the  laser  field  due  to  coupling  between  the  IS  and  2Po  states. 
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Consider  the  laser  perturbed  Schrodinger  equation  eq.(4) 


(?[ 


3  2  •  3  . 


r,(e 


_  —  t(uit+5w) 


+ 


V(r)V 


(38) 


The  laser  induced  perturbation  \e0  •  rt,  provides  a  dipole  coupling  between  previously 
orthogonal  field  free  atomic  states,  a  coupling  that  contains  a  time-dependence  elwt. 
Because  of  this  time-dependent  coupling,  the  field  free  states  are  no  longer  eigenfunctions 
of  the  Schrodinger  equation,  and  the  Schrodinger  equation  is  then  no  longer  time- 
independent.  In  section  2  it  was  shown  that  the  Floquet  Description  results  in  a  new 
orthogonal  basis  set,  composed  of  dressed  atomic  states  whi  an  be  treated  in  a  time- 
independent  manner.  These  three  properties  (orthogonality,  eigen  functions  and  time- 
independence)  of  the  Floquet  dressed  atomic  states,  allow  traditional,  field  free,  time- 
independent  scattering  theories  to  be  applied  to  laser  assisted  collisions.  However  the 
price  to  be  paid  for  these  three  properties  is  the  introduction  of  the  infinite  set  of  periodic 
solutions.  The  transition  between  the  finite  set  of  field  free  atomic  states,  and  the  infinite 
set  of  dressed  atomic  states,  is  achieved  by  projection.  This  has  been  introduced  in  eq.(33) 
of  section  3,  where  the  “measured”  transition  amplitude,  taken  with  respect  to  the  field  free 
basis  set,  is  written  in  terms  of  the  “dressed”  transition  amplitude,  taken  with  respect  to 
the  dressed  basis  set.  Because  the  dressed  atomic  states  have  a  periodic  time-dependence 
it  is  necessary,  in  order  to  achieve  this  projection,  to  write  the  field  free  atomic  states  as 
a  combination  of  field  free  atomic  states  and  periodic  states.  In  the  expansion 


j<ftq)  =  ,nw<i  4>a,n)  (39) 

n 

of  the  field  free  atomic  state  into  a  field  free  state  and  periodic  state  combination,  An  is 
non-zero  only  when  n  =  0.  Thus  the  field  free  atomic  state  j<pQ)  is  equivalent  to  the  n  =  0 
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Floquet  state  |<£a,0).  Using  this  with  the  above  set  of  dressed  atomic  states  eq.(37),  the 
new  dressed  basis  set  (in  this  first  order  approximation),  for  a  hydrogen  atom  originally 
in  the  2S  field  free  state,  for  example,  will  have  a  non-zero  population  density  only  in  the 
dressed  atomic  states  labelled  3,  4  and  7.  Conversely,  any  population  density  in  the  dressed 
atomic  states  3,  4  and  7  of  the  dressed  basis  set,  will  each  contribute  to  the  population 
density  of  the  2S  state  of  the  field  free  basis  set.  In  order  for  population  densities  in  the 
dressed  states  3  and  7  to  contribute  to  the  population  density  in  the  field  free  state  2S,  the 
laser  field  is  required  to  provide  (or  absorb)  a  non-integer  number  of  laser  photon  energies. 
This  is  because  the  2S  and  2P  field-free  energies  are,  in  fact,  not  degenerate,  and  because 
Es  =  E?,  —  hu  and  E^  =  E$  +  hu,  as  asserted  by  Floquet  theory.  Although  it  is  not 
immediately  obvious,  the  same  problem  is  associated  with  the  dressed  state  4.  Due  to  2S- 
nP  (n>2)  couplings,  the  energy  of  state  4  shifts  with  changing  electric  field  strengths  so  that 
j ELj  ^  Eis-  As  mentioned  in  section  3,  production  of  these  non  integer  photons,  is  known 
as  the  Raman  effect,  and  as  can  be  seen,  a  Raman  photon  will  be  involved  whenever  the 
“dressed”  transition  amplitudes  are  projected  on  to  the  “measured”  transition  amplitudes. 
This  is  considered  an  atom-field  effect  and  will  not  involve  the  projectile  directly.  For 
collisions  within  the  dressed  basis  set,  all  energy  differences  for  atomic  transitions,  are 
accounted  for  by  changes  in  the  wave  number  of  the  projectile,  plus  the  absorption  or 
emission  of  an  integer  number  of  laser  photons. 

To  show  the  effect  of  the  Floquet  approach  on  a  laser  assisted  collision,  the  1S-2S 
excitation  in  an  electron- hydrogen  collision 

H{IS)  +  t~  +  Nfiu)  —  H(2S)  4  f’ t  Mtua  (40) 
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will  be  considered  for  first  order  dressed  atomic  states.  Table  1  shows  the  “dressed” 
transition  elements  and  their  probability  amplitudes  when  projected  onto  the  “measured” 
transition  amplitude  Tis2S-  These  values  can  then  be  used  in  a  first  order  Born- Wave 
approximation.  It  is  interesting  to  note  that  the  3  “dressed”  transition  amplitude  elements 
for  T31  are  the  same  as  the  3  “dressed”  transition  amplitude  elements  for  T71 ,  with  identical 
changes  in  the  projectile  energy.  In  fact  the  transition  amplitude  elements  for  T31  are 
duplicated  for  all  transitions  1  — *  P  where  P  represents  all  states 

\P)  =  A\2P0,n)  +  Bei6“\2S,n  +  l)  +  Be-i6“\2S,n-l)  (41) 

This  feature  of  the  Floquet  dressed  atomic  states  is  very  useful  for  the  closely  coupled 
calculations  required  in  the  Multichannel  Eikonal  Treatment  of  the  collision,  where  the 
infinite  set  of  periodic  states  tend  to  make  such  calculations  appear  at  first  sight  prohibitive. 
Consider  two  members  of  the  set  P  of  states  defined  by  eq.(41), 


|  ifip,  1),  energy  level  =  Qp  +  hu> 

and 

IV’p,0),  energy  level  =  Qp 

In  a  transition  between  these  two  states,  the  change  in  projectile  energy  due  solely  to 
the  difference  in  energy  levels  is  ?Lu>,  ie  ~  Then  (V’P, l'|V’p> 0)  =  0, 

since  a  photon  must  be  absorbed  from  the  laser  field  in  the  transition,  a  photon  which 
neither  nor  |v’p,0)  can  absorb  yet  still  preserve  their  energy  levels.  If  however  each 

state  is  associated  with  a  projectile  of  the  same  energy,  forming  a  collisional  wavefunction 
=  il'pcxp(ik  •  R),  then,  for  a  collisional  transition  from  atomic  state  \u'p,0)  to 
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the  projectile  can  absorb  the  laser  photon  as  is  required,  and  regain  the  energy  hu>  which 
it  lost  in  making  the  transition  to  the  higher  energy  level.  Hence  the  projection  between 
the  two  collisional  wave  functions  is 


1|  exp(— ik/  •  R)  exp(ik,-  •  R)  \tpp,0)  =  1  (42) 


ie  the  wave  functions  are  no  longer  orthogonal  when 


J1.2  _ 
2 //  “ 


1  i,2  t, ,  _  1  u 2  _ 

27>1  ~  ~  2p  pO  ~ 


-Li-2 

2p  » 


(43) 


This  result  can  be  extended  to  all  other  members  of  set  P,  with  the  same  results  for  other 
sets  of  periodically  repeating  dressed  atomic  states.  This  feature  of  Floquet  dressed  states 
makes  it  impossible  to  define  the  probability  amplitude  of  being  in  a  given  dressed  atomic 
state  when  combined  with  a  projectile  wave  function.  Or  described  another  way  the  system 
exists  in  a  state  of  flux  between  all  periodic  states  of  the  same  form  associated  with  a  given 
projectile  energy..  This  necessitates  a  summation  of  the  probability  amplitudes  over  what 
are  termed  “grouped”  states  (Smith  1991).  Hence  the  entire  infinite  set  of  periodic  Floquet 
states  can  be  included  exactly  within  the  closely  coupled  calculation.  Although  essential 
for  the  Multichannel  Eikonal  Treatment  these  “grouped”  states  can  be  ignored  in  the  first 
Born  approximation. 

Inserting  the  contents  of  Table  (1)  into  eq.(33),  which  relates  the  “measured” 
transition  amplitude  to  the  “dressed”  transition  amplitude,  T\S2S  becomes 


Tists  =  4t2W  =  £*|s) 

+  ABV,S2P'C-'‘-  i(±k]  =  ±kls  -  M 
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+  ABV^e*-  K±k)  =  £tjs  +  M 
+  ABVtS2Pe{e“‘  +  e-“")  S(±k)  =  i^Po) 

+  b2vIS2S(i  +  r'2‘~)H±kj  =  &*»,  -  M 

+  B2  W(ei2‘-  +  1)  l(±k)  =  i^Po  +  ku>)  (44) 

In  contrast  the  laser  assisted  differential  cross  section  can  be  calculated  using  the  traditional 
time-dependent  Perturbation  theory  (Bayfield  1979).  For  an  oscillatory  perturbing 
potential 

H(r,t)  =  — ica  •  r(eia,V6“  -  e~iloi c~i6u)  (45) 

the  perturbed  or  dressed  wave  function  is 

OO 

<  =  EE  <>«(')  +t  («) 

n=0  /3 

where  is  a  summation  over  the  unperturbed  atomic  states  and  tp  is  the  energy  level 
of  state  To  a  first  order  approximation  this  gives  just  three  dressed  states  for  the  laser 
perturbed  IS,  2S  and  2Po  field  free  states 


<Pis  =  \fos)e-Ulsi/h  (47a) 

<t>2S  =  (  Dj  ms)  +  D2\4>2Po)eiutetS-  +  D2\4>2P<>)eriui e~^  \e-^t/h 
4>m  =  (  D,\<p2Po}  +  D2\4>2s)eiwiel6w  +  D2\4>2s)e~iuJt (47c) 

where  Dj  =  N,  D2  =  —i(<Ps\t0  ■  rj 0a)~r  and  N  is  a  normalization  constant.  Instead  of 
projecting  the  dressed  states  onto  field  free  states,  this  form  of  approach  assumes  that  the 
dressed  states  collapse  onto  their  field  free  source  (Sharmaand  Mohan  1990).  Taking  this 
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type  of  approach  the  T1S2S  transition  amplitude  becomes 


Has  =  DtKsis  =  £&) 

+  DA\s2P0'-“-  f(£k}  -  £*5S  -  M 

+  i>2  W„*“-  =  £&  +  M  (46) 

Comparing  the  Floquet  dressed  atomic  states  of  eq.f37)  with  the  Perturbative  dressed 
atomic  states  of  eq.(47),  it  can  be  seen  that  each  individual  dressed  state  has  a  very 
similar  form  under  both  treatments.  But  the  Perturbative  approach  does  not  predict  the 
infinite  set  of  periodic  states.  This  lack  of  an  infinite  set  of  states  makes  a  great  difference  in 
the  “measured”  transition  amplitude,  as  can  be  seen  by  comparing  the  Floquet  transition 
amplitude  of  eq.(44)  with  the  Perturbative  transition  amplitude  of  eq.(48).  Not  only  does 
the  Floquet  transition  amplitude  have  more  terms,  or  more  paths  for  the  transition,  it  also 
has  a  Su  dependence  that  will  be  preserved  in  the  cross  sections,  whilst  the  f>w  dependence  of 
the  Perturbative  transition  amplitude  will  clearly  cancel  out  in  the  resulting  cross  sections. 

By  first  predicting  the  infinite  set  of  states,  and  then  by  writing  the  '‘measured” 
transition  amplitude  in  terms  of  the  “dressed”  transition  amplitude,  this  Floquet  model 
provides  two  distinct  areas  of  possible  photon  absorption  and  emission,  which  can  be 
denoted  by  the  mechanisms  (A ,/iu>)  and  (e,A,/iu>).  (A,fuu)  refers  to  the  projection  of  the 
dressed  atomic  states  onto  the  field  free  atomic  states,  and  gives  rise  to  laser  photon 
absorption  and  emission  and  Raman  photon  emission.  (e,A,h.u>)  refers  to  the  transitions 
between  dressed  atomic  states  which  requires  the  influence  of  the  projectile,  and  which  only 
gives  rise  to  laser  photon  absorption  and  emission.  There  is  a  third  area  of  possible  photon 
absorption  and  emission,  which  is  omitted  in  this  model  that  ignores  the  laser  perturbation 
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of  the  projectile.  This  can  be  denoted  by  the  mechanism  (e,fi.u>),  which  refers  to  photon 
absorption  and  emission  by  the  projectile  before  and  after  the  collision.-  When  the  full 
laser  perturbation  is  included,  (section  6)  there  is  no  indication  of  such  photon  absorption 
and  emission,  so  that  the  multiphoton  events,  associated  with  the  projectile,  occur  during 
transitions  between  dressed  states  only,  or  is  described  by  (e,A,/iw).  Hence  (A,hu>)  and 
(e,A,ku>)  adequately  describe  the  mechanisms  and  sequence  of  all  multiphoton  events.  It 
can  therefore  be  seen  that  this  model  provides  a  concise  description  of  the  projectile  atom 
collision  in  a  laser  field,  where  the  interactions  (A ,hu>)  and  (e,A,fiu>),  as  recently  described 
(Newell  1990)  are  included  exactly. 
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5.RESULTS:  e~  +  H{1S)  +  Nhu> 


e~  +H(2S,2Po)  +  Mhu> 


State-to-state  cross  sections  for  the  1S-2S  and  lS-2Po  transitions,  in  laser  assisted 
e~-H(lS)  collisions,  are  presented  in  both  the  Multichannel  Eikonal  Treatment  and  the 
Born- Wave  approximation,  as  a  function  of  impact  energy  and  laser  field  intensity.  All 
cross  sections  shown  are  a  summation  over  the  state-to-state  cross  sections,  for  all  possible 
final  projectile  energies  consistent  with  a  given  initial  relative  energy.  Changes  in  the 
final  projectile  energy  arise  from  photon  absorption  and  emission.  The  laser  considered  is 
a  monotonic,  plane  polarized  CO2  laser  (hui  =  O.llTeV).  with  the  polarization  direction 
parallel  to  the  initial  projectile  velocity.  The  restricted  basis  set  -  IS,  2S,  2Po,  and  2P-ti 
-  is  used  for  the  dressing  of  the  states,  which  in  turn  provide  the  basis  set  used  within 
the  Multichannel  Eikonal  Treatment.  The  axis  of  quantization  is  along  the  direction  of 
incidence  of  the  projectile  electron,  and  as  stated  above  is  also  taken  as  the  direction  of 
polarization  of  the  laser  field. 

Cross  sections  vs  field  strength  are  presented  in  Fig. 6  and  Fig. 7,  for  Floquet  dressed 
hydrogen  atoms  whore  6 ^  =  0°.  The  cross  sections  are  compared  over  a  range  of  orders 
of  approximation  for  the  dressing  of  the  atomic  states,  it  can  be  seen  that  there  are 
very  striking  differences  between  the  different  levels  of  approximation,  especially  between 
the  first  order  cross  section  (which  fails  to  exhibit  a  second  stationary  point)  and  the 
higher  orders.  It  is  also  apparent  for  both  the  Born- Wave  approximation  (Fig. 6)  and 
the  Multichannel  Eikonal  Treatment  (Fig. 7)  that  at  least  a  third  order  approximation 
in  the  Floquet  dressing  is  needed  in  order  to  reach  convergence  over  the  range  of  field 
strengths  studied.  The  divergence  from  convergence,  for  each  order  of  approximation, 
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occurs  at  successively  higher  field  strengths.  Cross  sections  vs  field  strength  are  presented 
in  Fig. 8  and  Fig.9  for  Perturbatively  dressed  hydrogen  atoms.  In  contrast  to  the  Floquet 
dressed  atomic  states,  the  Perturbative  dressed  atomic  states  produce  cross  sections  with 
very  little  significant  difference  between  the  successive  orders  of  approximation.  Although 
convergence  is  not  even  satisfactorily  reached  by  a  fifth  order  approximation  -  a  fifth  order 
cross  section  is  nearly  identical  to  a  third  order  cross  section,  but  a  fourth  order  cross 
section  is  distinctly  different  -  a  reasonable  qualitative  result  can  be  obtained  from  a  first 
order  approximation. 

As  the  field  strength  e0  rises,  cross  sections  for  both  the  Floquet  dressed  atomic  states 
(Fig. 6  and  Fig. 7)  and  the  Perturbative  dressed  atomic  states  (Fig. 8  and  Fig.9)  show  an 
initial  rise  for  the  1S-2S  cross  sections  and  an  initial  fall  for  the  lS-2Po  cross  sections.  This 
initial  trend  is  due  to  the  dressed  atomic  states  S  and  P,  which  show  increasing  amounts  of 
each  others  zero  field  characteristics.  The  percentage  2Po  component  in  each  of  the  dressed 
atomic  states  S  and  P,  presented  in  Fig.'J  0  as  a  function  of  field  strength,  illustrates  this  in 
a  very  visual  manner  for  third  order  Floquet  dressed  and  third  order  Perturbative  dressed 
states,  (a)  and  (b)  respectively.  The  plots  of  Fig. 10  exhibit  the  play  off  between  the 
dressed  states  S  and  P,  with  each  simultaneously  giving  up  more  of  its  intrinsic  zero  field 
character,  and  adopting  more  of  the  other  character.  For  both  approaches  to  the  dressed 
atomic  states,  this  character  interchange  continues  past  the  point  where  the  dressed  state 
S  has  a  greater  2Po  characteristic  than  the  dressed  state  P.  It  is  to  be  expected  that  the 
scattering  cross  sections  would  rise  and  fall  according  to  the  character  of  the  dressed  atomic 
states,  ie,  the  1S-2S  cross  section,  is  expected  to  rise  past  the  falling  lS-2Po  cross  section 
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as  e0  increases.  This  is  the  case  for  the  cross  sections  produced  by  Perturbatively  dressed 
atomic  states,  where  the  plots  of  Fig. 8  and  Fig. 9  exhibit  a  monotonic  rise  and  fall  in  the 
1S-2S  and  1S-2P0  cross  sections.  But  it  is  clearly  not  the  case  for  cross  sections  produced 
by  Floquet  dressed  atomic  states.  The  plots  of  Fig.6  and  Fig. 7  show  cross  sections,  where 
there  are  distinct  peaks  and  dips,  which  seem  to  bear  no  relation  to  the  character  of  the 
dressed  states.  These  peaks  and  dips  arise  from  the  contributions  to  the  state-to-state 
cross  sections,  from  the  infinite  set  of  periodic  states.  The  first  order  cross  sections,  of 
Fig.6  and  Fig.7,  demonstrate  very  clearly  the  contribution  to  the  cross  section  via  the  first 
order  terms  only.  As  the  field  strength  rises,  the  contribution  to  the  cross  section  via  the 
first  order  terms  increases,  and  these  contributions,  whilst  initially  augmenting,  eventually 
begin  to  curtail  the  rising  change  in  character  sties,  thereby  causing  the  single  stationary 
point.  For  the  second  order  cross  sections,  contributions  are  allowed  from  both  first  and 
second  order  terms.  As  the  field  strength  rises,  the  rising  contribution  to  the  cross  sections 
from  second  order  terms  is  now  added  to  the  first  order  contributions.  The  second  order 
terms  initially  curtail,  but  then  augment  the  rising  change  in  characteristics,  and  eventually 
managed  to  dominate  the  first  order  contributions  to  cause  the  second  stationary  point. 
This  distinct  difference  between  the  first  and  second  order  approximations  to  the  cross 
sections,  emphasizes  the  role  played  by  the  higher  order  terms  through  the  infinite  set 
of  periodic  dressed  states.  Since  the  infinite  set  of  periodic  states  is  not  predicted  in  a 
Perturbative  dressing  of  the  atom,  contributions  from  higher  order  terms  cannot  occur. 
Thus  no  stationay  points  can  be  expected,  for  the  Perturbative  cross  sections,  as  the  cross 
sections  follow  the  character  of  the  dressed  atomic  states. 
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It  has  already  been  explained  in  section  3  that  any  experiment  would  measure  cross 
sections  over  a  range  of  phase  factors  Sw.  Cross  sections  are  presented  in  Fig.ll  and  Fig. 12, 
compared  over  a  range  of  phase  factors  from  8W  =  0°  to  90°.  Cross  sections  for  phase  factors 
in  the  range  90°  to  360°  are  a  repeat  of  those  shown.  A  model  attempting  to  predict 
experimental  results,  therefore  requires  averaging  over  the  reduced  range  of  0°  to  90° 
alone.  The  cross  sections  pertain  to  Floquet  dressed  hydrogen  atoms,  and  are  calculated  by 
the  Born- Wave  approximation  (Fig.ll)  and  the  Multichannel  Eikonal  Treatment  (Fig.12). 
Cross  sections  for  Perturbative  dressed  atomic  states  are  not  presented,  because  of  their 
phase  independence.  The  above  results  (Fig. 6  and  Fig. 7),  which  illustrate  convergence 
in  the  cross  sections  for  rising  orders  of  atomic  dressing,  justify  the  use  of  a  third  order 
approximation  in  the  following  cross  sections.  The  Born-Wave  approximation  (Fig.ll) 
and  the  Multichannel  Eikonal  Treatment  (Fig.12)  both  show  a  range  of  different  responses 
of  the  cross  section  to  the  changing  field  strength.  These  range  from  the  two  distinct 
stationary  points  exhibited  by  cross  sections  for  =  0° ,  and  the  single  stationary  point  of 
the  cross  sections  for  8^  —  90°.  For  8U  =  90°  there  are  no  contributions  to  the  cross  sections 
from  first  or  third  order  terms,  so  that  the  single  observed  stationary  point  is  dependent 
on  the  second  order  terms  only.  For  phases  between  0°  and  90°  the  contributions  from  the 
different  orders  vary,  giving  rise  to  the  range  of  cross  sections  between  the  two  extremes 
described.  Both  the  Born- Wave  and  the  Multichannel  Eikonal  results  are  presented  in 
Fig.  13  when  the  cross  sections  are  averaged  over  the  required  range  of  phases,  from  8 ^  = 
0  to  90' .  It  is  immediately  noticed,  while  the  dressed  atomic  states  S  and  P  practically 
exchange  zero  field  characteristics  as  the  electric  field  strength  rises,  that  the  averaged 
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cross  sections  for  the  1S-2S  and  lS-2Po  transitions  never  cross.  Also  the  field  strength  at 
which  the  two  cross  sections  are  identical,  is  the  same  as  the  field  strength  at  which  the 
dressed  atomic  states  S  and  P  each  contain  a  50:50  combination  of  2S  and  2Po  field  free 
states.  This  non-crossing  of  the  averaged  cross  section  in  Fig.13,  is  due  to  the  projection 
of  the  “dressed”  transition  amplitudes  onto  the  “measured”  transition  amplitudes.  If  the 
dressed  atomic  states  S  and  P  are  dominated  by  each  others  zero  field  characteristics,  they 
will  then  contribute  more  to  the  state-to-state  cross  section  of  the  character  that  dominates 
them,  via  the  infinite  set  of  periodic  states,  than  to  the  state-to-state  cross  sections  of  their 
own  zero  field  character.  Hence  a  dressed  state  S,  with  a  predominately  2Po  character, 
contributes  mostly  to  the  lS-2Po  cross  section,  rather  than  the  1S-2S  cross  section.  This 
then  leads  to  viewing  the  phase  dependent  state-to-state  cross  sections,  as  being  variations 
from  these  phase  averaged  cross  sections.  These  variations  occur  due  to  constructive  and 
destructive  interference.  When  the  cross  sections  are  averaged,  the  effect  of  the  interference 
is  cancelled  out.  To  contrast  with  this  non-crossing  of  the  cross  sections  in  Fig.13,  cross 
sections  vs  frequency  are  presented  in  Fig.  14  for  third  order  Perturbative  dressed  atomic 
states.  For  the  Perturbative  dressed  atomic  states,  the  1S-2S  and  the  lS-2Po  cross  sections 
do  cross  in  Fig. 14,  with  the  cross  sections  again  being  identical  at  the  same  field  strength 
that  the  dressed  states  S  and  P  each  contain  a  50:50  mixture  of  2S  and  2Po  field  free  states. 
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6.LASER  PERTURBATION  OF  THE  PROJECTILE 

The  full  effect  of  the  laser  field  is  now  considered  by  applying  the  laser  perturbation 
to  both  the  projectile  and  the  target  atom.  The  perturbed  Hamiltonian  is 

H  =  -kP*~£P*-A/C+2Z TWO? 

+  £i(pri+A/c)’  +  V(r)  +  V(R,r)  (49a) 

i 

=  hPR  -  •  A/c  +  +  Bp  +  V(R,  r)  (496) 

in  the  center  of  mass  frame  of  reference,  where  the  Floquet  Hamiltonian  for  the  atom 

Hr  =  Y,  J(Pw  +  A/c)2  +  F(r) 

i 

has  already  been  considered  in  section  2.  The  charge  and  the  mass  of  the  projectile  are  Z 
and  mv.  For  heavy  particle  collisions,  rnp  ^  p  the  reduced  mass  of  the  collisional  system. 
This  difference  occurs  because  the  laser  perturbation  must  be  applied  in  the  laboratory 
frame  of  reference,  and  the  resulting  Hamiltonian  is  then  transformed  to  the  center  of  mass 
frame  of  reference.  The  laser  perturbations  can  be  handled  by  treating  the  projectile  and 
the  atom  as  isolated  systems.  The  internal  atomic  portion  has  already  been  considered  via 
the  solution  of  idi<p  —  Hf4>  in  section  2.  Thus  a  solution  for  the  projectile  Hamiltonian 


HP  = 


1  p2 
2  “ 


m  p 


pR  '  AA  +  Kj(AZA) 


(50) 


only  is  needed.  The  wave  function  which  satifies  the  Schrodinger  equation  ?’4^  =  Hp$, 
is  represented  by 


$  =  exp 


i^P#-R+^  J  PR  •  (A0/c)cos(u;t  +  6u..)  -  a  j  dr 


(51) 
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where  a  is  a  constant  phase  shift  that  cancels  out  in  the  collisional  S-matrix.  In  the  original 
paper  by  Volkov  (1935)  the  solution  $  is  expressed  in  the  general  form 


$  =  exp  Pp  ■  R.+  ^-P r  ■  ( A0 / u>c)  sin(wf + ^ )  -  a  j 


•Za 


(52) 


where  Z0  is  dependent  on  P#,  A0,  w,  mp,  8W  and  x.  If  it  is  assumed  that  the  only 
dependence  of  Z0  on  x  can  be  acknowledged  by  varying  8W>  then  8W  can  be  easily  varied 
for  any  value  of  x,  to  give 


u>x  +  8W  =  0  and  Z0  —  1 


for  use  in  eq.(52)  giving  rise  to  the  usual  Volkov  dressed  wave  function  for  a  laser  perturbed 
plane  wave 


In  deriving  the  wave  function  eq.(53),  the  assumption  has  been  made  that  there  are  no 
transient  effects  arising  from  turn-on  or  turn-off  of  the  laser,  as  is  reasonable  in  this  soft- 
photon  weak-field  regime.  Having  set  ujx  4-  8 u  and  Z0~  1,  it  must  be  realized  that  8 w  has 
now  taken  on  a  definite  and  distinct  value  for  each  collision.  If  /  =  x  is  defined  as  being 
the  time  at  which  the  point  of  closest  approach  for  the  collision  is  reached,  then  8^  is  the 
phase  shift  of  the  laser  field  at  the  point  of  closest  approach,  as  shown  in  Fig.l .-  In  any 
experiment,  collisions  will  occur  across  a  wide  range  of  space  and  time.  Thus  8 w  will  not 
be  in  general  a  constant  for  all  collisions.  As  mentioned  in  section  3,  a  theoretical  model 
that  purports  to  predict  experimental  results,  must  include  an  average  of  the  collisions 
over  a  range  of  8 w  from  0°  to  360°. 
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This  dressed  projectile  wave  function  can  now  be  combined  with  the  dressed  atomic 
wave  functions  from  section  2,  and  using  the  techniques  of  section  3  the  laser  perturbed 
collisional  S-matrix 

Spq{i)  =  ($p| $q)6pq  +  i2ir  Tpq  6(QP  +  +  (n+x-m)hw  -  Qq  -  (54) 

can  be  derived.  This  S-matrix  is  of  exactly  the  same  form  as  eq.(33),  with  the  new  “dressed” 
transition  amplitude 

=  £  E  £  E  £  a?daw  (*»  imi>  urn  (55) 

n  m  a  fi  x 

The  Bessel  function  argument  D  is  ^(k,  —  k/)  •  (A 0/u?c),  and  as  usual  arises  from  a 
Fourier  expansion  of  the  Volkov  dressing  of  the  projectile.  The  new  conservation  of  energy 
requirement  is  now 

5Zk2f  =  Ytkl-(n  +  x-m)hu}  (5fi) 

where  the  projectile  contribution  x,  to  the  number  of  photons  absorbed  and  emitted  during 
transitions  between  the  dressed  wave  functions,  is  also  included.  The  “measured”  transition 
amplitude  T# Q  can  now  be  written  in  terms  of  this  new  “dressed”  transition  amplitude  Tpg, 
as  outlined  in  section  3. 

Cross  sections  vs  field  strength  are  presented  in  Fig. 15  and  Fig. 16,  for  third  order 
Floquet  dressed  hydrogen  atoms,  in  a  Born- Wave  approximation  of  laser  assisted  e~- 
H(1S)  collisions.  Cross  sections  with  the  laser  perturbation  included  for  both  the  atom 
and  the  projectile,  are  compared  with  cross  sections  that  neglect  the  laser  perturbation 
of  the  projectile.  It  has  already  been  shown  that  the  Volkov  dressing  produces  a  very 
significant  effect  on  individual  differentia]  and  integral  cross  sections,  with  a  specific  final 
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projectile  energy  (Weinsgarthofer  et  al  1979,  Byron  et  al  1987,  Jetzke  et  al  1987,  Sarkar 
and  Chakraborty  1988).  When  the  state-to-state  transition  cross  sections  are  however 
summed  over  an  increasingly  wider  range  of  projectile  energies,  corresponding  to  higher 
order  multiphoton  events,  the  fully  perturbed  collisional  cross  sections  converge  on  the  cross 
sections  that  include  a  laser  perturbation  of  the  atom  only.  This  is  clearly  evident  from 
the  cross  sections  displayed  in  Fig. 15  and  Fig. 16.  As  would  be  expected,  this  convergence 
improves  as  the  initial  impact  energy  increases.  This  study  has  only  been  conducted 
using  the  Born- Wave  approximation,  due  to  the  difficulties  associated  with  producing 
a  practical  numerical  technique  that  can  handle  the  Volkov  dressed  states  within  the 
Multichannel  Eikonal  Prescription  in  a  reasonable  amount  of  computer  time.  Despite 
the  practical  problems  there  is  no  difficulty  in  principle  from  using  the  new  transition 
amplitude,  and  there  is  nothing  to  indicate  that  the  convergence  exhibited  by  the  Born- 
Wave  approximation  would  not  be  reproduced  by  the  Multichannel  Eikonal  Treatment  as 
well. 
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7.CONCLUSION 


This  work  has  studied  laser  assisted  collisions,  in  a  soft-photon  weak-field  regime, 
for  intermediate  energy  projectiles.  Within  this  regime,  the  most  appropriate  method 
of  incorporating  the  laser  perturbation  of  the  atom,  has  been  found  to  be  the  Floquet 
prescription.  This  Floquet  approach,  with  its  predicted  infinite  set  of  periodic  states, 
is  considered,  by  the  authors,  to  be  superior  to  the  more  traditional  time-dependent 
Perturbation  theory  with  which  it  was  compared,  because  it  provides  a  very  concise 
description  of  projectile-atom  collisions  in  a  laser  field.  Unlike  the  time- dependent 
Perturbation  theory,  the  Floquet  prescription  exactly  describes  the  location  during  the 
collisional  sequence  of  events  at  which  photons  are  absorbed  and  emitted.  From  section 
4,  there  are  two  distinct  photon  absorption  and  emission  processes,  allowed  within  the 
Floquet  approach.  The  first  being  absorption  and  emission  by  the  projectile,  during  the 
projectile-atom  interaction  (e,A,7iu;),  and  the  second  being  absorption  and  emission  by 
the  atom,  well  before  and  well  after  the  projectile- atom  interaction,  (A,fuo).  The  laser 
perturbation  (e,fiu>)  of  the  projectile  was  also  acknowledged  in  this  work,  in  the  form  of 
Volkov  dressed  projectile  wave  functions.  It  has  however  been  shown  that  these  projectile 
states  have  very  little  effect  on  the  state-to-state  cross  sections,  when  summed  over  all 
possible  final  projectile  energies  consistant  with  a  given  projectile  energy.  These  dressed 
wave  functions  also  gave  no  indication  of  allowing  photon  absorption  or  emission  by  the 
projectile,  before  and  after  the  projectile-atom  interaction. 

This  work  ha:-,  also  shown  that  the  Floquet  treatment  provides  cross  sections  with  some 
very  distinct,  interesting  and  new  features,  the  most  important  of  which  are  summarized 
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below: 


(a) :  When  there  is  significant  coupling  between  the  atomic  states,  the  Floquet 
approach  requires  at  least  a  third  order  approximation  in  the  dressed  atomic  states,  to 
attain  convergence  in  the  state-to-state  cross  sections  over  the  full  range  of  field  strengths 
studied  in  this  work.  This  convergence  in  the  higher  order  dressed  states  is  critical,  since  the 
first  order  approximation  does  not  even  yield  a  qualitative  indication  of  the  overall  variation 
of  the  converged  cross  sections  with  field  strength,  and  a  second  order  approximation  is 
often  not  much  better.  This  indicates  that  the  contribution  to  the  cross  sections  from  the 
higher  terms  in  the  atomic  dressed  states  are  important. 

(b) :  The  Floquet  cross  sections  display  a  laser  phase  dependence  of  6^.  It  has  also 
been  shown  that  for  comparison  with  experiment,  it  is  necessary  to  take  an  average  of  the 
state-to-state  cross  sections  over  this  phase,  from  6^  =  0°  to  360°. 

(c) :  For  a  third  order  approximation  the  Floquet  dressed  atomic  states  can  be  written 
in  the  general  form. 

| S)  =  A\2S,n)  +  Bei6v  |2Po,n  +  l)  +  Be~'6*  \2P0,v  -1) 

+  Cet26“\2S,n+2)  +  Ce~t26u\2S,n-2) 

+  DelU"  |2P0 ,  n  +  3)  +  |2P0 ,  n  - 3)  (57) 

and 

| P)  =  A\2Pb,n)  +  B(,b^  !25,u  +  l'  +  j2S,?j-l) 

-r  |2P0,«+2>  +  Ct~l2b"  j2Po,n-2) 

•f  Dtl36“  |25,n  -f  3)  ■+  De~l36*  |2S,n-3)  (57) 
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If  these  dressed  atomic  states  S  and  P  are  plotted  against  field  strength  then  the 
dressed  states  practically  exchange  intrinsic  zero  field  characteristics  as  the  field  strength 
rises.  Eventually  the  dressed  state  S  will  have  a  greater  2Po  character  than  the  dressed  state 
P.  Despite  this  character  interchange  of  the  dressed  states,  the  1S-2S  state-to-state  cross 
section  never  rises  aoove  the  lS-2Po  state-to-state  cross  section  when  the  phase  averaged 
1S-2S  and  lS-2Po  cross  sections  are  plotted  against  field  strength.  This  is  demonstrated 
in  Fig.17  for  electron  hydrogen  collisions  using  Floquet  dressed  atomic  states.  It  should 
be  noticed  that  the  1S-2S  and  lS-2Po  cross  sections  have  the  same  value,  at  precisely  the 
same  electric  field  strength  that  the  dressed  states  S  and  P  have  a  50:50  composition  of  2S 
and  2Po  field  free  states. 

Where  ever  possible,  numerical  calculations  are  based  on  both  the  Multichannel 
Eikonal  Treatment  and  the  Born- Wave  approximation.  Both  treatments  produce  similar 
qualitative  variations  to  the  cross  sections  with  electric  field  strength,  and  any  quantitative 
difference  tends  to  orginate  from  the  difference  between  the  field  free  scattering  theories 
adopted. 
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FIGURE  AND  TABLE  CAPTIONS 


Figure  (1):  Vector  potentials  _  A0cos (u>i  +  Sw)  and . A0cosu ;t,  for  a 

monotonic  plane  polarized  laser  field  of  frequency  u>,  positioned  in  time  relative  to  the 
orbit’s  pericenter  at  time  t  =  0. 

Figure  (2):  A  restricted  two  state  Floquet  matrix  for  the  2S  and  2Po  states  of  a  hydrogen 
atom  in  a  laser  field,  is  shown  for  Floquet  N-  n=  -1,  0,  and  1,  where  an  unrestricted  matrix 
would  have  Floquet  N-  up  to  ±oc.  The  matrix  is  diagonalized  to  determine  the  “Quasi” 
energies  Qq  and  the  dressed  atomic  states 

Figure  (3  a-c):  Radial  distribution  functions  of  the  dressed  atomic  states  S  and  P,  for 
a  third  order  Floquet  dressed  hydrogen  atom,  are  compared  for  increasing  electric  field 
strengths  of  a  monotonic  CO2  laser  field  of  photon  energy  hu  =  0.0043au. 

Figure  (4):  The  percentage  2Po  component  vs  field  strength,  of  the  dressed  atomic  states 

S -  and  P _ ,  are  compared  for  a  hydrogen  atom  in  a  laser  field  of  photon 

energy  Tiu>  =  0.0043au. 

Figure  (5):  Energy  level  diagram,  E  —  Eo  +  hu>,  of  the  hydrogen  dressed  atomic  states  S 
and  P  at  a  photon  energy  of  hu>  =  0.0043au. 

Figure  (6):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Pu  transitions,  (a)  and  (b) 
respectively,  in  e~-H(lS)  collisions  in  a  laser  field,  using  Floquet  dressed  hydrogen  atoms 
and  a  Born-Wave  approximation,  are  compared  over  a  range  of  orders  of  approximation 
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for  the  Floquet  dressing  ( _ =  1st, .  =  2nd, _ =  3rd  and _ 

=  4th  order)  at  a  photon  energy  huj  =  0.0043au  and  =  0°.  The  laser  perturbation  is 
included  for  the  atom  only.. 

Figure  (7):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  e~-H(lS)  collisions  in  a  laser  field,  using  Floquet  dressed  hydrogen 
atoms  and  a  Multichannel  Eikonal  Treatment,  are  compared  over  a  range  of  orders  of 

approximation  for  the  Floquet  dressing  ( _ =  1st, .  =  2nd, _ =  3rd 

and _  =  4th  order)  at  a  photon  energy  =  0.0043au  and  =0°.  The  laser 

perturbation  is  included  for  the  atom  only. 

Figure  (8):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Po  transitions,  (a) 
and  (b)  respectively,  in  e~-H(lS)  collisions  in  a  laser  field,  using  Perturbative  dressed 
hydrogen  atoms  and  a  Born- Wave  approximation,  are  compared  over  a  range  of  orders  of 

approximation  for  the  Floquet  dressing  ( _ =  1st, _ =  2nd, .  = 

3rd, _ —  4th  and  _  =  5th  order)  at  a  photon  energy  hui  =  0.0043au.  The 

laser  perturbation  is  included  for  the  atom  only. 

Figure  (9):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e_-H(lS)  collisions  using  Perturbative  dressed  hydrogen 
atoms  and  a  Multichannel  Eikonal  Treatment,  are  compared  over  a  range  of  orders  of 

approximation  for  the  Floquet  dressing  ( _ =  1st,  =  2nd, .  = 

3rd, _ _ =  4th  and _ =  5th  order)  at  a  photon  energy  —  0.0043au,  The 

laser  perturbation  is  included  for  the  atom  only. 
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Figure  (10):  The  percentage  2P<>  component  vs  field  strength,  for  the  dressed  states 


S  _  and  P - ,  are  compared  for  both  Floquet  and  Perturbative  dressed 

hydrogen  atoms,  (a)  and  (b)  respectively,  in  a  laser  field  of  photon  energy  hu  =  0.0043au. 

Figure  (11):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e~-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  and  a  Born- Wave  approximation,  are  compared  over  a  range  of  phases  8W 

( _  =0°, _ =  22.5°, .  =  45°, _ =  67.5°  and _ = 

90°)  at  a  photon  energy  hu  =  0.0043au.  The  laser  perturbation  is  included  for  the  atom 
only. 

Figure  (12):  Cross  sections  vs  field  strength  for  the  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e~-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  and  a  Multichannel  Eikonal  Treatment,  are  compared  over  a  range  of 

phases  8W  ( _ -  0°, _ =  22.5°, .  =  45°, _ =  67.5°  and _ 

- =  90°)  at  a  photon  energy  huo  —  0.0043au.  The  laser  perturbation  is  included  for 

the  atom  only. 

Figure  (13):  Cross  sections  vs  field  strength  for  the  1S-2S _ and  lS-2Po _ 

—  transitions  in  laser  assisted  e“-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  are  compared  in  the  Born-Wave  approximation  and  the  Multichannel 
Eikonal  Treatment,  (a)  and  (b)  respectively,  at  a  photon  energy  hue  =  0.0043au.  The  cross 
sections  are  an  average  of  phase  dependent  cross  sections  for  —  0‘  to  90%  and  the  la^er 
perturbation  is  included  for  the  atom  only. 
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Figure  (14):  Cross  sections  vs  field  strength  for  the  1S-2S _  and  lS-2Po _ 

_  transitions  in  laser  assisted  e--H(lS)  collisions  using  third  order  Perturbative  dressed 
hydrogen  atoms  are  compared  in  the  Born- Wave  approximation  and  the  Multichannel 
Eikonal  Treatment,  (a)  and  (b)  respectively,  at  a  photon  energy  hu>  —  0.0043au.  The  laser 
perturbation  is  included  for  the  atom  only. 

Figure  (15):'  Cross  sections  vs  field  strength  for  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e~-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  and  a  Born- Wave  approximation  are  compared  with  the  laser  perturbation 

included  for  the  atom  only  ( _  )  and  the  laser  perturbation  included  for  both  the 

atom  and  the  projectile,  where  the  later  cross  sections  are  a  summation  over  a  range  of 

final  projectile  energies  ( _ _ =  ±20/)u>, .  =  ±25/iu>  and _ =  ±30 hun 

about  the  field-free  final  projectile  energy).  The  phase  6 ^  —  0°  and  the  photon  energy 
hu>  =  0.0043au. 

Figure  (16):.  Cross  sections  vs  field  strength  for  1S-2S  and  lS-2Po  transitions,  (a)  and 
(b)  respectively,  in  laser  assisted  e_-H(lS)  collisions  using  third  order  Floquet  dressed 
hydrogen  atoms  and  a  Born- Wave  approximation  are  compared  with  the  laser  perturbation 

included  for  the  atom  only  ( _  )  and  the  laser  perturbation  included  for  both  the 

atom  and  the  projectile,  where  the  later  cross  sections  are  a  summation  over  a  range  of 

final  projectile  energies  ( —  ±20/iu;,  =  z.25hu:  and _ _  =  ±30 ha: 

about  the  field-free  final  projectile  energy).  The  phase  6^.  =  90f  and  the  photon  energy 
hu>  —  0.0043au. 
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Figure  (17):  Cross  sections  in  e~-H(lS)  collisions  are  compared  against  the  2Po  component 
of  the  dressed  atomic  states,  (a)  and  (b)  respectively,  over  a  range  of  electric  field  strengths 
for  third  order  Floquet  dressed  hydrogen  atoms  in  a  laser  field.  The  cross  sections  1S-2S 

_  and  lS-2Po _ are  an  average  of  the  phase  dependent  cross  sections  over 

a  range  of  phases  from  8W  =  0°  to  90°,  using  a  Multichannel  Eikonal  Ttreatment.  For  the 

percentage  2P0  component _  represent  the  dressed  state  S  and _ represent 

the  dressed  state  P. 

Table  (1):  “Dressed”  transition  amplitude  elements  Jpg,  final  projectile  energies 
and  coefficients  of  the  “dressed”  transition  amplitude  Jps  within  the  “measured”  transition 
amplitude  JiS2S- 
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Table  (1):  “Dressed’'  transition  amplitude  elements  Tpq.  final  projectile  energies 
and  coefficients  of  the  “dressed”  transition  amplitude  Tpq  within  the  “measured”  transition 
amplitude  T1S2S- 


“Dressed” 

Final 

Projection  operator 

transition  amplitude 

projectile 

probability  amplitude 

matrix  elements  Tpq 

energy 

from  Tpq  onto  Txsjs 

Tn  =  AV\s2s 

£*/  =  frls 

A 

T«  =  BVis,P0e-iS“ 

&kf  =  &*2S  + 

A 

T<  1  =  BVls,P0e+it* 

&*/  =  &*a5  - 

A 

T)i  =  AVisjft 

2^*7 = 

Be+i6* 

r3i  =  BVisase-*- 

2^/  =  £*22*,  +  *w 

Be+li * 

Tji  =  BVlS2Se+iS~ 

Be+l6- 

T71  =  AV\S2Pq 

Be~'6* 

Tn  =  BVXSiSe 

jfr*/aS&*2ft+,lW 
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Figure(l):  Vector  potentials  -  A0cos(»*t 

monotonic  plane  polarized  laser  field  of  frequency 
orbit's  pericenter  at  time  t  =  0.- 


and  A0cosu.-t.  for  a 

.  positioned  in  time  relative  to  the 


FIGURE  1:  SMITH  AND  FLANNERY 
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Figure(2):  A  restricted  two  state  Floquet  matrix  fur  the  2S  and  2P(,  states  of  a  hydrogen 
atom  in  a  laser  field,  is  shown  for  Floquet  ,\-  n=  -1.0.  and  1.  where  an  unrestricted  matrix 
would  have  Floquet  N-  up  to  -oc.  The  matrix  is  diagonalized  to  determine  the  ‘’Quasi" 
energies  Qq  and  the  dressed  atomic  states  . 
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Figure(4):  The  percentage  2P0  component  vs  field  strength,  of  the  dressed  atomic  states 

S - —  and  P - -  are  compared  for  a  hydrogen  atom  in  a  laser  field  of  photon 

energy  /?«;  =  0.0043au. 
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Transport- Collisional  Master  Equations  for  Termolecular  Recombination 

M.  R.  Flannery 
School  of  Physics 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332-0430 

Abstract 

Sets  of  Transport-Collisional  Master  Equations  are  developed  for  the  microscopic  distribution  n(R ,  E ,  L ) 
of  pairs  over  internal  separation  iZ,  energy  E  and  orbital  angular  momentum  L  of  (A-B)  pairs  in  a  background 
gas  M  of  variable  density.  Expressions  are  also  provided  for  the  rate  of  recombination  of  A  and  B  as  a 
function  of  gas  density.  Analytical  solutions  for  the  pair  distributions  n  and  microscopic  probabilities  for 
recombination  are  obtained  in  the  classical  absorption  limit.  They  pertain  to  exact  (A-B)  trajectories  under 
general  symmetric  interaction  V^fl)  between  A  and  B  and  are  applied  to  ion-ion  and  electron-ion  collisional 
recombination  in  a  gas.  A  classical  variational  method  is  also  presented. 


1*  Introduction 

In  previous  papers1'4,  the  rate  a  for  termolecular  recombination, 

Or 

A  +  B  +  M  Zl  AB  +  M  (1.1) 

kd 

of  species  A  and  B  (assumed  structureless)  in  a  low  density  gas  M  was  considered  by  exact  (quasi¬ 
steady-state1  and  variational2)  treatments  and  by  various  approximate  methods  based  on  energy-diffusion8, 
bottleneck1  and  electrical  analogues4.  The  exact  treatments1’2  involve  the  solution  of  either  a  one¬ 
dimensional  set  of  integral  equations  for  the  relaxation  of  the  (A-B)  pair  distribution  n(E)  in  relative  energy 
E  or  a  variational  determination  of  n(E)  with  respect  to  external  parameters  varied  to  yield  the  minimum 
rate  a.  Both  procedures1,2  are  equivalent  to  a  quasi-steady-state  distribution  n(E)  among  highly  excited 
levels  E.  Ion-atom  association5  at  low  gas  density  involves  the  solution  of  a  two-dimensional  set  of  integral 
equations  in  E  and  in  relative  (orbital)  angular  momentum  L  of  the  associating  (A-B)  pair. 

The  aim  of  this  paper  is  to  develop  the  appropriate  generalization  to  all  gas  densities,  when  there  is 
non- equilibrium  not  only  in  the  E  and  L  degrees  of  freedom  but  also  in  the  (A-B)  internal  separation  R.This 
inhibited  relaxation  arises  from  the  increased  difficulty  of  transport  of  A  towards  B  by  diffusional-drift 
through  the  gas  M.  The  two  particle  phase-space  distribution  n(R,  p;f)  of  (A-B)  Tiair:  with  respect  to  R 
and  relative  (orbital)  momentum  p  =  is  then  governed  by  the  microscopic  transport  equation6, 

d  dn 

— n(R,p;*)=  —  +  v-  VRn-  VRF  •  Vpn  (1.2a) 

dn 

=  ¥  +  [VH-(nv))p-VRy.i'Tpn]R  (1.26) 

where  ^(R)  is  the  energy  of  mutual  interaction  between  A  and  B  with  reduced  mass  This  transport 

rate  is  then  set  equal  to  the  input-output  rate  for  the  net  collisional  production  at  fixed  R  of  state  p  from 
all  states  p'  of  the  pair  (A-B).  A  Boltzmann  equation  for  n(R,p;f)  is  then  obtained6  and  is  valid  for  dilute 
concentrations  of  reactants  A  and  B  in  a  gas  bath  M  of  general  gas  density  N  so  that  (A-B)  collisions  with 
the  gas  M  provide  the  dominant  state  changing  mechanism.  Naumann7  has  recently  reproduced  the  same 
basic  equation6  via  an  alternative  stochastic  route. 
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The  next  development8  was  the  transformation  of  (1.2)  for  symmetric  F(R)  to  the  more  natural  form, 


-n(R,J EuL^t) 


dnt  1  d 


[R2n(R,  Ei,  L*;i)vR]B"L2 


(1.3) 


for  the  resulting  R-symmetric  (A-B)  pair  distribution  nt(R)  =  n(R,Et,I2)  over  R,  Et  and  I?.  The 
radial  speed  dR/dt  is  vr  so  that  the  microscopic  radial  current  of  pairs  expanding  and  contracting  at  R  is 
jt  =  R2ntVR  per  unit  interval  dEtdL*>  which  is  conserved  and  which  depends  only  on  Et  under  equilibrium 
in  R,  Ex  and  If  (Appendix  A).  The  collisional  RHS  of  the  Boltzmann  equation6  may  now  be  written  in  the 
collisional  input-output  form, 


In(R,  £„!?;<)  =  =  -JdEj  j  dL)  [n^R^R)  -  nf{R)v}l{R)]  (1.4) 

where  i /tj  is  the  frequency  per  unit  dEfdL 2  of  state  changing  EULX  —►  E/,I/  transitions  produced  by 
collision  between  the  (A-B)  pairs  in  state  i  and  the  gas  species  M  and  the  integrations  are  over  all  Ej  and 
1/  accessible  at  a  given  R. 

The  above  two  developments6,8  of  (1.2)  and  (1.3)  now  facilitate  the  constn  ~tion  in  this  paper  of  two 
sets  of  Transport-Collisional  Master  Equations  for  the  pair  microscopic  distribution  n(R,  E)  I2)  and  the 
I2-integrated  distribution  n(R,  E).  The  sets  essentially  represent  three-dimensional  (R,  E,  I2)  and  two- 
dimensional  (R,  E)  integral  equations  respectively.  The  two-dimensional  (R,  E)  set  is  however  composed 
of  two  separate  sets,  each  valid  in  exclusive  regions  of  (R,  E)  space  which  become  coupled  via  boundary 
conditions  at  the  (R,  R)-interface.  The  full  sets  are  in  general  difficult  to  solve  by  customary  numerical 
algorithms. 

In  this  paper  analytical  solutions  are  obtained  in  the  absorption  or  classical  limit,  when  the  back-coupling 
terms  Ufx  in  (1.4)  are  consistently  neglected.  Analytical  expressions  are  then  derived  for  the  probability 
Pa(R0)  E,  I2)  for  association  of  (£,  I2)-pairs  with  internal  separations  R  <  R0  and  for  the  various  I2- 
and  R-averaged  probabilities,  PA(R0yE)  and  R^(Ro)  respectively,  as  a  function  of  gas  density.  Not  only 
does  this  analytical  route  provide  further  detailed  insight  into  recombination  but  also  the  strategy  therein 
suggest  procedures  valuable  for  eventual  detailed  numerical  solution  of  the  full  Transport-Collisional  Master 
Equations. 

In  this  paper  the  recombining  species  A  and  B  are  assumed  structureless  and  the  recombination  proceeds 
via  collisional  and  transport  relaxation  in  the  relative  (A-B)  coordinates  R,  E  and  I2  in  the  (atomic  or 
molecular)  gas.  Even  with  this  simplification  the  resulting  set  of  three-dimensional  transport-collisional 
equations  are  quite  complicated.  Generalization  is  however  required  to  cover  the  possibility  of  recombination 
process  to  cover  the  possibility  of  recombination  process 


Xe+  -f  Cl~  +  -Ye  -  XeCl  +  2Xe  (1.5) 

where  the  new  two-body  tidal  mechanism  introduced  by  Bates  and  Morgan10, 

{Xe+{i  hlJt)-cr}Ei  { Xet{vf,Jf)-Cr}E ,  XeCU  Xe  (1.6) 

is  mediated  by  collisions  with  the  third  bodies  M  and  was  shown  by  their  recent  computer  simulation10 
to  be  necessary  in  order  to  explain  the  observed11  production  of  XeCl  rather  than  Ae2C7  as  implied  by 
(1.1).  Although  the  ‘tw'o-body’  rate  of  (1.6)  involving  orbital  — »  internal  energy  conversion  can  in  principle 
be  added  directly  to  the  collisional  RHS  of  (1.4),  the  interaction  F(R)  between  -  CP  is  however 
orientation  dependent  and  the  transport  side  (1.3)  must  be  so  modified.  Inclusion  of  this  aspect  within  the 
present  theory  is  feasible,  but  implementation  involves  an  additional  order  of  complexity  to  what  is  already 
a  multilayered  structure.  The  detailed  inclusion  of  orientation  dependent  (A-B)  interactions  1^(R)  will  be 
considered  at  a  later  stage. 

Since  the  time-dependent  set  of  equations  based  on  (1.3)  and  (1.4)  can  be  reduced  to  a  time-independent 
set  via  a  quasi-steady  state  (QSS)  criterion  for  highly  excited  bound  (A-B)  pairs,  some  background  discussion 
of  QSS  is  useful. 
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1.1  Quasi-Steady-State  (QS5)  Approximation  : 

Association  (1.1)  which  proceeds  at  forward  rate  a(cm3$-1)  is  coupled  with  the  reverse  process  of 
dissociation  which  occurs  at  frequency  fc^s”1).  Both  processes  are  naturally  time  dependent  and  the 
measured  rate  constants  a  and  kd  satisfy 

dN-Adfl)  =  -<*NA(i)NB(i)  +  kdNAB{i)  (1.7) 

where  Na,  Nb  and  Nab  are  the  time  varying  concentrations  (cm-3)  of  reactant  and  product  species.  The 
rate  constants  a  and  kd  can  however  be  determined  by  a  time  independent  procedure1,9  in  which  association 
and  dissociation  can  be  treated  separately. 

Association  emerges1’0  in  the  time  independent  picture  as  if  equilibrium  concentrations  Na  and  Nb  of 
the  dissociated  species  A  and  B  with  relative  energies  E  in  the  range  0  <  E,  <  oo  (which  defines  a  reactant 
block  C)  are  collisionally  transferred  via  an  intermediate  block  £  of  highly  excited  states  in  the  energy  range 
0>£/>“Stoa  product  block  S  of  fully  stabilized  molecular  levels  in  the  energy  range  -S  >  E  >  -D 
where  -D  is  the  lowest  energy  level  of  the  molecule  AB  (cf.  Fig.  1  of  ref.  1).  The  energy  level  -S  lies  within1 
at  most  20 kT  below  the  dissociation  limit.  The  intermediate  block  £  is  in  quasi-steady  state  (QSS)  since 
collisions  occur  in  a  time  scale  much  shorter  than  the  characteristic  time  for  overall  recombination  and  the 
stabilized  product  block  S  is  considered  to  have  zero  population  in  the  time-independent  treatment.  The  C 
block  is  depleted  and  the  S  block  is  filled  at  a  rate  ckNaNb  via  the  steady  state  conduit  of  £-block  levels. 

Dissociation  emerges1,0  in  the  time-independent  picture  as  if  a  thermal  concentration  Nab  of  pairs  in 
the  S  block  are  collisionally  transferred  through  block  £  to  block  C  maintained  at  zero  population.  The  rate 
constants  so  deduced  satisfy  the  detailed  balance  relation, 

ocNaNb  =  kdNAB •  (1*8) 

The  QSS  procedure  is  equivalent  to  seeking  the  lowest  eigenvalue  Ao  =  kd  in  the  expansion  of  the  time- 
dependent  distribution  n,(t)  of  level  i  in  terms  of  the  basis  set  with  j=l,  2,  ...  which  define  various 

timescales  of  relaxation.  After  initial  transients  described  by  the  rapid  time  decay  dependencies  exp(-Aif), 
exp(-A20>  etc*  the  internal  degrees  of  freedom  then  quickly  relax  into  a  quasi-steady-state  which  is  slowly 
decaying  as  exp(— fc^t).  It  is  this  final  relaxation  which  is  of  interest  here. 


2.1  Basic  Theory:  Exact  Rate  in  Collisional  and  Transport  Representations 

On  assuming  that  (A-B)  pairs  in  the  block  £  of  highly  excited  bound  levels  are  in  quasi-steady-state, 
the  rate  of  collisional  termolecular  recombination  is  the  collisional  net  current, 


olNaNb 


,00  tLl  /  rE  rL),n(R)  \ 

=  dEj  dL\  /  dR  (  /  dEf  dl)  [n,(R)i/t/(J?)  -  n,(R)i//t(R)]  (2.1) 

Je  Jo  Jr;  \Jv(r )  Jo  ) 


'which  is  constant  across  an  arbitrary  level  E  =  -  |  E  |  embedded  in  block  £.  When  E  =  0,  (2.1)  is 
simply  the  net  collisional  rate  of  production  of  bound  pairs,  or  depletion  rate  (1.8)  of  dissociated  pairs  and  is 
appropriate  for  pure  (A-B)  Coulombic  attraction  V(R)  =  -e2/i?  which  does  not  support  any  bound  levels 
with  positive  energies. 

Here  n,(R)  =  n(R,  EUL*)  and  dL7x  dR  is  the  number  density  of  (A-B)  pairs  with  internal 

relative  energy  El)  internal  relative  angular  momentum  squared  if,  and  internal  separation  R  in  the 
interval  dEx  dlf  dR  about  (It,lf,R).  The  frequency  of  collisions  with  third  bodies  for  (EX)L7  — ► 
Ef  +dEf)L)  ■f  (iZj)  transitions  at  fixed  (A~H)  separation  R.  is  (R^dE f dl^ .  Also  denote  the  orbit’s 
pericenter  (— )  and  apocenter  (+)  or  turning  points  of  radial  motion  under  the  effective  interaction 


Vt(R)  =  V(R)  + 


£ 

2A 1abR2 


(2.2) 


The  maximum  value  Ltm  of  angular  momentum  L  =  R  x  p,  accessible  at  fixed  Et  and  R  is  Rp,  and  is 
determined  also  from  E,  =  V,(R)  so  that 
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L\m{R)  =  2 MABR2[E,  -  V(R))  =  R2p f  =  M2ABR2vf(R),  (2.3) 

where  pt  is  the  momentum  for  relative  (A-B)  motion- with  speed  vt.  The  maximum  L 2  accessible  overall  at 
given  Et  is  either  infinite  for  dissociated  pairs  ( Et  >  0)  or  is  lf0  given  by  (2.3)  evaluated  at  the  radius  Rq 
for  bound  ( Ex  <  0)  circular  orbits  (where  dVt/dR  =  0).  In  thermodynamic  equilibrium  at  temperature  T 
the  distribution  of  pairs  over  R,  Et  and  if  is  (Appendix  A), 


nt(R)  =»i(R,J5Mif)  = 


NANBe-B''kT  2tt 


™  -  (2irMABkT)W  RHr 

where  the  radial  speed  of  relative  (A-B)  motion  is  determined  from  the  energy  relation, 

E,  =  \mabv\  +  F(£)  +  JL-  =  \mabv2(R)  +  V(R).  (2.5) 

The  one-way  equilibrium  collisional  rates, 

Clf(R)  =  ni(R)utf(R)  =  nf(R)vft(R)  =  Cft(R)  (2.6) 

satisfy  the  principle  of  detailed  balance.  The  recombination  rate  (2.1)  then  vanishes  under  equilibrium 
conditions.  The  rate  (2.1)  with  E  =  0  may  now  be  rewritten  in  terms  of  the  fractional  distribution, 

Pi(R)  =  nt(R)/ni(R),  (2.7) 

which  represents  the  departure  of  i-pairs  from  thermodynamic  equilibrium,  as 

aNANB  =  r  dE,  r  dL f  [°°  dR  [°  dEf  /*  *  dL)  [p,(fl)  -  pf{R)}  Cif(R).  (2.8) 

Jo  Jo  Jr;  Jv(r)  Jo 

This  rate  may  alternatively  be  rewitten  in  terms  of, 

=  (29) 

introduced  here  as  the  probability  for  subsequent  collisional  stabilization  of  an  f-pair  colhsionally  produced 
in  a  bound  level  f  =  (Ej,  L2)  from  a  dissociated  i-pair.  Then, 

rOO  *  oo  »oo  p0  rL‘fm 

aNANB  =  /  dE,  /  dL2  dR  dE}  /  n,(R)u,}(R)Pf,{R)dL) ,  (2.10a) 

Jo  Jo  Jr;  Jv(R)  Jo 

in  the  collisional  representation ,  which  in  turn  can  be  rewritten  as, 

~  -  r°°  n  rR{)  r°  f1**'"  c 

aNANB  =  lim  /  dE,  /  dL 2  /  dR  /  dEj  /  n,{R)v,f(R)P?t{R)dL2f,  (2.106) 

Ro-**>Jo  Jo  Jr;  Jv{r)  Jo 

where  p20  =  2MAB(Et  -  V(Ro))>  Provided  Ro  is  large  enough  so  that  equilibrium  in  (Ex,  if)  is  collisionally 
maintained  then  the  distribution  pt(R)  =  Pj{R)  =  p{R)  for  R  >  Ro  is  independent  of  (EUL2)  and  Pjx  of 
(2,9)  then  vanishes  for  R  >  Rq,  As  Rq  — >  oo,  then  p(R)  — ♦  p^,  the  normalized  distribution  maintained  to 
preserve  steady-state. 

From  (1.3)  and  (1.4)  the  microscopic  steady-state  distributions  satisfy  the  collisional-transport  set  of 
Master  equations8 


I  *'*  =  -  f~R  dEf  dL)  [nfWM*)  -  (2.11) 
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when  n*  (R)  and  nf  (iZ)  are  the  respective  distributions  of  pairs  which  are  radially  expanding  (+)  or  radially 
contracting  (~)  at  R,  ie.  n+  and  n"  are  respectively  characterized  by  postive  (vr  >  0)  and  negative  (vr  <0) 
radial  speeds  ur.  These  distributions  are  coupled  via  the  boundary  conditions, 

n~(R  — *  oo)  =  nf  (R) 
n;(R~)  =  nt(R;) 

at  infinity  and  the  pericenter  JZf  for  dissociated  states  i,  and  by 

nJ(RJ)  =  n+(Rj) 
nj(R+)  =  *+(J2+) 

at  the  apses  R±  of  bound  orbits  f.  Within  each  set,  the  collision  frequency  vxj  couples  the  distributions 
nf  with  all  ( Ef,L2f )  state-distributions.  The  two  sets  of integro-diflerential  coupled  equations  (2.11)  for  nf 
and  nf  are  coupled  via  their  boundary  values  (2.12)  and  (2.13).  Simple  as  (2.11)  appears,  this  general  class 
of  partial  integro-differential  equations  involving  boundary- value  conditions  are  notoriously  difficult  to  solve. 
The  Volterra  type  equation  (2.11)  has  as  yet  not  been  solved  theoretically  or  numerically.  Simpler  versions 
of  (2.11)  are  currently  receiving  much  attention  in  the  mathematical  literature12. 

Addition  of  (2.11)  for  n*  and  nf  provides  the  microscopic  continuity  equation  in  terms  of  the  microscopic 
net  radial  current, 

MR)  =  J(R,Ei,L2,)  =  [#+(*)  -  n-(J?)]u*  (2.14) 

across  R  as, 


(2.13) 


(2.12) 


_LJL 

R 2  dR 


r°°  rLt<» 

[R2«  -  n~)  |  \]B  L,  =  -  /  dEf  /  dl)  [«,(fl)M*)  ~  MRWMR)}  (2-15«) 

”  *  JV(R)  Jo 


or  as, 

1  <9  r00  rL\m 

WdR  Ir2j'^  =  ~jv[R)dEi  J0  dL * M*) -Pf(R)]Ct/(R)>  (2.156) 

which  is  expressed  in  general  terms  by, 

V-Jt  =  n<(fl)  f  dE}  f  ,m  dL)vxj{R)PhR).  (2.15c) 

Jv(R)  Jo 

Collisional  depletion  and  production  terms  involve  the  summed  distribution  nt  =  (n*  +  nf )  which 
appears  on  the  right  hand  sides  of  (2.15),  while  the  difference  (nt+  -  nf )  appears  on  the  transport  side. 
Upon  integration  of  (2.15)  between  JZf  and  oo,  the  rate  (2.1)  or  (2.8)  across  the  dissociation  limit  E  =  0 
may  then  be  expressed  in  the  equivalent  transport  representation  as, 

^  [Rlvl 

aNANB  =  iim  /  dE,  /  dl\  [4*  J^n,"  (*oH  pf  ( Ro)P,A{Ro)  (2.16) 

Ro—cc  J0  J0 

where  the  one-way  equilibrium  fiux  (Appendix  A) 


47rR2nf  vr  =  47r2Ar,4  AT^e*'£?i^T/(27rAfy4BtT)3/2 
across  R  is,  by  (2.4)  independent  of  R  and  h\.  Also  the  normalized  distributions  are, 

-(R)  =  0±(RE  ^  =  iH) 

P.  ( R )  -  P  (£.  k,  k )  -  ^ 


(2.17) 

(2.18) 
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and 


Mm- 1  "*(*>-,  <?(*> 


(2.19) 


is  introduced  here  as  the  probability  for  association  within  a  sphere  of  radius  R  of  all  contracting  i-pairs  which 
enter  the  R- sphere.  Of  the  radial  microscopic  flux  4?r R2n~  |  vr  |  at  R  only  a  fraction  P^i?)  eventually 
becomes  associated  as  exhibited  by  (2.16). 

The  set  of  Master  Equations  (2.11)  yield,  for  the  normalized  distributions  (2.18)  to 


JLJL 

R2  dR 


(jflMw)-,/" 

'  '  JV(R)  JO 


dL ) 


Vf  (*)-/£(*)]  c}(Jl) 


(2.20) 


where  the  one-way  equilibrium  current  jf  =  n*vx,  and  where  the  forward  and  backward  one-way  equilibrium 
collisional  rates  between  levels  i  B  (Et1 1%)  and  f  are, 


<%(*)  =  nf(R)u,}(R)  =  nf(R)vft(R)  =  cft(R) 
which  satisfy  detailed  balance.Then  (2.20)  with  (2.17)  reduces  to, 


dR 


J  B.,L\ 


[pf(*)-p±(*)]  (utf(R)/vR) 


This  set  is  solved  according  to  boundary  conditions. 


(2.21) 


(2.22) 


and 


P~{R~*  oo)  =  1; 


p~(R,E,  -*  oo)  =  1 


(2.23) 


p;(K)  =  pt(R7) 

pt(Rt)  =  p;(R?) 


(2.24) 


which  in  effect  serve  to  couple  both  (±)  sets  of  integro-differential  equations  (2.22).  Since  Xf  =  (2 mE,)p2 
where  p  is  the  impact  parameter,  then  the  transport  representation  (2.16)  of  the  rate  can  be  rewritten  as 
the  (R0  — *  oo)  limit  of 


a=  1  j  ee~(  dc  J  2irp  dp  p~  (R0)PtA(Ro) 

where  f  =  E,/kT,  and  where 


(2.25) 


pi  =  Rl[  1  -  V(R0)/Et)  (2.251 

is  the  maximum  impact  parameter  accessible  for  formation  of  pairs  with  R  <  Ro  at  energy  Et . 

Many  mathematical  complexities  exist  in  determining  the  theoretical  and  numerical  solution  of  the 
general  class  of  equations  subject  to  boundary  (not  initial  !)  conditions,  as  represented  by  (2.22)-(2.24). 
While  investigations  into  appropriate  numerical  procedures  are  currently  in  progress,  it  is  important  now  to 
contruct  simplified  models  suggested  by  the  respective  collisional  and  transport  forms  (2.10)  and  (2.16)  for 
the  rates  and  reported  in  the  following  sections  with  the  aim  that  new  insight  and  conclusions  will  emerge, 
without  the  necessity  of  detailed  solution.  The  models  will  then  in  turn  suggest  various  procedures  valuable 
for  eventual  detailed  numerical  solution. 


3.  Improved  Classical  Treatment 

Treatments  which  include  forward  collisional  couplings  vxj  and  neglect  backward  collisional  couplings 
Ujt  are  termed  classical  in  the  sense  that  the  set  of  coupled  equations  (2.11)  or  (2.20)  are  then  replaced  by 
the  two  sets  of  uncoupled  equations  in  nf  and  n“  which  contain  only  absorption  (loss)  terms  and  which  can 
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be  solved  in  closed  analytical  form.  Classical  solutions  in  the  sense  as  defined  here  therefore  only  include 
collisional  absorption  from  dissociated  states  (i?t,£f). 

3.1  Analytical  Rate  with  Collisional  Absorption.  Physical  Assumption: 

Neglect  the  redissociation  terms  nfV}x  which  collisionally  populate  continuum  level  i  in  (2.1)  from  the 
bound  states  f  for  R  <  Ro.  Thus  Pft(R)  the  stabilization  probability  (2.9),  is  unity  for  R  <  Ro .  Assume 
that  thermodynamic  equilibrium  in  ExyL\  (but  not  in  R)  is  collisionally  maintained  at  fixed  R  for  pairs 
with  R  >  Ro,  so  that  pt(R)  =  pf(R)  for  R  >  Ro  is  independent  of  (Rt,lf).  The  stabilization  probability 
Pfi{R  >  Ro)  °f  (2.9)  then  vanishes.  Thus  Ro  is  assumed  to  be  sufficiently  small  that  the  net  effect  of 
collisions  at  R  <  Rq  is  pure  absorption,  depleting  the  dissocation  channels,  and  yet  is  sufficiently  large  that 
the  net  (input-output)  effect  of  collisions  at  R  >  Ro  is  to  maintain  thermodynamic  equilibrium  in  Et  and 
if,  whereby  preventing  any  stabilization  from  occurring  outside  Ro- 

The  total  frequency  fcr  collisional  destruction  of  continuum  states  i  in  (2.11)  may  be  taken  for  small 
R  <  Ro  as  the  total  collisional  frequency  in  (2.10)  for  direct  formation  of  all  bound  levels,  ie. 

*,(R)=  /  dEf  i/,/(R  )dL)*l  dEf  vif(R)dL )  (3.1) 

JV(R)  Jo  Jv(R)  Jo 

In  this  approximation  the  bound  pairs  are  now  fully  stabilized  against  redissociation  for  R  <  R0  and  are 
therefore  considered  as  the  recombined  products.  The  recombination  rate  (2.9b)  then  reduces  to, 

-  -  rRoPo  rRo  v(R\ 

kNaNb  =  a(R0)n{R0)  =  /  dEt  dtf%  /  [4 *R2n,(R)  |  vR  |]  dR ,  (3.2) 

Jo  Jo  Jr;  \vR\ 

in  the  collisional  representation .  Even  though  both  representations,  collisional  (2.10)  and  transport  (2.16), 
are  equivalent,  physical  insight  and  algebraic  development  is  facilitated  more  readily  via  the  transport  form. 


3,2  Analytical  Solution  for  R  <  Ro:  Nonequilibrium  in  E%)  L\  and  R. 

On  ignoring  therefore,  the  backward  rates  n*  Vfx  for  collisional  repopulation  of  level  i,  each  set  of  coupled 
Master  equations  (2.11)  reduces  to 


=  -nf(R)MR)<  K<Ro  (3.3) 

for  pure  absorption  within  R  <  Rq.  Here  and  below,  vr  is  now  used  to  denote  the  positive  radial  speed 
]  vr  |.  Eqns.  (3.3)  are  coupled  by  boundary  conditions  (2.12)  and  (2.13).  Since  pf  =  n± /nf ,  (3.3)  is  then 


Qp?(RY 

dR 


e„l; 


=  =F 


MR) 

VR 


P\ 


R  <  Ro 


(3.4) 


which  follows  also  from  (2.22).  When  integrated  subject  to  boundary  conditions  p"(Ro)  at  Ro  and 
P~{RX  )  =  p+(R~)  at  the  pericenter  R~  (3.4)  provides  the  distributions 


p;(R)  = 


n,(R) 


p~  (R0)  exp 


dR 


for  pairs  contracting  at  R,  and 


(3.5) 


p,+  (fl)=  =  pr(#0exp 

n,  (R) 


«(R) 


dR 


vr 


>  R  <  Ro, 


(3.6) 
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for  pairs  expanding  at  R,  respectively.  It  is  now  simpler  to  insert  these  distributions  into  the  two-dimensional 
transport  representation  (2.16)  rather  than  into  the  equivalent  three-dimensional  collisional  representation 
(3.2).  In  either  case  however  the  rate  reduces  to, 


where, 


olN^Nb  = 


dL\  [4^n,-(iZo)t;fi]  P;(R0)Pa(R0) 


(3.7) 


pi  =  2m[Ei  -  V{R0)\  (3.8) 

and  where  the  probability  (2.19)  for  association  within  Ro  of  ( R  <  R0)  pairs  now  takes  the  simple  form 

PA(Ro)  =  1  -  exp  |  j  (3.9) 

for  absorption  within  the  segment  RX<R<  Ro  of  the  trajectory  enclosed  by  the  sphere  of  radius  Ro-  The 
integrand, 


dR  =  Vi(R)dt  =  ^  =  ^di,  (3.10) 

VR  At  At 

is  the  elemental  probability  that  A-B  collides  with  M  during  time  interval  dt,  or  within  trajectory  element 
dsx  of  their  actual  (EX)  It2)-orbit  s ,  for  A-B  relative  motion  at  speed  vt.:  Hence, 


p,a(Ro)  =  1  - 


pjjR  o) 
p7(Ro) 


(3.11) 


is  simply  the  probability  for  (A-B)-M  collisions  with  microscopic  free  path  length  A,  =  vxjvx  towards  any 
collision  within  the  portion  of  the  orbit  enclosed  by  a  sphere  of  radius  R0. 


Binary  Decomposition:  When  the  (A-B)-M  collisions  at  a  given  pair  separation  R(A-B)  can  be 
decomposed  as  binary  (A-M  and  B-M)  collisions  (an  excellent  decomposition  for  ionic  species  A*  -  B~  in 
a  neutral  gas),  then  ux  =  vxA  +  uxB  so  that  A~l  =  A^1  -f  A~J.  Then  (3.11)  factors  as, 

P*(Ro)  =  l- SASB  (3.12) 

which  decomposes  as, 


Pa(R0)  =  PA(R0,  Ka)  +  P,A(Ro,  Kb)  -  PA(Ro,  Ka)Pa(Ro .  Kb), 

where  the  probability  of  survival  against  either  (A-M)  or  (B-M)  collisions  within  the  segment  is 


$a,b  =  exp  <j>^  ds./A.^sj. 


(3.13) 


(3.14) 


This  decomposition  holds  only  at  this  ( Ex  L2)  microscopic  level  (cf.  §3.3).  The  corresponding  probability 
for  association  arising  from  an  individual  collision  within  the  trajectory  is 


Pxa{Ro,Ka>b)  =  1  -  $a,b  (3.15) 

The  above  bmary  assumption  leads  quite  naturally  to  the  decomposition  in  (3.13)  of  the  (EX,L2) 
microscopic  probabilities  for  individual  collisions.  Since  the  sum  Pa(Aa)  -f  PA( A#)  includes  the  probability 
P/^A^P/^Aa)  for  simultaneous  collisions  twice,  the  simultaneous  probability  must  be  subtracted  as  in 
(3.13).  Note  that  the  above  solutions  (3.5)  and  (3.6)  of  (3,3)  for  the  normalized  distributions  p±  and 
(3.11)  for  the  microscopic  probability  of  association  show  quite  directly  that  these  JJt,if-pairs  are  not  in 
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equilibrium  with  respect  to  J?x  in  this  classical  approximation  since  pf  are  clearly  dependent  on  L.  The 
distribution  p~  ( Ro )  of  pairs  contracting  at  Rq  may  be  in  i2-equilibrium. 


3.3  Averaged  Association  Rates  and  Probabilities 
With  the  aid  of  Appendix  A,  the  rate  (3.7), 

__  _  pOO  fRoPa 

a NaNb  =  /  dE  dL 2  [4 zR2n^{R0,E,L2)vR]  p'{R0,E,L2)PA{R0iE,L2)  (3.16) 

Jo  Jo 

may  be  rewritten  as, 


a  =  xv  [°°  €e‘<  de  P;{R0,Eyp)PA{Ro',E,p)dp7  (3.17) 

Jo  Jo 

where  p  is  the  impact  parameter  associated  with  the  ( Et ,  Lt)~  trajectory.  Here  L 2  =  (2 MabE)p2)  e  =  E/kT) 
v  is  the  mean  speed  (SkT/nMAs)1^  and  p\  =  -  V(Ro)/E],  On  performing  the  if  and  E-integrations, 

then  with  the  aid  of  Appendix  A, 


a  =  (NaNb)'1  j“  dE  [47 :Rln~  {R0,  E)^]p~  (Ro,  E)Pa(R0 ,  E) 

(3.18) 

=  JT vj  ee~'  de  pi  p~{Ro,E)Pa{Ro,E) 

(3.19) 

=  rp2maxvp-  (Ro)Pa(Ro)-,  pL x  =  Rll  1  ~  ^(*o)/*T] 

(3.20) 

where  the  integrated  densities  and  distributions  are  given  by 

r*lrl  ,  , 

n~(Ro,  E)=  n~  (R0,  E,  L*)dL 2  =  p~  (Ro,  E)n~  (R0,  E) 

Jo 

(3.21) 

and, 

n~(R0)=  [  n~(R0)E)dE  =  p~(Ro)n~(R0). 

Jo 

(3.22) 

Also, 

xplp-  ( Ro ,  E)Pa(Ro,  E)  =  n  P  p~  (Ro,  E,  p)PA(R0;  E,  L*)dp 2 

Jo 

(3.23) 

and 

np2maxp-(R0)PA(Ro)=  fee-'  de  np20p- (R0,E)PA(R0,E), 

Jo 

(3.24) 

define  the  appropriate  if  and  (if,£t)  averages,  PA(Ro,E)  and  PA(Ro),  respectively,  of  the  original 
microscopic  association  probability  (2.19). 

In  the  absorption  limit, 


Pa(Ro,  E,  L2)  =  1  —  exp 


(3.25) 


and  for  general  interactions  V(R),  analytical  expressions  for  Pa(Rq,  E)  and  Pa(Rq)  can  be  deduced  under 
various  equilibrium  conditions. 

For  R0  sufficiently  large,  the  distribution  of  pairs  contracting  at  Ro  is  in  equilibrium  with  respect  to  L. 
Then 
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»r(R  ,e,l*)  =  -+ 


n"(R ,E,L2)  _ 


S"(R,  E ) 


n~(R,  E) 


(3.26) 


so  that  p  (R0y  E ,  L2)  =  p  (Ro,E)  is  now  independent  of  L.  The  X2-averaged  probability  in  (3.23)  becomes 

If  in  addition  the  incoming  JZo-pairs  are  also  in  equilibrium  with  respect  to  E  so  that  p~(R0>E)  = 
p~(i2o)*  The  E-averaged  or  macroscopic  probability  is  then, 


pA(Ro)=-^~  denplPA(Ro,E)  =  l~^\.  (3.28) 

" Pmax  Jo  n 

Although  the  microscopic  absorption  probability  (3.25)  satisfies  the  rule  (3.13)  for  decomposition 
into  probabilities  for  individual  A  and  B  collisions  with  the  gas,  the  trajectory  p  and  energy  E-averaged 
probabilities  (3.27)  and  (3.28)  do  not  satisfy  this  decomposition.  That  is, 


Pa(R0,  A;  E)  *  1  -  (1  -  P,)(l  -  P3)  -  Pi  +  P2  -  Pi P3  (3.29) 

where  Pi  =  Pa(R0,\a',E)  and  P2  =  Pa(R0)Xb',E).  The  equality  has  always  been  tacitly  assumed 
in  previous  macroscopic  treatments13”15.  The  breakdown  is  directly  apparent  from  comparison  of  the 
corresponding  expr  ssions  derived  explicitly  for  straight-line  trajectories  (E  >  V(/Z)  for  weak  interactions). 
Here,  the  Thomson  straight-line  (Eyp)  microscopic  probability18  is 

Psl(Ro\  A;  E  -  oo,  p)  =  1  -  exp{-2(i?o  -  PJ)1/J/A}  (3.30) 

with  p-average14, 


psl(X>  E  00 )  =  psl(X)  =  1  —  2^2  [1  ~  e~2X^  +  2^))i  X  =  Ro/X  (3.31) 

which  does  not  rigorously  satisfy  the  decomposition  rule  (3.13).  The  decomposition  (3.13)  however  is  satisified 
numerically  fairly  closely.  While  the  Thomson  microscopic  probability  (3.30)  depends  on  L  (via  p),  the  I2- 
average  (3.31)  is  independent  of  E  since  E  >  V(iZ).  The  Thomson  distributions  within  i?o  are  therefore 
in  E-equilibrium  but  not  in  X2-  equilibrium.  The  distributions  which  yield  (3.25)  are  in  general  not  in 
(E}  X2)-equilibrium. 

For  ion-neutral  (hard-sphere)  collisions  with  constant  (speed  independent)  path  length  A,  -  A  then 
(3.25)  reduces  to, 


Pa{Ro ,  E,  L2)  =  1  -  exp  [-£(JZo,  E,  I2)/A]  (3.32) 

where  C  is  the  length  of  the  segment  of  the  (E,L2)- trajectory  enclosed  by  the  sphere.  For  (ion-induced 
dipole)  collisions  at  constant  frequency  ux  =  v  then 

Pa(Ro,  Et  L2)  =  1  —  exp(— i/T)  (3.33) 

where  T  is  the  transit  time.  Figure  l(a,b)  illustrate  the  variation  with  A'  =  Ro/X  of  the  p2  and  (p2,i?)- 
averaged  values  (3.27)  and  (3.28)  of  the  microscopic  probability  (3.32)  associated  with  ion-ion  recombination 
under  Coulombic  attraction  V(R)  —  ~e2 /R  =  —  ( RtjR)kT  within  a  characteristic  radius  i?o  =  0.408i?e  (cf. 
§4.4  and  ref.  4).  Various  analytical  expressions  for  the  microscopic  segment  length  C  of  the  (£\  L)-trajectory 
have  been  derived  in  Appendix  B  (cf.  (B6),  (B9),  (B15)  and  (B19)),  together  with  corresponding  L2  or  p2 
averages  (B31)  for  general  interactions  V(iZ). 

As  E  increases,  the  p2-averaged  probability  (3.27)  decreases  monotonically  (see  Fig.  1)  from  the  parabolic 
(E  <  V(R))  envelope, 
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(3.34) 


PA{X  =  Ro/X,  E  -  0)  =  2  [y/y/Tv]  V  exp  (-2X(1  -  y2)1/2)  y  dy, 

obtained  from  (B15b)  for  £  in  (3.32)  with  (B22),  to  the  straight  line  (E  >  V(iZ))  envelope  (3.31).  This 
monotonic  decrease  between  the  two  limiting  envelopes  essentially  arises  from  the  decrease  of  the  p2-averaged 
segment  C(E>p]Ro)  of  the  enclosed  trajectory  from  2Ro  (cf.  B23)  for  parabolic  motion  (cf.  B23)  to  |fZo 
for  rectilinear  motion.  The  envelopes  increase  initially  with  density  N  ^  A"1  as  the  collision  probabilities 
2 R0/X  and  |iZ0/A  in  the  parabolic  and  straight  L:ne  limits,  respectively.  For  intermediate  energies  E  the 
gradients  are  confined  to  within  these  limits  (cf.  Fig.  la,  inset).  Fig.  (lb)  illustrates  that  the  E-averaged 
probability  (3.28)  is  graphically  indistinquishable  (to  within  three  significant  figures)  from  PA(Ro,E  =  kT) 
and  is  much  closer  to  the  parabolic  limit  (3.34),  than  to  the  straight-line  limit  (3.31). 

The  use  of  (B19)  for  £  in  (3.27)  is  universal  in  that  the  variation  of  PA{Ro,  E)  with  Ro  at  a  fixed  gas 
density  N,  in  addition  to  the  above  variation  with  N  at  a  fixed  Ro,  can  be  illustrated  also  by  Fig.  1  simply 
by  regarding  the  selected  values  of  the  normalized  energy  a s  values  of  e'  =  roe,  and  the  absissa  X  as  linear 
in  Ro.  As  R0  =  r0Rt  increases  from  zero,  PA  follows  for  all  e  the  initial  increase  of  the  e'  =  0  parabolic 
envelope  and  eventually  falls  onto  the  straight-line  envelope.  The  iZo-variation  oi  P  is  shown  in  Fig.  (2a)  for 
E  -  kT  at  one  density  (/  =  Re/X  =  1)  and  in  Fig.  (2b)  at  several  densities  ranging  from  high  to  low. 

There  are  as  yet  no  analytical  results  for  the  averaged  probabilities  (3.27)  for  Pa(Rq,E)  and  (3.28) 
for  PA(Ro)  for  general  V(R),  not  even  for  the  Coulomb  interaction.  Figs.  (la,b)  are  based  on  numerical 
/^-integrations  over  Elliptic  Integrals  (Appendix  B).  New  analytical  expressions  now  will  be  developed  for 
general  V(R)  for  the  following  two  cases:  (a)  The  low-gas  density  limit  (when  PA  ->  £/At)  in  §4.1  and  (b) 
The  case  of  pairs  in  L2 -equilibrium  (§6.1). 


4.  Low  Density  Limit:  Classical  Analytical  Rates  and  Probabilities  Under  Exact  (A-B) 
Trajectories 

The  association  probability  (3.11)  tends,  as  N  -♦  0  ie.  as  A t/Rt  — ♦  oo,  where  R<  is  a  characteristic 
length  given  by  the  outermost  root  of  |  V(Re)  1=  kT ,  to 

PA(R0,Et,L2)^  (4.1) 

which  is  linear  in  N.  The  transport  form  (3.16)  of  the  classical  rate  then  tends  to, 

c*lou.NaNb  =  f  dEt  dL?[4*Rln;{Ro,EuL2t)vR]fR  ^  (4.2) 

which  provides  the  required  linear  increase  with  N,  Since  the  collision  frequency, 

(43) 

At  At 

is  already  linear  in  N  via  An  the  low-density  limit  of  the  colhsional  form  is  obtained  from  (3.2)  by  simply 
setting  nt(JZ)  =  n,(iZ),  the  zero-order  approximation  for  all  dissociated  pairs,  to  give  the  low  density  N  linear 
limit, 


poo  pH  pft\  i 

aL(m,NANB  =  dEt  /  rfif  /  [47rfl2n,(#,Et,I?)](r,/A,)di? 

Jo  Jo  Jr~ 

which  is  identical  with  (4.2)  since  dt  -  dR/vR  =  dst/v,.  Since  the  one-way  equilibrium  flux, 


4it2NaNb  exp (~E,/kT) 


iltR  "•  E"  L' ^  VR  1  (2itMABkT)V2 

across  a  sphere  of  radius  R  is  independent  of  R  and  Lx  and  since 


(4.4) 


(4.5) 


11 


L7  —  (2  mEt)p2 

in  terms  of  the  impact  parameter  p,  then  (4.2)  is  also 


(4.6) 


qLou!  — 


'  8fcr 
*Mab 


(e~e  de 


ds , 

X 


(4.7) 


where  e  =  Et/kT  is  the  normalized  energy.  The  maximum  impact  parameter  po  accessible  for  formation  of 
the  i?o-complexes  at  energy  Ex  is  given  by, 


Po  —  Rlpl/(2mEt)  —  - (4-8) 

Even  when  A,  is  independent  of  R  and  p,  (4.7)  involves  knowledge  of  the  length  C(Eyp)  of  the  enclosed 
trajectory  as  a  function  of  E  and  p.  For  pure  Coulombic  attraction  £(E,p)  can  be  provided  (Appendix  B)  in 
terms  of  incomplete  Elliptic  Integrals  of  the  first  and  second  kind,  but  the  double  (£,  p)-integration  remaining 
in  (4.7)  appears  intractable.  A  more  elegant  approach  is  based  on  the  recognition  that  the  innermost  double 
integral  of  (4.7)  is  simply, 


(4.9) 


This  identity  is  proven  directly  in  Appendix  B  (B32-B36),  but  follows  most  conveniently  by  reversing  the 
order  of  the  ( L7x ,  ^-integrations  in  the  collisional  form  (4.4)  of  the  rate.  Upon  assuming  A,  is  independent 
only  of  I,,  then  (4.4)  gives 


ociowNaNb  =  f  dEt  f  4*R7vx(dR/\t)  f 
Jo  Jo  Jo 

The  Lf-integrated  equilibrium  distribution  here  is  then  the  standard  Maxwell-Boltzmann  distribution, 
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nx(R,  Ex,  L7)dL7 . 


(4.10) 


nt(R,  Et)dEx  = 


E'  1  NaNb  exp(~Et/kT)d(E,/kT) 


v/tt  [  kT 


(4.11) 


such  that 


nt(R,  Ex)vtdEt 


'  8 kT  ’ 

1/2 

*MAb. 

V(R) 


Et 


NaNb€c  €  de, 


with  the  result  that  (4.10)  then  reduces  to 


(4.12) 


OIL  o 


8 kT  11/2  r°°  .  ,  , 

— — —  /  ee  de  4rr  i?2 

f,  v(*)] 

■kMab\  Jo  Jo 

Et  J 

dR 


J  A,(J?,£,) 


(4.13) 


Comparison  with  the  transport  form  (4.7)  provides  the  valuable  relation  (4.9)  valid  for  all  (curved) 
trajectories  under  general  V(R)  when  the  radius  Rq  is  independent  of  Lx,  but  not  necessarily  of  Et.  When 
A,  and  R o  are  both  independent  of  e  then  integration  of  (4.13)  over  e  yields, 


&  Lou>  — 


SkT  11/2  fRo 


kMab  J 


j; 


4t tR2 


1  - 


V(R) 


kT 


dR/X(R) 


(4.14) 


Both  (4.13)  and  (4.14)  are  exact  (new)  classical  results  for  all  trajectories  under  general  V(R)  valid  when 
the  path  length  A,  is  or  is  not  dependent  of  J?,,  respectively.  The  customary  classical  (Thomson  straight-line) 
result  follows  from  the  direct  use  of  rectilinear  paths  in  (4.7)  or  by  setting  Ex  >  V(R)  in  (4.13). 
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4.1  Averaged  Association  Probabilties 
The  rate  (3.18)  may  also  be  written  as, 


aLowNANB  =  J™  PA(E,-,  R0)dE,  (4.15) 

where  PA(Et,  i?o),  the  microscopic  probability  (4.1)  averaged  over  p  (or  If),  is  the  probability  for  association 
of  an  equilibrium  number  of  impairs  contracting  across  By  comparison  with  (4.12)  and  (4.13),  the 
probability  at  low  gas  density  is  given  exactly  by, 

_  i  vi mi  <ir  _ 


V(*)l 

dR 

E,  J 

WE,) 

where  p\  is  given  by  (4.8).  Since  the  ^-averaged  flux 


J  ^nt(R0, Et)vi(R0)4irRl  dEt  =  np2mazv 


and  where  t>  is  the  mean  relative  speed  (SkT/xMAs)1^)  then  (4,15)  reduces  simply  to, 

gLou’NaNb  ~  Trp^ axFPA(ito),  (4.19) 

where  PA  is  the  (Et}  p)-averaged  of  the  microscopic  probability  (3.28)  given  exactly  by, 

Pa(Ro)  =  -1-  4jri?2 fl  -  W\dR/X(R)  (4.20) 

*Pma*  JO  L  kI 

when  Rq  and  At  are  assumed  independent  of  Et.  The  probability  (4.16)  and  its  £, -average  (4.20)  hold  for  any 
curved  trajectory  under  V(R).  Expressions  (4.16)  and  (4.20)  are  new  analytical  results  for  the  probabilities 
under  general  V(R). 

4.2  Ion-Ion  Recombination  Under  Coulombic  Attraction: 

For  hyperbolic  motion  under  Coulombic  attraction  VC{R )  -  -e2/#>  the  £, -microscopic  probability 
(4.16)  reduces  to, 

where  both  R o  and  A,  may  be  functions  of  E,.  The  macroscopic  probability  (4.20)  is, 


4  Ro  f,  3  VC(R0) 


PA(R o)  =  ~  i- 


]/H 


v  u'  3  A  1  2  kT  J/  [  kT  J 

for  constant  Ro  and  A.  The  rate  (4.14)  or  (4.19)  is  exactly, 

&Low(Rq)  =  *p2maxvPA(Ro)  =  J  1  - 

for  Coulomb  attraction.  This  is  appropriate  for  termolecular  ion-ion  recombination, 

A+  +  B~  +  M  -  AB  +  M 

in  a  low  density  gas  M.  For  straight-line  trajectories,  Thomson  obtained  the  rate, 
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(4.25) 


QJJ*  —  irR0vPT  (J£q) 

where  the  Thomson  (straight-line)  probability  is 


_  4 

T  “  3T 


(4.26) 


The  corresponding  probability  (4.22)  and  cross  section  (^p^occ)  for  the  correct  hyperbolic  trajectory 
are  greater  than  the  corresponding  straight-line  Thomson  probability  P?  and  cross  section  itR^  by  the 
enhancement  factors, 


and 


Fv 


(4.27  a) 


^c  =  Pma*/^0  = 


(4.276) 


respectively,  where  Re  =  e7/kT  is  the  natural  unit  of  length.  The  probability  factor  Fr  remains  bounded, 
increasing  from  1  for  large  R0  »  Rt  to  1.5  in  the  limit  of  small  R0  <  The  focusing  factor  Fc  is  unlimited 
increasing  from  1  to  Re/Ro  over  the  same  range  of  Ro*  The  two  factors  however  combine  to  give  the  amplified 
ratio, 


w«r= (i +ff;)-{!: 


W'/Ro) 


Ro  ^  Rt 
Ro  <  Rt 


(4.28) 


For  typical  Ro  %  Rt/2  characteristic4  of  ion-ion  recombination  this  trajectory-correction  for  a  is  quite  large 

M). 

Consideration  of  the  inner  integral  of  (4.7)  in  isolation  presumedly  hindered  previous  efforts  to  obtain 
the  correct  generalisation  of  Thomson’s  straight-line  result  (4.26)  to  actual  hyperbolic  trajectories.  It  can  be 
seen  however  that  the  cross  section  focusing  factor  Fc  to  the  cross  section  is  the  main  correction.  Natanson’s 
result15, 


<*N  =  *Pmox  x’Pt{Ro)  (4.29) 

which  accounts  for  Fc  but  retains  the  Thomson  straight-line  association  probability  (4.26)  therefore  remains 
quite  accurate  by  being  lower  than  the  correct  rate  (4.23)  by  the  Fr-correction  of  between  1  and  1.5  . 


4.3  Electron-Ion  Recombination 

For  termolecular  recombination  electron-ion  recombination, 


e~  +A+  +  +  (4.30) 

the  kinetic  energy  of  the  electron  of  mass  m  changes  by  an  amount  (T,  -  T/)  =  (2m/M)Tt  after  an  elastic 
(low-energy,  isotropic)  collision  with  the  gas  atom  B  of  mass  M.  The  change  in  internal  energy  Ex  of  the 
(e“  -  A+)  pair  is  2?,  -  Ef  =  T,  -  T/  =  (2 m/M)[Et  -  F(jf?)j.  Hence  the  pair  is  bound  (Ef  <  0)  provided16, 


Et< 


2m  e2 

A 7r 


(4.31) 


where  Em  rather  than  oo,  is  now’  the  upper  limit  to  E%  in  (4.13). 

On  rearranging  the  order  of  integrations  to  reflect  this  R-dependent  limit  to  EX)  the  rate  (4.13)  for 
constant  A,  is  then, 


<*L  = 


8  kT 


kMab 


nl/2,  rR»  r6RJR  / 

}  A  l  l  ( 


Re 


ee 


de 


(4.32) 
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(4.33) 


which,  with  the  aid  of  e  (6i2«/*)  ^  1  —  6RC/R  for  small  <5  =  (2 m/M),  yields, 


q5,  =  8?r 


m  RoR\ 

m~t 


SkT 


ll/2 


xMab  J 


where  Re  =  e2  /kT.,  When  Ro  =  §i2e  the  Thomson  radius,  this  classical  rate  aeL  agrees  exactly  with 
Pitaevsku’s  rate17  derived  from  a  Fokker-  Planck  analysis  of  diffusion  in  energy  space.  A  previous  classical 
treatment10  of  (4.30)  established  the  important  result  that  ae  varies  linearly  with  Ro  (as  above)  rather 
than  Rq  as  in  (4.23)  for  ion-ion  recombination  (4.24).  The  formulation16  adopted  however  a  less  rigorous 
weighting  procedure  which  differed  from  (4.13)  by  a  factor  of  2/3  (in  the  inner  integral)  and  which  therefore 
resulted  in  a  rate  § ap  rather  than  ap  as  found  here. 

The  main  difficulty  in  applying  (4.23)  and  (4.33)  to  the  collisional  recombination  processes  (4.24)  and 
(4.30)  is  that  the  trapping  radius  Ro  is  uncertain.  The  rates  increase  without  limit  as  R$  and  R0  for  each  case, 
respectively.  The  radius  may  be  assigned  de-facto  by  comparison  with  results4  of  more  elaborate  theories 
based  on  accurate  numerical  solution  of  the  full  Collisional  Input-Output  Master  Equation.  It  would  however 
be  advantageous  if  a  classical  treatment  were  formulated  where  R0  appears  as  a  variational  parameter.  Such 
a  treatment  is  now  developed  in  the  following  section. 


4.4  Classical  Variational  Treatment 

Let  the  bound  level  V'(iZo)  separate  the  reaction  zone  from  the  product  zone.  The  one-way  rate  across 
this  transition  state  from  (2.1)  with  E  =  V(i20)  is, 

_  _  r°°  rRlvl  rRt 

av(Ro)NANB  =  /  dE,  dL ?  /  n,(R)ul(R)dR  (4.34) 

Jv(R0)  Jo  Jr ~ 

where  the  frequency  of  collisional  transitions  across  the  boundary  is 

rv(R„)  .i‘,H 

v't{R,Ro)=  /  dE,  u„{R)dL  (4.35) 

Jv{R)  JO 

and  where  R+  is  the  apocenter  for  bounded  motion  (0  >  Ex  >  V'(Jff))  and  is  Ro  for  dissociated  pairs. 
The  contribution  ac  to  (4.34)  from  the  continuum  states  is  given  by  (4.13)  with  vx  replaced  by  v[  and  is 
supplemented  by  the  additional  contribution  ap  from  the  bound  states  between  0  and  V(Ro)  to  give, 


ay  =  Olc  +  Ck£ 


where, 


t°  rR°  dR 

aB  ~  /  dEl  I  [4*i?2nt(R,r>,]  v>/p  W\ 
Jv(Ro)  Jo 


'V(Ro) 

which  with  the  aid  of  (4.12)  for  bound  states  reduces  to 

p-Ro 


-W-jOo  *1  ’“'"’(Tff1'1)? 


(4.36) 


(4.37) 


(4.38) 


where  p(c)  =  n,/n,.  «o  =  ~V(R0)/kT  and  A|(f0)  =  is  in  principle  a  function  of  co.  At  low  gas  densities 
p  is  independent  of  R  due  to  R-equilibrium. 

For  Coulomb  attraction  and  constant  A(,  then  the  contribution  to  (4.34)  from  the  continuum  states  is, 


<±c(Ro)  =  (*  +  2€°)  >(-Re/Ro 

as  before  in  (4.23).  When  p( f)  =  1  and  the  contribution  from  bound  levels  is, 


(4.39) 
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(4.40) 


<*b{Ro)  =  g  *£<>•£( 
The  summed  rate  (4.36)  is  therefore 


1+2€0 


av{Ro)  =  \*Roj; 


1  - 


Vc(i?o) 


2  kT 


,  3 

l  +  2«o 


exp(-y(iio)/*T) 


as  a  function  of  #0i  or  is 


(4.41) 


«v(e0)  = 


(4.42) 


when  expressed  as  a  function  of  the  binding  energy  €o  =  e2 /{RokT)  of  the  transition  level.  Since  upward 
transitions  past  this  level  V“(i?o)  are  neglected  by  (4.41)  and  since  nx  is  set  equal  to  nx  which  is  an  upper 
limit  for  the  distribution  of  bound  levels  in  the  range  0  >  Et  >  V(J?0)>  the  rate  (4.41)  is  therefore  an  upper 
limit.  On  assuming  that  the  variation  of  A'(eo)  in  (4.42)  with  €o  is  much  slower  than  the  remaining  eo- 
dependence,  (4.42)  displays  a  minimum  at  Cq  =  </6  which  corresponds  to  R £  =  0.408RC.  This  level  to  then 
acts  as  a  bottleneck  to  the  recombination  as  in  the  bottleneck  method  of  Byron  et  al.  This  value  compares 
very  favorably  with  exact  assignment  of  Ro  (cf.  Fig.  2  of  reference  4). 

The  bound  levels  between  V(Rq)  2.45&T  and  the  dissociated  limit  are  in  general  not  in  energy 
equilibrium  as  assumed  in  (4.40)  and  also  in  the  original  bottleneck  method1®.  On  setting  a  trial  (variational) 
non-equilibrium  distribution, 


p(«)  = 


exp(-f/e0) 


which  is  physically  realistic4 ,  in  (4.38)  then  (4.34)  yields, 


(4.43) 


4  V 
av  =  -*Rtjj 


11  (7c0  -  3eg)  -  (2eg  -  c0  -  5)exp(t0  -  1) 

A  2eo  2e0(e0  -  l)3 


(4.44) 


which  displays  a  minimum  at  ej  =  2.55  which  corresponds  to  R*0  =  0.392J2e.  The  choice  of  R*0  is  therefore 
rather  insensitive  to  p(e).  The  ratios  of  (4.42)  and  (4.44)  to  |tt.R*(F/A')  are  however  1.75  and  1.56, 
respectively. 

Since  A(  =  vx/u[  is  in  general  different  from  the  mean  path  A,  for  production  of  bound  from  dissociated 
pairs,  and  is  in  general  unknown,  the  present  strategy  is  to  assign  the  unknown  radius  Ro  in  §4.2  and  §4.3 
as  the  above  variationally  determined  R*0 .  Hence  (4.23)  yields, 


aL{Rm0)  =  0.3178^ffJ?J  (4.45) 

This  result  is  only  within  a  factor  of  2  higher  than  previous  numerical  results4  of  the  exact  treatment 
in  which  the  quasi-steady  state  solutions  (n/  of  a  collisional  input-output  integral  equation  for  bound-pairs 
Ef  is  used  in  equation  (2.11)  of  §2)  over  a  wide  range  (0.1  <  a  <  1)  of  mass  parameters, 


a  =  Mg/(2MX  +  Mg)  (4.46) 

for  recombination  of  ions  of  equal  mass  Mx  in  a  gas  of  mass  Mg . 

It  is  also  interesting  to  note  that  agreement  of  (4.45)  with  the  Thomson  straight-line  rate-equation  (4.25) 
with  Ro  =  |  Re  but  with  v  set  equal  to  the  customary  RMS  value  of  1.0854  v , 

aT  =  0.3216-Ti^Y  (4-47) 

3  A 

must  be  regarded  as  fortuitous  in  that  the  Coulomb  focusing  factor  of  4.676  (neglected  in  (4.45))  offsets 
the  considerable  reduction  (Rq/Rt)z  ~  0.23  of  the  reaction  volume.  It  is  also  interesting  to  note  that  the 
variational  value  R0  =  0.408Re  agrees  closely  with  the  value  R n  =  0.417Rf  in  Natanson’s  rate, 
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<xN  =  it. R2nv  1  -  Pt(rn)  (4.48) 

which  differs  from  (4.19)  only  in  the  use  in  P?  of  Thomson’s  straight-  line  probability  (4.26)  rather  than 
the  exact  probability  (4.22)  associated  with  hyperbolic  orbits. 

5.  Angular-Momentum  Integration:  if -Equilibrium,  (R.  f?)- Nonequilibrium 

5.1  Master  Equations 

Integration  of  the  basic  set  (2,11)  of  Master  Equations  for  nf(R)El,L^)  over  the  full  range  0  <tf< 
R2p f  =  2 mR2[Et  -  F(R)j  =  Xj^R)  of  if  accessible  for  given  R  and  Et  is  facilitated  by  Leibnitz’s  rule 
which  provides, 

["  F55  I*’"- < * 1 « {f  ls’"‘ <*>  I  **  'I 

r  ^9  r  2  ■ 

-nf(R,&,.&<)[^^  (5.1) 

for  differentiation  with  respect  to  R  of  an  integral  with  R-dependent  limits.  The  distribution  n*(R,  Et)  X^J 
in  (5.1)  is  evaluated  for  that  maximum  angular  momentum  i^(R)  =  R2pf (R)  associated  with  the  ( Et ,  X^t)- 
orbit  which  just  reaches  R  at  either  turning  point  R,f.  Since  contracting  pairs  n~  are  transformed  into 
expanding  pairs  nt+  at  the  pericenter  and  expanding  n*  pairs  are  transformed  into  contracting  pairs  at  the 
apocenter,  the  R-space  must  be  divided  into  two  Regions  I  and  II  which  exclusively  contain  the  percenters 
R~  and  the  apocenters  Rt+,  respectively  for  a  given  Ey  and  all  accessible  X,.  In  Region  I  sources  of 
nf(R~)  originate  from  n“(R~),  while  in  Region  II  sources  of  n~(Rt+)  originate  from  nt+(J?t+).  The  effective 
interaction  (2,2)  associated  with  L2mt  is, 

Mfi)  =  V(R)  +  ^  (M 

Region  I  which  contains  the  pericenters  of  all  Xf-orbits  with  a  given  Et  has  range  R~0  <  R  <  Rb  and  Region 
II  which  contains  the  apocenters  is  Rb  <  R  <  R+0  where  the  boundary  radius  Rb  is  the  radius  of  the  circular 
orbit  given  by  the  minimum  of  Vm%}  ie.  by  the  zero  of 

dVm,  =  dV  2(E,-V)_  1  dL*nt 

dR  dR  R  2 M  -  ABR 2  dR  [  '  ‘ 

and  Rf0  are  the  turning  points  associated  with  the  most  penetrating  X,  =  0  orbit.  Thus  Rb  =  e7/2  \  E%  \ 
for  a  bound  orbit  under  Coulombic  attraction.  Region  II  does  not  exist  for  dissociated  pairs.  The  boundary 
conditions  (2.12)  and  (2.13)  can  therefore  be  incorporated  into  (5.1)  by  setting  the  distribution  evaluated  at 
the  turning  points  to  be, 

n*/D  p  j2  X  _  fn;(R,E„R2pJ)\  Region  I  (dVm,/dR  <  0)  .  . 

n,  ( R ,  Lmt)  -  |  n+(J2)  ^  jjyj  Region  n  {dVmt/dR  >  o)  W 

since  n“  is  the  precursor  (source)  of  n*  in  the  pericenter  Region  I  and  n*  is  the  precursor  of  n~  in  the 
apocenter  Region  II  when  it  exists,  integration  of  (2,10)  with  the  aid  of  (5.1)  and  (5.4)  therefore  yields, 

=  * [”■'  i  **  I  ]  ”  [wiE‘  t'  dL>  t iL-  ['?<*>  -  <4W]  <*R>  <5-5> 

where  C,f(R)  —  n,utj(R)  is  the  one-way  equilibrium  collisional  rate  and  where, 
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R2p2 

j±(R)  =  j±{R,  Ei )  =  f  '  n,(R,  Eu  if)  \  vR  \  dl\ 
Jo 


(5.6) 


is  the  current  per  unit  dKdEt  across  the  R-sphere. 

Addition  of  the  two  equations  in  (5.5)  eliminates  n*  and  shows  that  the  2?t -microscopic  net  current 


(5.7) 


Jt(R)  =  J(RtEt)=[j+{R)-j~(R)] 
across  R  satisfies  the  microscopic  continuity  equation, 

1  d  ,  ^  ,  r°°  tL™s  «  fL'n' 

WdR  =  -  Jm)  dEi  l  dLf  l  M*)  -  P/(*)]  C„(R)  (5.8) 

in  general.  The  net  effect  on  the  net  current  (5.7)  of  the  conversions  at  the  turning  points  are  null.  For 
(JE'j,  £f)-equilibrium,  pt(R)  =  Pf{R)  =  p(R)  and  (5.8)  then  implies  constant  flux  4tt R2Jt(R)  across  the 
R-sphere. 

^-Equilibrium:  For  equilibrium  in  i2,  then 


n^R,^,!,2) 


Ht(R,  Eu  if)1 


n,(R,  2?t) 


(5.9) 


nt(R, Et)  J 

which  is  (Appendix  A), 

n;(R,  Eti  Li )  =  ^  2R2p  pr  )  n‘ (^>  )  (5.10) 

where  pt  and  pr  are  the  relative  and  radial  momenta,  and  where  the  if-integrated  distribution  is 

n<(R,  Et)  =  f  n,(R,  E„  I?)dl?  (5.11) 

Jo 

The  normalized  distribution  pt(R,Ex,L f)  of  (2.7)  is  then  independent  of  if.  The  current  (5.6)  is  then, 

;?(*.£.)  =  J#*(R ,E,)vt  =  (I^)nf(R,£„I^)  I  VR  I .  (5.12) 

Under  If-equilibrium  (5.5)  therefore  reduces  to  the  set, 


1  8 


±RidR  t'R2j‘  ~±L2 


1 

raiLi 

mx 

ijL 

.  dR  . 

E, 


j;(R,Et)  +  J  [pj(i2)  -  pf(R)]  C?f{R)dEh  R<Rh  (5.13) 


for  R  in  Region  1  (where  dI?mJdR  >  0),  and  to, 


1  Q  1  QT  2  poo 

±FdR  lR2£(R'E'^  =  =F  jr  ~ef  E  tf(R’E>)  +  J  [pfW-ptW }  C*f(R)dEj,  R>Rt  (5.14) 


for  R  in  Region  II  (where  dL2mJdR  <  0).  The  one-way  equilibrium  rate  (per  unit  dE,dEj )  for  Ex  — ♦  Ej 
collisional  transitions  is, 

fR2T:  fR-v) 

C±(R)=  dlf  /  n^lf(E>,L2t;E},Ly,R)dL)  =  nt(R<E>)utJ(El,Ef,R)  (5.15) 
Jo  Jo 

Note  that  d£^t(i?)  at  constant  Ex  in  (5.13)  is  simply  the  increase  (if  positive)  or  decrease  (if  negative)  in  the 
number  of  orbits  (states)  L^t  with  their  turning  points  within  the  range  between  R  and  R  -f  dR.  In  Region 
I  which  contains  the  pericenters,  this  number  increases  as  R  is  increased  and  in  Region  II  which  contains 
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the  apocenters  this  number  decreases  as  R  increases.  Thus  dL^/L^t  (5.13)  and  (5.14)  is  the  fractional 
increase  (or  decrease)  in  the  current  ^(R,  Et)  across  the  R-sphere  as  R  is  increased  to  R  +  dR.  Equation 
(5.13)  therefore  incorporates  directly  the  boundary  conditions  (2.12)  and  (2.13)  by  including  the  appropriate 
source/sink  terms.  Contracting  (— )  pairs  which  disappear  as  (— )  pairs  reappear  as  sources  of  expanding 
(-f)  pairs  in  the  pericenter  region.  Equation  (5.14)  acknowledges  that  (+)  pairs  are  sources  of  (— )  pairs  in 
apocenter  Region  II.  For  Coulombic  attraction,  Region  I  comprises  all  R  for  Et  >  0  and  0  <  R  <  Ri  for 
Et  <  0  while  Region  II  comprises  the  range  Ri  <  R  <  2 Rb  for  Et  <  0.  The  boundary  radius  Rt,  =  e2/2  \  Et  j 
is  the  radius  of  the  circular  orbit  of  energy  Et. 

The  net  radial  microscopic  current, 

MR)  =  m .  £.)  =  -  n;)v ,  =  \nf{p+  -  p;)v,  =  [Pt+(R)  -  P;{R)\  jf  (5.16) 

satisfies  the  equation 


~  [R'MR))^  =  l^{p,(R)-p,(R))Ctf(R)dE}, 


(5.17) 


since  pt  =  (p+  +  pt~)/ 2  and  Ct/  =  2 Cfj(R).  The  net  effect  of  the  geometrical  (+)  — ►  (-)  conversions  in 
(5.13)  and  (5.14)  to  the  net  current  Jt  is  obviously  null,  in  accord  with  (5.17). 

The  rate  (2.1)  can  then  be  written  in  terms  of  the  association  probability 


PA(r  RV  l  Pt(R) 

in  the  transport  representation  as, 

aNANB=]im  [  dE,\i*R2n;(R)^]p;  {R)P,A{R) 

R—  oo  Jq  l  2-1 

where  i  labels  £t-state  quantities.  In  terms  of  the  probability, 

v  n.t'.y  p,(R) 


(5.18) 


(5.19) 


(5.20) 


for  subsequent  collisional  stabilization  of  Ej -pairs  originally  produced  collisionally  from  a  dissociated  J?,-pair, 
the  collisional  representation  of  the  rate  is, 


c*NaNb  =  lim  [  dE,  f  dR  f  nt{R)vtf{R)Pft(R)dE , 

R-°°J0  Jo  Jv(R) 

exp  (-Et/kT) 


In  terms  of  the  one-way  equilibrium  flux 

4*RJjf {K,E,)  =  4 xR2  (^(R.S.Jy)  = 


(2rr  MABkT)V2 


per  unit  interval  dEt  across  the  R-sphere,  then  the  source-sink  terms  of  (5.13)  and  (5.14)  are, 


1 

\dLL] 

.  dR  . 

(5.21) 


(5.22) 


(5.23) 


The  set  (5.13)  and  (5.14)  then  reduce  to  the  following  equations  in  Region  I  and  Region  II,  respectively. 

±*4  (*’?)  =  /“  (5.24) 


(**?<?<*))  *  =  [m  ["? 1 <*>  -  «*(*>]  CWVE>-  <5-25> 
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When  (— )  and  (+)  pairs  are  simultaneously  in  i£t-equilibrium  then  pf(R)  =  pf{R)  =  p*(.R),  and  the 
collisional  RHS  of  (5.24)  and  (5.25)  vanishes.  The  LHS  then  predicts  constant  flux  47 rR2jf[p+(R)  -  p“(P)] 
for  the  exterior  R  >  Ro  region.  The  above  set  (5.24)  and  (5,25)  for  I2-equilibrium  is  identical  with  that 
deduced  previously19  from  conservation  considerations. 

Classical  solutions  are  now  developed  in  §6  under  I2-equilibrium  in  the  interior  (R  <  Ro)  region. 
Dissociated  pairs  belong  only  to  pericenter  Region  II  and  classical  (full  absorption)  solutions  to  (5.24)  will 
be  provided  in  §6.  Bound  pairs  with  energy  E  <Eh-  V(R)  -f  (1/2 )RdV/dR  and  R  >  Rbi  the  zero  of  (5.3), 
are  in  Region  II  governed  by  (5.25),  which  is  also  solved  in  §6.2  in  the  absorption  limit. 


0.  I2-Equilibrium:  Analytical  Expressions  for  Classical  Association  Rates  and 
Probabilities  for  a  General  Interaction  V'(P) 

Assume  as  in  §3.1  that  Pj?t,  the  probability  (5.20)  for  stabilization  of  the  bound  state  f  collisionally 
produced  from  dissociated  state  i  is  unity  for  R  <  R0  and  is  zero  for  R  >  P0*  The  rate  (5.21)  in  the 
collisional  representation  is  then, 


where, 


aNANB  ~ 


nt(R)ut(R)dR 


Vx 


ulf(R)dEj 


is  the  frequency  for  collisional  production  of  bound  levels. 


(6.1) 


(6.2) 


6.1  Pericenter  Region  I  {E  >  Eb  -  V{R)  +  \%\ R  <  Rb  for  E  >  0  and  all  R):  Microscopic 
Distribution  of  Dissociated  and  Bound  Pairs. 

The  absorption  form  of  (5.24)  is  then, 


where  for  small  R, 


2 uJjR) 
vt(R) 


p7(R)  = 


(6.3) 


v% 


vij{R)dEj  ^  f  vx}{R)dEf 

JV{R) 


(6.4) 


is  close  to  the  total  collision  frequency  to  all  levels.  On  integrating  (6.3)  between  R  and  R0 ,  the  solution  is, 


p7(r)  -  p7  (ro)  exP 


p;{Ro)Ps(RoxR), 


(6.5) 


which  in  effect  defines  PJ  (Po>  P),  the  probability  of  survival  against  collision  of  (-)  pairs  from  Ro  to  P.  By 
comparison  with  (3.5),  the  main  consequence  of  if -equilibrium  on  the  distribution  is  to  simply  replace  the 
radial  speed  vr  by  its  I2-average  of  ~vx. 

The  outgoing  distribution  satisfies  (5.24)  which  reduces  in  the  collisional  absorption  limit  to, 


8  [W]t^(jiW)4!*y]»:w 


or 


(6.6a) 


which  is, 


i?yP,v«/A- 


=  ^(Ry)p:{R)e 


d 

Jr 
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(6.6b) 


The  above  equations  (6.3)  and  (6.6)  also  follow  directly  from  Li- integration  of  the  corresponding  (Et>  In¬ 
equations  (3.3)  over  all  If,  use  of  Leibnitz’s  integration  rule  and  setting  pf(R>E<R2p2)  =  p~(R,  Et)  so  as 
to  acknowledge  conversion  at  pericenter  R~.,  The  solution  of  (6.6b)  can  be  written  as, 

P,+(S)=  [5  p;(R7)u>-(R;)p+(R;,s)dR;  (6.7) 

Jo 

for  S  in  the  range  0  <  5  <  i?o  of  Region  I.  The  physics  of  this  result  is  that  the  source  of  outgoing  pairs 
are  the  incoming  pairs  which  undergo  conversion  at  each  pericenter  R~  between  R  =  0  and  S.  The  outgoing 
pairs  formed  at  R”  then  survive  to  S  with  probability 


psiR7's)  =  exP 


(6.8) 


The  fractional  weight  of  incoming  pairs  with  pericenters  R”  within  the  interval  dR~  about  R“  is, 


w;{r;)<ir: 


SVs 


(6.9) 


under  X2-equilibrium.  Pairs  with  anguiar  momentum  Lmax  —  Sp,  have  a  relative  orbit  which  just  touches 
the  S-sphere  at  S  which  in  turn  encloses  all  the  pericenters  R~  between  0  and  S. 

The  probability  of  association, 


PtA(Et,R0)  =  l- 


riVUt  q) 
p;(e%,Ro)' 


within  Ro ,  with  (6.5)  and  (6.7),  yields 


(6.10) 


PtA(Et,  Ro)  =  l-  ~  r  Ps(Ro,  R; K  (K)Ps(K  >  Ro)dR; .  (6.11) 

-^oPo  Jo 

The  physics  of  this  result  is  that  the  incoming  (-)  pairs  which  arrive  at  the  pericenter  with  probability 
RJ  (R0,  Rt  )  and  weight  (R~ )  per  unit  dR~  are  converted  at  each  R~  to  outgoing  (+)  pairs  which  survive 
to  Ro  with  a  probability  P^(R^}R0)  With  the  above  definitions  of  P$  and  integration  by  parts  then 
yields, 


P,A(E„  Ro)  =  -3  fR°  4nR2  (l- 
*Po  Jo  \ 


T) 


exp 


-Tt(Ro-R) 


dR/X„ 


(6.12) 


where  the  maximum  impact  parameter  po  which  has  its  pericenter  on  Ro  is  given  by  p2  =  Rq[1  -  l7(Ro)/It] 
and  where  A,  may  depend  on  Et .  In  the  limit  of  low  gas  densities  Ro/A,  — ►  0  and  (6.12)  reproduces  the 
previous  low-density  limit  (4.16). 

For  energy  independent  A,,  the  macroscopic  probability  (3.28)  with  (6.12)  yields, 


^(i?o) 


1 

^Pmcix 


ym 

kT 


'dR' 

A 


(6.13) 


Expressions  (6.12)  and  (6.13)  pertain  to  general  interactions  V(R)  and  are  valid  when  the  (Et,  If)  pairs  with 
R  <  Ro  are  in  if  equilibrium  and  in  (E%,  lf)-equilibrium,  respectively.  The  corresponding  distributions 
(6.5)  and  (6.7)  at  Ro  associated  with  PA  are  not  in  R-equilibrium. 

On  recalling  the  original  complexity  entailed  with  the  multidimensional  integrations  (3.23)  and  (3.24) 
involving  the  microscopic  distributions  (3.5)  and  (3.6)  and  the  probabilities  (3.25),  both  expressions  (6.12) 
and  (6.13)  represent  a  significant  algebraic  reduction  applicable  to  general  V(R), 


6.2  Apocenter  Region  II  {E  <  £*  =  F(R)*f  |^;R  ^  R*>):  Microscopic  Distribution  of  Bound 
Pairs. 
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Equation  (5.25)  reduces  in  the  absorption  limit  to, 


for  pf  and  to 


(6.14) 


—  [ R2P  p,  (R)  exp(-2J2/A)]  =  - 


pf(R)exp(-2R/X) 


for  pt  .  The  solution  of  (6.15)  subject  to  pre-assigned  px  (i?o)  is, 


(6.15) 


P."(5)  =  -0sp;  (Ro)Ps(Ro,  s)  +  Js  p?{RW(Rt)P;(Rt,S)dR?  (6.16) 

for  S  in  the  range  Rb  <  S  <  Ro.  The  first  term  on  the  RHS  is  the  direct  contribution  from  incoming  iZcr pairs 
which  survive  against  collision  to  S  with  probability, 


PJ  {Ro,S)  =  exp 


(6.17) 


and  RIpI/SVs  is  the  fractional  number  of  instates  accessible  under  i2-equilibrium.  The  second  term  is  the 
contribution  from  apocenter  (~f)  — »  (-)  conversions  at  i?t+  which  survive  against  collision  back  to  5  <  Rt 
with  probability  Ps  .  The  fractional  weight  of  apocenter  pairs  within  interval  dR+  is, 


»t(Rt)dRt  = 


SJp|(S)  • 


The  solution  of  (6.14)  subject  to  p,+  (Rb)  at  the  boundary  between  Regions  I  and  II  is, 


(6.18) 


p?{R?)  =  p+(Rt>)P+{Rb,Rt),  Regionll  (6.19) 

where  (6.8)  provides  the  probability  for  survival  of  outgoing  pairs  to  R*  in  the  range  Rb  <  R*  <  Ro-  The 
distribution  p^  at  Rb  is  obtained  from  the  pericenter  Region  I  solution  (6.7)  which  gives, 

pt(Ri)  =  p7 (Rb)  fR'  ps(Rb> (R;)P^(R7,Rb)dR;t  (6.20) 

Jo 

which  is  by  now  self  explanatory.  On  using  (6.16)  for  p^{Rb)  then  in  terms  of  the  assigned  pt~(i?o)j 


Pi+  {R' } = { §ipr  {Ro)ps  {R°' Ri) + C p' {R'  )tt,*+  {R')ps (R'  - Rb)dR' } 


X  |  J*h  p- (Rb<  r~ X  (R;  )p+  (r;  ,  Rb)dR; 


Pf(Rb,  i?t+),  Regionll  (6.21) 


The  notation  adopted  here  facilitates  the  recognition  of  the  various  sequences  involved  in  the  generation 
of  p*  from  p~  in  Region  II,  It  also  illustrates  quite  dramatically  that  pt+  is  the  solution  of  an  integral 
equation  in  R-space.  The  use  of  the  absorption  limit  adopted  in  this  section  §6  essentially  eliminates  the 
integral  equation  in  E-space  involved  with  the  original  Master  set  of  equations  (5.24)  and  (5.25).  The  full 
solution  of  the  Master  equations  involves  integral  equations  in  E  and  R  space,  which  are  difficult  to  solve 
numerically. 

The  distribution  p~(R)  in  Region  I  which  originates  from  p~{Ro)  in  Region  II  is, 

Pi  {R)  =  P,~  {Rb)Ps  ( Rb ,  R)\  Regionl  (6.22) 
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6.3  Application:  Coulomb  Association  Probabilities  for  Pericenter  Region  I. 

Dissociated  pairs  with  all  R  are  all  in  Region  I.  For  Coulombic  attraction  V(R)  =  -e2/JZ,  and  constant 
free-pat  h-length  A,  (6.12)  yields, 


P*(Ro,  E)  =  i  -  jij  [(1  -  e"4* )  (  jr  [x2  -  X’]  -  l)  +  AX 

(6.23) 

where  the  parameters 

X  =  R0/X-,  X(e)  =  Po(e)/Ro 

(6.24) 

are  related  by 

X2(R0,e,X)=  l  +  j±  X2 

(6.25) 

in  terms  of  the  natural  units  Re  =  e7/kT  and  e  =  E/kT  of  length  and  energy. 
The  low-energy  parabolic  ( E  — ►  0)  limit  to  (6.12)  yields 


1).  1-A&-.H*) -.{»/*■  XxZl  (««> 

and  the  high-energy  rectilinear  (E  — >  oo)  limit  is, 

PA(Ro ,  E  -  oo)  =  1  -  -L  [iX  -  1  -  e"4X]  -  (  fR°/X'  (6-27) 

The  initial  gradients  of  PA  are  proportional  to  the  enclosed  parabolic  and  chord  segments  2Ro  and  |i?oi 
respectively.  The  above  limits  (6.26)  and  (6.27)  are  the  envelopes  to  the  family  of  plots  of  (6.23)  versus  gas 
denisty  for  various  energies  e.  The  E-integrated  probability  (6.13)  is, 

Pa(R0)  =  PA{Ro,kT)  (6.28) 

and  is  identical  with  the  probability  (6.23)  evaluated  at  e  =  1.  The  above  probabilities  derived  under  the 
assumption  of  X2- equilibrium  will  not  in  general  agree  with  the  corresponding  results  (3.27),  (3.31)  and 
(3.34)  based  on  non- equilibrium  distributions  (3.5)  and  (3.6),  except  in  the  low  density  (A”  — ♦  0)  limit  when 
the  gradients  of  all  the  envelopes  are  equal  to  2X  and  |A”,  respectively.  Figure  4  provides  the  X-variation 
of  the  E-averaged  envelopes  (6.26)  and  (6,27),  and  probabilities  (6.23)  for  various  e.  The  probabilities  are  in 
general  lower  than  the  corresponding  probabilities  of  fig.  (la)  associated  with  i2-nonequilibrium  since  not 
all  collisions  result  in  absorption  (cf.  fig.  5). 

All  collisions  with  (Et,L*)~ pairs  are  not  expected  to  result  in  stabilised  recombination,  eg.  those 
involving  only  angular-momentum  redistributions  in  ( EUL f  — »  f„lj)  transitions  can  be  excluded.  These 
collisions  which  involve  only  angular  momentum  changes  are  expected  to  promote  equilibrium  in  X2  rather 
than  absorption.  In  this  sense  therefore  the  classical  (absorption)  solutions  (6.12)  and  (6.13)  associated  with 
X2-equilibrium  represent  a  more  natural  outcome  than  those  (3.5)  and  (3.6)  associated  with  full  absorption 
of  all  if  levels  upon  collision.  The  probabilities  (6.12)  or  (6.23)  and  (6.13)  or  (6.28)  are  therefore  expected 
to  provide  a  representation  more  physically  correct  than  (3.27)  and  (3,31)  ie.  the  lower  curve  of  fig.  5  is  to 
be  preferred. 

7.  Solution  as  a  Function  of  Gas  Density:  Reaction  and  Transport  Rates 
The  rate  (2.14)  is  exactly, 


qlNaNb  —  lim 

—  oc 


J*00  /»J?yp0 

/  dE,  I  dLl[4*Rln;(Ro)vR}p;{Ro)P,A(Ro) 
Jo  Jo 


(7.1) 
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where  the  microscopic  distributions  pf  implicit  in  PA  (=  1  —  pf  / p~ )  are  solutions  of  the  transport-collisional 
Master  Equations  (2,20)  for  general  non-equilibrium  in  (Lf,i?i,i?)  coordinates,  or  of  (5.24)  and  (5.25)  for 
non-equilibrium  only  in  Ex  and  R.  This  procedure  provides  in  principle  the  frll  variation  of  a  with  gas  density 
N. 

The  pairs  with  large  R  are  generally  in  (Ei,L*)  equilibrium  but  not  in  R  equilibrium.  Further  insight 
and  development  may  however  be  achieved  by  address  to  the  macroscopic  continuity  equation.  The  net 
microscopic  current  J(R,  i?,  if)  and  the  net  If -integrated  current  J(R,EX)  satisfy  microscopic  equations 
of  continuity  (2.15b)  and  (5.8),  respectively.  On  integrating  (5.8)  which  holds  quite  generally  for  non- 
equilibrium  in  R,  Et  and  If,  over  all  Et  accessible  at  a  given  R  and  use  of  Leibnitz’s  rule, 


/“ ,  hh  = vis  {*’£/<*  w } + 


-  0 


is  zero  due  to  the  null  effect  of  all  the  collisional  redistribution  terms, 


r  dEt  r  dE,[pf(R)-pi(R)}Ctf{R)  =  0. 

JV(R)  JV(R) 


(7.2) 


(7.3) 


Since  the  net  current  J(Rl)Et  =  V{R))  =  5(11*  -  n~)t>t  vanishes  at  the  outermost  turning  point  Rt 


given  by  Ex  =  V’(Ji),  (7.2)  provides  the  macroscopic  continuity  equation, 

VJ(R)  =  0 

which  is  always  valid  for  the  net  macroscopic  radial  current 


J(R)=  r  J,(R)dE,=  r  dEi\v.+{R,El)-n-(R,Et)}^, 
JV(R)  Jv(R)  1 

rL  m, 

=  dE ,  dl?  [n+  (R,  E„  if)  -  n~  (R,  E„  if)]  vr, 

JV{R)  JO 

/•OO 

-/  dE,  /  dlf  [p+(i?)-prw]  (nfvR). 

JV(R)  Jo 


(7.4) 

(7.5) 

(7.6) 

(7.7) 


Assumption:  For  large  R  >  Rq ,  the  pairs  are  in  (Et)  lf)-equilibrium  ie.  p±(R)  =  P±(R)  irrespective  of  i. 
Hence  (7.7)  yields, 

Jl(R)  =  j+(R)-]-(R)=  l-NANB'xp(-V(R)/kT)v\p+(R)-p-(R)}t  R>  R0  (7.8) 

when  the  bound  levels  in  the  range  [0  -  V(R)}  are  included  in  the  integration,  or 

F(p+(J?)-p-(fl)]  (7.9) 

when  they  are  excluded.  For  V  <  kT>  then  J i  — ►  J2.  Under  steady-state  conditions 

4t tR2J(R  >  Ro)/NaNb  =  MR  -  00)  =  nRl  (l  -  1^1)  vPa(R0)p~ {Ro)  (7.10) 

where  p(R)  is  the  non-equilibrium  normalized  distribution  n(R)/h(R)  of  expanding  and  contracting  pairs 
and  is  related  to  p~{R)  by 


MR)  =  7  NaNb 

4 


kT 


p{R)=\m±tm= 


l--PA(R)\p-(R) 


(7.11) 
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when  RA(Ro)  is  the  macroscopic  probability  for  association  within  R.  Since  the  incoming  pairs  are  generated 
at  infinity  with  a  Boltzmann  distribution  p~(R  — >  oo)  — >  1  and  p(R  — ►  oo)  -+  1  -  }>PA(<x>).  The  summed 
distribution  p(oo)  is  therefore  not  in  E-equilibrium.  Also  p“(oo)  >  p+(oo)  and  a  net  inward  current  (7.8) 
exists  at  infinity. 


7.1  Non- Equilibrium  p(R):  The  Transport  Rate 

The  basic  Master  Equations  (2.20)  followed  from  (1.2)  which  was  originally  derived  from  a  Boltzmann 
equation  in  (R,  p)  phase  space.  In  the  absence  of  a  sink  at  R,  it  was  also  shown6  from  the  same  Boltzmann 
equation  that  (7.5)  the  macroscopic  current  J  can  be  represented  by  the  standard  diffusion-drift  expression, 


3(R)  =  -DVn(R)  - 


— VV(R) 
€ 


n(R),  R  >  R0 


(7.12) 


where  n(R)  is  the  non-equilibrium  (A-B)  pair  distribution  which  is  affected  by  the  strength  of  sinks  within 
R  <  Ro.  The  macroscopic  collisional  coefficients  D  and  K  for  relative  diffusion  and  drift  of  A  towards  B 
through  a  gas  M  are  given  respectively  by  D  =  Da  +  Db  and  K  =  Ka  4*  Kb  where  Da,b  and  Ka,b  are  the 
diffusion  coefficients  and  mobilities  for  the  individual  species  A  and  B  in  M,  These  coefficients  at  low  field 
strengths  (as  Coulombic  at  large  R)  are  related  by  KkT  =  De  or  by  Ke  =  DRt  where  Rt  is  the  natural 
unit  e2  /kT  so  that 


m 

NaNb 


~D'* 


The  solution  of  (7.13)  between  Ro  and  oo  yields, 


(7.13) 


where  the  transport  rate, 


p{R0)  =  p( oo)  (1  -  a/aTR{Ro)) 


(7.14) 


aTR{Ro)  —  4xD 


/  txp{V{R)/kT)R-2dR 
Jro 


(7.15) 


is  the  rate  of  production  by  diffusional-drift  of  E0-pairs  from  pairs  originally  with  infinite  separation. 
Once  p(R0)  can  be  furnished  from  (7.1)  in  terms  of  a,  then  (7.14)  can  be  solved  for  a. 


7.2  The  Reaction  Rate 

We  now  choose  Ro  in  (7.1)  to  be  finite,  yet  sufficiently  large  that  no  sinks  exists  outside  Ro  (eg. 
Ro  >  b(e2/kT)  for  a  Coulomb  interaction4  at  low  gas  densities).  Both  (7.4)  and  (7.13)  are  then  valid.  The 
effect  of  non-reactive  collisions  outside  Ro  will  now  be  acknowledged  via  appropriate  choice  of  p~  (Ro)  in 
(7.1).  The  component  species  A  and  B  are  assumed  to  be  in  collisional  equilibrium  with  their  mutual  field 
V{R)  at  these  large  R  >  Ro.  The  kinetic  energy  gained  from  V(R)  by  the  contracting  pairs  as  R  is  decreased 
is  lost  upon  collision  with  the  gas  species  M.  Following  each  collision,  the  kinetic  energy  T  at  each  R  >  Ro 
is  therefore  reduced  to  the  original  kinetic  energy  Ex  =  Too  of  generation  of  pairs  with  infinite  R.  Assuming 
that  the  last  collision  suffered  by  the  Repair  occurs  at  (R,  +  A),  where  A  is  the  macroscopic  mean  free  path, 
then  the  pre-collision  kinetic  energy  of  the  R-pairs  is, 

T'(R)  =  Too  +  [V(R  +  A)  -  V{R) }  (7.16) 

which  represents  an  increase  of  <5F(R)  =  [^(R  +  A)  -  ^(R)]  over  its  original  thermal  value  Too  resulting 
from  the  previous  collision.  Prior  to  the  next  collision  at  R  the  internal  energy  is  then, 

E'X(R)  =  T'(R)  4-  V(R)  =  Too  -r  V{R  4-  A).  (7.17) 

At  low  gas  densities  A  — +  oc,  and  Et  of  course  remains  conserved  for  all  R  >  Ro  while  T,  increases 
from  Toe  =  Et  at  infinity  to  [Et  -  Vr(R)]  at  R.  The  internal  energy  of  pairs  contracting  at  Ro  is 
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E[  =  Too-  |  V{Ro  +  A)  |  rather  than  Ei  =  T*,  as  it  would  be  in  the  absence  of  diifusional-drift  collisions 
outside  Ro. 

The  degradation  of  Et  from  its  value  at  infinity  therefore  causes  Et  in  (7.1)  to  be  replaced  by  Hence 
(7.1)  yields, 


<*  =  <*rn{Ro)  [p{Ro)/p{oo)] 


=  (naNb)  1  r  dE,  r°V°dL\  \4irI%Z(Bo)vR]  PA(R0 ) 

V  '  JV[R0+X)  JO 


(7.18) 

(7.19) 


This  defines  the  rate  of  reaction  oirn(Ro)  within  Ro  which  equals  a  if  full  equilibrium,  p(Ro)  =  1,  were 
to  exist.  On  inserting  (7.18)  into  (7.13)  then  the  overall  rate  is, 

a  =  (7.20) 

<*rn{Ro)  +  <*tr{Ro) 

This  decomposition  (7.20)  of  a  into  reaction  and  transport  rates  was  known  in  recombination  since 
the  macroscopic  result  of  Natanson16  was  recast20  in  that  form  (which  also  can  be  deduced  directly  from 
Debye-Smoluchowski  theory21).  What  is  new  here  is  the  following  expression  to  be  used  in  (7.20)  for  the 
reaction  rate, 


<*Rn(Ro)  =  ^  [a'RN(R0)p-(Ro)}  (7.21) 

where  the  macroscopic  equivalent  (3.20)  of  the  integral  in  (7.19)  is, 

<4jv(*o)  =  *p*maxT>txp(-V(R0  +  X)/kT)PA(R0)  (7.22) 

where, 

r?  -  P2 

Pmaz  —  ^0 

With  the  use  of  (7.11)  in  (7.21)  then, 

/1-Woo)\ 

<*Rn{Ro)  =  TpA^j  J  aRNiRo)  (7-24) 

since  p~(oo)  — *  1.  Previous  treatments  had  tacitly  assumed  that  olrx  in  (7.20)  was  given  by  a'RN  the 
microscopic  rate  (7.22),  The  factor  in  (7,24)  which  essentially  arises  from  the  differences  between  the  values 
p~(J?o)/p’(oo)  and  p(R0)/p(o o),  tends  to  unity  in  the  limits  of  low  and  high  gas  densities  when  PA  — +  0 
and  1,  respectively. 


i  V{R0)  |  V(R0  +  \) 


kT 


kT 


(7.23) 


8.  Summary 

Basic  sets  of  Coupled  Collisional-Transport  Master  Equations  have  therefore  been  developed  in  this 
paper  for  the  microscopic  distributions  ni(R,  £\  L2)  of  expanding  (+)  and  contracting  (— )  pairs  (A-B)  in 
a  gas  M  of  variable  density  with  respect  to  their  (A-B)  internal  separation  R,  internal  (orbital)  energy  E 
and  orbital  angular  momentum  L,  as  well  as  for  the  corresponding  leveraged  distributions  n±(R,  E), 
Expressions  have  been  introduced  for  the  microscopic  probability  PA(R)  for  association  of  (E%>  if)-pairs 
with  separations  <  R  and  for  the  probability  P^(R)  °f  multicollisional  stabilization  of  (Ej ,  i2)-pairs  with 
separation  R  formed  collisionally  from  (J?,,  Xf)-pairs.  The  recombination  rate  has  then  been  expressed  in 
the  collisional  form  and  the  transport  form  which  is  more  amenable  to  physical  insight  and  to  algebraic 
development. 

Each  set  has  been  solved  exactly  in  closed  analytical  form  in  the  (classical)  limit  of  collisional 
absorption  where  backward  collisional  couplings  Vft  are  neglected.  The  corresponding  probabilities  PtA{R) 
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for  association  have  also  been  determined.  In  particular  the  multidimensional  integrations  inherent  in  PA  for 
orbits  under  general  V(iZ),  including  Coulombic  attraction,  have  been  reduced  to  new  and  simple  analytical 
results  for  association  (a)  at  low  gas  density,  and  (b)  at  all  gas  densities  for  pairs  which  are  maintained  in 
X2-equilibiium  by  collision. 

A  classical  Variational  treatment  for  which  the  highly  excited  bound  pairs  need  not  be  in  E-equilibrium 
has  also  been  introduced,  in  order  to  provide  assignment  to  the  reaction  radius  Ro  adopted  in  the  above 
classical  treatments. 

Expression  for  the  rate  a  as  a  function  of  gas  density  in  terms  of  the  reaction  and  transport  rates 
<xrh  and  arn,  respectively  ,  has  been  proposed.  It  differs  from  previous  well  known  formula  in  the  form 
of  the  macroscopic  reaction  rate  ax#.  Finally,  useful  expressions  for  the  microscopic  path  length  in  E1 X2- 
Coulombic  Orbit  enclosed  by  a  sphere  of  radius  R  and  simple  analytical  expressions  for  the  X2-averaged 
enclosed  length  over  all  trajectories  under  a  general  V(i2)  have  been  derived  in  Appendix  B. 
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Appendix  A;  One-Way  Equilibrium  Microscopic  Flux  and  Distributions 

These  distributions  can  be  calculated  directly21  from  the  relevant  statistical  formula28  or  from  the  simple 
differentiation  of  action  integrals24.  It  is  convenient  here  to  summarise  the  key  results.  M  is  now  the  reduced 
mass  of  the  (A-B)  pair,  and  with  internal  separation,  energy  and  orbital  angular  momentum  R,  E  and  L, 
respectively.  E  =  ekT  and  X2  =  2 MEp2\  where  p  is  the  impact  parameter  for  dissociated  A-B  pairs  with 
energy  E,  with  radial  speed  =  dR/dt  and  mean  speed  v  =  (SkT/irMAs)1^- 

1.  (Ef  L)~Micro$copic  Current  and  Flux: 

7±(i?,£\I2)  =  n±(R,i?1I>* 

47rJ?2j±(i?)  E,L*)dEdL7  =  exp(-J S/kT)dEdL2 

=  Ft d(  c?(xp2) 

2.  E-Microscopic  Current  and  Flux: 

f  (R,  E)  =  I '  i*  (R,  E,  I2)dl2  =  -tm*  (R,  E) 

) 

=  np2(R)vf  e~(  de 


4 irR7j*  (R,  E)dE  =  ve  e~f  denR‘ 


V(R) 


3,  Macroscopic  Current  and  Flux: 
The  current  is 


~j±(R)=  [  r(R,E)dE=-ve-V(~R),kT 
Jv(R)  4 

for  all  bound  and  dissociated  pairs,  and  is 


~j*(R)=  r~j±(R,E)dE=  ir 


v(Ry 

kT 


for  dissociated  pairs.  Respective  fluxes  are, 
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1- 


v(R) 


kT 


4xJ?2J±  (i?)  =  7riZ2Fexp(-y(i?)/fcT) 

4  (*)  =  7rPmax^;  Pmat  =  ** 

(£?,  L )-Microscopic  Distribution: 

4jrR2n±(R,  E,  L^dRdEdL*  =  HeM~E/kT)dRdEdL 2 

5,,  E-Microscopic  Distribution: 

itR7n±{R,E)dRdE='^-(e-t  dt(l- ~^)dR 
=  47i\R2di2 


£-  y 


JfcT 


* 


=  (4^2dE)  [G&b(«,*)*] 


where  (7^  is  the  one-way  Maxwell- Boltzmann  velocity  distribution  in  above  [  ]. 
6,  Macroscopic  Distribution: 


n±(/?)  =  /  n±(R,  £)d£  =  J  exp  {-V{R)/kT) 

JV(R)  2 

for  available  bound  and  dissociated  pairs  and, 

nJ(R)  =  rn±(R,E)dE 

Jo 

=  iexp(-y(i?)/*r)  {l  -  -  -^X1/2exp(-y/*T)]  } 

for  dissociated  pairs,  where  the  error  or  probability  integral  is, 

in  terms  of  A"  =  V(R)/kT . 


Appendix  B:  Microscopic  Path  Length  £(Etp\Ro)  of  (2?,p)-Coulombic  Orbit  Enclosed  by 
Sphere  of  Radius  i?o  and  p-averaged  Path  Length  for  General  V(J2)  . 

The  enclosed  length  may  be  determined  either  from 


C  = 


which  implies  knowledge  of  the  relevant  orbit  0(J?),  or  from 


m 


rR  °  t>,  rR,> 

C(E ,  L)  R0)  =  2  — dR  =  2 

Jr,  vR  Jr, 


y/E  -  V(R)dR 


V-E- V(fl)-L2/2M/?2’ 


(B  2) 


where  the  pericenter  R,  is  the  root  of  the  denominator.  For  a  hyperbolic  orbit  under  Coulomb  attraction 
V(R)  =  -e7/R,  (B.2)  reduces  to, 


tRu  y/{R  T  a)2  -  a2  _  /“*“+0 

a(t-l)  ^(fl  +  a)2  —  a2£2  /at  Vx2  -  a2c2 


(B3) 
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where  the  eccentricity  e  and  semi-major  axis  a  are  given  by, 

■’-•♦issr1*"*  “=e!/2£- 

The  maximum  impact  parameter  for  the  enclosed  family  of  trajectories  is  given  by, 


m = «s(i  -  '^p1) 


Integral  (B3),  with  a  singularity  at  the  lower  limit,  may  be  evaluated  numerically  by  standard  integration 
techniques  or  from  the  analytical  expression, 

C(E,  L]  R0)  =  -  - —  P(4>,  A)  -  aeE(<f>,  k)  +  {R0  +  a)  sin(<£)  (B6) 

ae 

in  terms  of  the  incomplete  elliptic  integrals  F  and  E  of  the  first  and  second  kinds  respectively,  with  angle  <j> 
and  modulus  k  given  by, 

(*o+«)>- «v_(Bo+ «)’ -(*:+«)’ 


7  •  1  ° 

k  =  sm  q  =  -  =  — - „  -,t 

£  (p7  +  a-y/2 

respectively.  The  deflection  angle  V’  is  (ir  -  2a).  On  eliminating  c  via  (B4)  then  (B6)  is, 


C(E,  p;  R0)  =  2 (p2  +  a2)1/3[A,2F(«>,  A)  -  E(<j>,  k)}  +  2 (pg  +  a2)1/2  sin 


.2  _j_  . 


where, 


A'2  as  1  -  A2  =  p2/(p2  +  o2). 
Since  the  pericenter  i?,  =  a(e  -  1)  satisfies 


(Rx  4-  a)2  =  p2  4-  a2  (#11) 

from  (B4),  then  (B9)  is  alternatively  recast  in  terms  of  energy  E  and  Rx  regarded  as  new  independent 
variables,  as 


C(E>  RX]R0)  =  2 (Rt  +  a)[k,2F(<f>x  k)  -  E{<f>}  k)]  +  2(i?0  +  a) sin  <t> 


where 


( J?,  +  a)2 

For  the  limiting  case  of  parabolic  motion  e  — *  l(fc  — ♦  l,fc#  — >  0),  E  — >  0  (a  — »  oo)  direct  use  of  the  orbit, 

R{0)  =r  2^/(1  4-  cos0),  Rt{L)  =  L2 /me7  (514) 

in  (Bl)  or  on  setting  E  =  0  in  (B2)  yields, 


pRo 

£(E  -  0,  Rt\  R0)  =  2  \R/(R  -  R,)  11/2 dfl, 


=  2iZ,  In  tan  -  (?r  4-  #o)  +  2fJosin  -0o, 
4  2 


(515a) 

(5156) 
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where  0o  =  0(£o).  This  is  also  the  corresponding  limit  of  (B12)  where  E(<j),  k)  must  be  expanded  to  first 
order  in  fc'2,  where  F(<f> ,  1)  =  In  tan  -  (x  +  0o)  and  where, 

sin7<f>  — ►  sin2  =  1  -  Rx/Rq,  (£16) 

such  that  0o  =  2<£o*  For  computational  purposes  the  probability  (3.11)  with  (B9)  for  C  is  recast  in  the 
dimensionless  units, 

P  =  p/po;  «'  =£f]^  =  6ro;  X-Ro/X]  Re  =  e2/kT  ( B17 ) 

as, 

PA(X ,  e')  =  f  [1  -  exp (-£(*,  p')/\))  dp'2-,  {BIS) 

Jo 

with, 


A 


=  2X 


«♦* 


1/2 


[  VV2  +  «'2 


F{4>,  k)  -  (, 0 12  +  al2)1/2E{4>, k)  +  (1  +  a'2)1/2  sin  <t> 


(319) 


where, 


a'2  =  (4eV  +  l)]"1 ;  sin2  =  1  -  p'2;  *2  =  a'2(p'2  +  a'2)"1.  (320) 

The  variation  of  this  probability  (B18)  with  the  density  parameter  X  gas  density)  and  reduced  energy 
parameter  e*  or  ro  is  displayed  in  Fig.  1,  Expression  (B19)  is  universal  in  that  it  simultaneously  provides 
variation  with  gas  density  X  and  Rq  (see  Fig.  2). 


Averaged  Enclosed  Length  C(E\  Ro)  Over  All  Accessible  Trajectories. 

At  low  gas  densities  the  p-averaged  collision  probability  (3.27)  for  constant  A  is  {C)p/ A  where  the  average 
enclosed  length  of  all  possible  trajectory  segments, 


C(E\  Ro)  =  (C{E,  p;  Ro))p  =  4  C(E' «  (321) 

Po  Jo 

enters  into  expression  (4.7)  for  the  low-density  rate  of  ion-ion  recombination.  Direct  use  of  (B9)  in  (B21) 
appears  cumbersome,  except  for  the  limiting  case  (£  — ♦  0,  a  — ►  oo,  e  — ►  1)  of  parabolic  motion.  Then  use  of 
the  parabolic  segment  (B15b)  in  (B21)  with, 


dp  _  2 (Rj  -f-  a)dRt  dR% 

Pp)  ’  (Ro  +  a)2  -a2 


(B22) 


in  this  parabolic  limit.  The  p2-averaged  value  (B21)  after  elementary  integration  of  each  term  of  (Bl5b) 
yields 


C(E  — >  0;  Ro)  =  (^  +  ^3o  =  230 
which  also  follows  from  (Bl5a)  by  noting  that, 


(323a) 


2  [R  o  [Ro  I  "d  o  [Ro  [R  i 

c(e - o, Ro) -  yo  dRl  ta^2S° 


(3236) 
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At  high  energies  E  >  V(5),  po  — >  Ro  and  the  p2-average  (B21)  of  £  =  (Ro  —  Po)1^2  ^  |Ro*  The  enclosed 
length  (B21)  (and  hence  (3.27))  will  therefore  decrease  from  2Ro  to  |5o  as  E  increases  for  all  gas  densities 
(cf.  gradient  slopes  of  inset  figure  la). 

For  the  case  of  electron-ion  recombination  the  infinite  limit  to  the  energy  E  in  (4.7)  is  replaced16  by 
Em  =  (2m/M)(e7 /R)  so  that  the  low  density  rate  (4.7)  is, 


^  pRo  pR-O  p€ni 

Oil  =  t-Ro  /  dRt  <j)  dst  /  e  e~e  de(l  +  Re/eRo)^ 
*  Jo  Jrx  Jo 


where  here  e  =  E/kT .  This  reduces  to  first  order  in  ( m/M )  to, 

*7  9™  r^l2  rR o 


v2 m  I*  e2  1  [Ro  r 
aL~  XM[kT\  Jo  * JRi 


which,  for  the  low-energy  parabolic  motion  involves, 


fRo  ds  1 

<p  —  =  4  In  tan  -  (tt  +  0O) » 

Jrx  R  4 


the  average  of  ^  over  each  trajectory  rather  than  the  path  length  (B15).  Hence  (B25)  reduces  to 


ml)  (  e7  V 

=  8irMX\kr) 


for  collisional  electron-ion  recombination,  in  exact  agreement  with  the  energy-diffusion  result  of  Pitaevskii17 
when  Ro  =  | Rt  the  Thomson  radius. 

For  the  general  case  of  ion-ion  recombination  the  averaged  path  length  (B21)  is  required  but  direct  use 
of  (B9)  appears  prohibitive.  On  using  (Bl  1)  however,  for  the  transformation  p—»  i?“  at  constant  E,  then 

1  fPo  PR  0  a  pRv  pRa 

C(E\  R0)  =  —  /  dp7  &  d$t  =  -j  /  (Rx  +  a )dRt  /  f(Rt  Rx)dR  (528) 

Po  Jo  Jrx  Po  Jo  Jrx 


where  for  a  Coulomb  interaction  14(5), 


/(*.*.)  = 


(5  +  a)7  -  a 2 
(R  -fa)2  -  (5,  +  a)2 


is  the  integrand  of  (B3)  expressed  in  terms  of  R%  rather  than  t ,  On  interchanging  the  order  of  integrations, 
then  (B28)  reduces  to, 


a  pRo  pR 

C{E ;  50)  =  -T  /  dR  (Rx  +  ct)/(5,  Rt)d5„ 

Po  »/o  Jo 


which  is  in  a  form  amenable  to  elementary  integration.  Thus, 


General  F(5):  The  result  (B31)  above  is  but  a  particular  case  of  the  following  result  for  general  V(R), 
On  interchanging  the  order  of  integrations  in  (B21)  with  (Bl)  then 

1  pPo  pRo  9  rRt>  fPm(R)  [n's  1 

C(E-,R0)  =  -y  dp2  f  ds,  =  -y  dR  -j-  dp\  (532) 

Po  </o  «/i2,  Po  Jo  Jo  dR 

The  trajectories  which  access  R  have  impact  parameters  p  in  the  range  0  <  p  <  pm(R)  where, 


31 


Since 


AW-4-zm} 


(£33) 


ds{  —  Pm{R) 

Jr\  “ 


from  (B2),  the  inner  integral  of  (B32)  is  simply  2 so  that, 


£(£;£„)  =  4/ pl(R)dR, 

Po  Jo 

to  give 


(534) 


(B35) 


C(E-,R0) 


1*2(1  -  V(Ro)/E) 


]  />(-¥)'* 


(J536) 


as  the  /^averaged  path  length  for  general  V(R),  This  valuable  expression  yields  (B31)  for  Coulombic 
attraction  is  applicable  to  all  V(R)  and  agrees  with  equation  (4.9)  of  text. 


References 

1.  M.  R.  Flannery  and  E.  J.  Mansky,  J.  Chem.  Phys.  88  4228  (1988). 

2.  M.  R.  Flannery,  J.  Chem.  Phys.  89  214  (1988). 

3.  M,  R.  Flannery,  J.  Chem.  Phys.  87  6847  (1987). 

4.  M.,  R.  Flannery  and  E.  J.  Mansky,  J.  Chem.  Phys.  80  4086  (1988). 

5.  D.  R.  Bates  and  C.  S.  McKibbin,  Proc.  Roy.  Soc.  (Lond.)  A  339  13 
(1974). 

6.  M.  R.  Flannery,  Phil.  Trans.  Roy.  Soc.  (Lond.)  A  304  447  (1982). 

7.  W.  Naumann,  Physica  A  150  627  (1988). 

8.  M.  R.  Flannery,  J.  Phys.  B:  At.  Mol.  Phys.  20  4929  (1987). 

9.  M.  R.  Flannery,  J.  Phys.  B:  At.  Mol.  Phys.  18  L839  (1985). 

10.  D.  R.  Bates  and  W.  L.  Morgan,  Phys.  Rev.  Lett.  04  2258  (1990). 

11.  S.  P.  Mezyk,  R.  Cooper  and  J.  Sherwell,  J.  Phys.  Chem.  93  8187 
(1989), 

12.  R.  J.  DiPerna  and  P.  L.  Lions,  Annals  of  Math.  130  321  (1989). 

13.  J,  J.  Thomson,  Phil.  Mag.  47  337  (1924). 

14.  L.  B.  Loeb,  Basic  Processes  of  Gaseous  Electronics)  The  Univ.  Cal. 

Press  (Berkeley)  1955,  Chap.  5 

15.  G.  L.  Natanson,  Sov.  Phys.-  Tech.  Phys.  4  1263  (1959). 

16.  D.  R.  Bates,  J.  Phys.  B:  At.  Mol.  Phys.  13  2587  (1980). 

17.  L.  P.  Pitaevskii,  Sov.  Phys.  -  JETP  15  919  (1962). 

18.  S.  B.  Byron,  R.  C.  Stabler  and  P.  I.  Bortz,  Phys.  Rev.  Lett.  8  376  (1962). 

19.  D.  R.  Bates  and  I.  Mendas,  Proc.  Roy.  Soc.  (Lond.)  A  359  275  (1978). 

20.  D.  R.  Bates  and  M.  R.  Flannery,  J.  Phys.  B:  At,  Mol.  Phys.  2  184  (1969). 

21.  M.  R.  Flannery  and  E.  J.  Mansky,  Chem.  Phys.  132  115  (1989). 

22.  D.  R.  Bates  and  C.  S.  McKibbin,  J.  Phys.  B:  At.  Mol.  Phys.  0  2485  (1973). 

23.  L.  D.  Landau  and  E.  M.  Lifshitz,  Statistical  Phy$ics}  Pergamon  Press  ,  3rd  Edition,  New  York  1980, 
Chap.  3 

24.  M.  R.  Flannery,  in  preparation 


32 


Figure  Captions 


Figure  la-b,  (a)  Probability  Pa(Rq)  E\  A),  equation  (3.27),  for  collision  of  an  ion  in  a  hyperbolic 
(Coulomb)  trajectory  with  the  gas  within  a  sphere  of  radius  Rq  =  0ARe  as  a  function  of  reduced  gas  density 
X  =  jR0/A  at  various  reduced  Coulomb-orbit  energies  e  =  E/kT\  (b)  Dotted  curve:  E-averaged  probability 
P^iZolA),  equation  (3.28),  and  PA(R0,kT;  A).  Solid  lines  are  parabolic  (e  =  0)  and  rectilinear  (e  — >  oo) 
envelopes  to  the  collision  probabilities,  equations  (3.34)  and  (3.31),  respectively.  Inset  figure:  E-variation 
between  envelopes  at  low  gas  density.  The  ion-pairs  here  are  not  in  .^-equilibrium. 

Figure  2.  Collision  Probability  PA(R0,E  =  kT\  A)  as  a  function  of  radius  Ro  for:  (a)  a  fixed  gas 
density  with  A  =  Rri  and  (b)  for  various  gas  densities  proportional  to  /.  Inset  figure:  P  initially  increases 
linearly  with  Ro  along  parabolic  envelope  and  with  further  increase  of  Rq  tends  ultimately  to  the  straight-line 
envelope. 

Figure  3*  The  recombination  rate  ay/(|fl\R;f~)  versus:  (a)  transition-state  energy  £o  =  E/kT)  or  (b) 
reaction  radius  r0  =  Ro/Re  =  lAo- 

Figure  4.  Analytic  Probabilities  PA(Ro,E)i  equation  (6.23)  for  Coulomb  Association  within  a  sphere 
of  radius  Ro  =  0ARe  as  a  function  of  reduced  gas  density  A"  =  i?o/A.  Parabolic  (equation  (6.26))  and 
rectilinear  (equation  (6.27))  limits  e  =  0  and  t  =  oo,  respectively.  Note  that  I2-equilibrium  is  assumed. 

Figure  5.  Comparison  of  E-averaged  association  probabilities  as  a  function  of  reduced  gas  density 
X  =  i?o/A  for  ions  with  (lower  curve)  and  without  (upper  curve)  I2-equilibrium. 
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Figure  la-b.  (a)  Probability  PA(Ro,  E\  A),  equation  (3.27),  for  collision  of  an  ion  in  a  hyperbolic 
(Coulomb)  trajectory  with  the  gas  within  a  sphere  of  radius  Rq  =  0ARe  as  a  function  of  reduced  gas  density 
A'  =  Ro/\  at  various  reduced  Coulomb-orbit  energies  e  =  E/kT ;  (b)  Dotted  curve:  E-averaged  probability 
PA(f?ojA),  equation  (3.28),  and  Pa{Rq>  kT;  A).  Solid  lines  are  parabolic  (s  =  0)  and  rectilinear  (e  —  oo) 
envelopes  to  the  collision  probabilities,  equations  (3.34)  and  (3.31),  respectively.  Inset  figure:  E- variation 
between  envelopes  at  low  gas  density.  The  ion-pairs  here  are  not  in  ^-equilibrium. 


Figure  la-b*  (a)  Probability  Pa(Rq<  E\  A),  equation  (3.27),  for  collision  of  an  ion  in  a  hyperbolic 
(Coulomb)  trajectory  with  the  gas  within  a  sphere  of  radius  Rq  =  0.4iZf  as  a  function  of  reduced  gas  densit\ 
A*  =  Ro/X  at  various  reduced  Coulomb-orbit  energies  e  =  E/kT\  (b)  Dotted  curve:  E-averaged  probability 
PA(Ro\X),  equation  (3.28),  and  PA(RoykT'1  X).  Solid  lines  are  parabolic  (s  =  0)  and  rectilinear  (s  — *  oc) 
envelopes  to  the  collision  probabilities,  equations  (3.34)  and  (3.31),  respectively.  Inset  figure:  E-variation 
between  envelopes  at  low  gas  density.  The  ion-pairs  here  are  not  in  ^-equilibrium 


COLLISION  PROBABILITY  P(X,f,£=l) 


Figure  2.  Collision  Probability  Pa(Rq,  E  =  kT;  A)  as  a  function  of  radius  Ro  for:  (a)  a  fixed  gas 
density  with  A  =  Rt,  and  (b)  for  various  gas  densities  proportional  to  /.  Inset  figure:  P  initially  increases 
linearly  with  Ro  along  parabolic  envelope  and  with  further  increase  of  Ro  tends  ultimately  to  the  straight-line 
envelope. 


COLLISION  PROBABILITY 
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Figure  2.  Collision  Probability  PA(R0,E  =  kT:  A)  as  a  function  of  radius  Ro  for:  (a)  a  fixed  gas 
density  with  A  =  Ret  and  (b)  for  various  gas  densities  proportioned  to  /,  Inset  figure:  P  initially  increases 
linearly  with  Ro  along  parabolic  envelope  and  with  further  increase  of  Ro  tends  ultimately  to  the  straight-line 
envelope. 


TRANSITION-STATE  ENERGY  (E/kT)  REACTION  RADIUS  (R0/Re) 


ASSOCIATION  PROBABILITY  P(X,£) 


Figure  4.  Analytical  Probabilities  Pa(Rq.  E),  equation  (6.23)  for  Coulomb  Association  within  a  sphere 
of  radius  R0  =  0.4R,  as  a  function  of  reduced  gas  density  X  =  *o/A  for  various  e.  Parabolic  (equation 
(6.26))  and  rectilinear  (equation  (6.27))  limits  c  =  0  and  s  =  oo,  respectively.  Note  that  L  -equilibrium  is 
assumed. 
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_Flgure  5;  Comparison  of  E-averaged  association  probabilities  as  a  function  of  reduced  gas  density 
A  —  Rq/  A  for  ions  with  (lower  curve)  and  without  (upper  curve)  Z^-equilibrium. 


